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I.  STATEMENT  OF  THE  PROBLEM  STUDIED 


In  simple  theories  of  penetration,  the  material  properties  of  target  and 
penetrator  are  often  represented  only  by  constant  characteristic  stresses,  as 
for  example  in  Tate  [1],  Although  this  approach  leads  to  results  that  are 
qualitatively  correct,  it  can  be  difficult  to  use  quantitatively.  Some  prob¬ 
lems  have  to  do  with  actual  deformations  in  target  and  penetrator  including 
lateral  motion,  and  others  are  associated  with  the  fact  that  the  plastic  flow 
stress  is  determined  only  by  the  deviatoric  components  of  stress  whereas  the 
spherical  or  pressure  component,  which  may  be  quite  large  ahead  of  the  pene¬ 
trator  and  contributes  significantly  to  the  retardation  of  the  penetrator,  is 
unrelated  to  flow  stress  (e.g.  see  Wright  [2]).  Wright  [2],  in  his  survey 
article  on  long  rod  penetrators,  ellucidated  vividly  some  of  the  problems  with 
the  existing  penetration  models.  In  another  extensive  review  article,  Backman 
and  Goldsmith  [3]  discussed  superbly  the  work  done  in  penetration  mechanics 
until  1977.  Jonas  and  Zukas  [4]  reviewed  various  analytical  methods  for  the 
study  of  kinetic  energy  projectile-armor  interaction  at  ordance  velocities  and 
placed  particular  emphasis  on  three-dimensional  numerical  simulation  of  perfo¬ 
ration.  Anderson  and  Bodner  [5]  have  recently  reviewed  the  status  of  the 
ballistic  impact  modeling.  A  penetration  model  that  is  not  too  difficult  to 
use  has  been  proposed  by  Ravid  and  Bodner  [6],  They  studied  the  penetration 
problem  by  presuming  a  kinematically  admissible  flow  field  in  the  target  and 
found  the  unknown  parameters  by  utilizing  an  upper  bound  theorem  of  plasticity 
modified  to  include  dynamics  effects. 

In  an  attempt  to  shed  some  light  on  questions  raised  by  Wright  [2],  Batra 
and  Wright  [7]  recently  studied  an  idealized  penetration  problem  that  simu¬ 
lates  the  following  situation.  Suppose  that  the  penetrator  is  in  the  interme¬ 
diate  stages  of  penetration  so  that  the  active  target/penetrator  interface  is 
at  least  one  or  two  penetrator  diameters  away  from  either  target  face,  and  the 


remaining  penetrator  is  much  longer  than  several  diameters  and  is  still  trav¬ 
eling  at  a  speed  close  to  its  striking  velocity.  This  situation  has  been 
idealized  as  follows.  It  is  assumed  that  the  rod  is  semi -  inf inite  in  length, 
the  target  is  infinite  with  a  semi- infinite  hole,  the  rate  of  penetration  and 
all  flow  fields  are  steady  as  seen  from  the  nose  of  the  penetrator,  and  that 
no  shear  stress  can  be  transmitted  across  the  target/penetrator  interface. 

This  last  assumption  is  justified  on  the  grounds  that  a  thin  layer  of  material 
at  the  interface  either  melts  or  is  severely  degraded  by  adiabatic  shear. 

These  idealizations  make  it  possible  to  decompose  the  penetration  problem  into 
two  parts  in  which  either  a  rigid  rod  penetrates  a  deformable  target  or  a 
deformable  rod  is  upset  at  the  bottom  of  a  hole  in  a  rigid  target.  Of  course, 
in  the  combined  case  the  contour  of  the  hole  is  unknown,  but  if  it  can  be 
chosen  so  that  normal  tractions  match  in  the  two  cases  along  the  entire  bound¬ 
ary  between  penetrator  and  target,  then  the  complete  solution  is  known  irres¬ 
pective  of  the  relative  motion  at  the  boundary.  Even  without  matching  the 
normal  tractions,  it  would  seem  that  valuable  qualitative  information  about 
the  flow  field  and  distribution  of  stresses  can  be  gained  if  the  chosen  con¬ 
tour  is  reasonably  close  to  that  found  in  experiments. 

We  have  continued  the  work  initiated  by  Batra  and  Wright  [7],  and  have 
studied  the  axisyrametric  thermomechanical  steady  state  deformations  of  a 
viscoplastic  target  being  penetrated  by  a  rigid  cylindrical  penetrator,  and  of 
a  viscoplastic  rod  striking  a  known  hemispherical  rigid  cavity.  The  effects 
of  heat  generated  due  to  plastic  deformations,  compressibility  of  the  mate¬ 
rial,  work-hardening,  strain-rate  hardening,  thermal  softening,  and  the  shape 
of  the  penetrator  nose  in  the  former  case  and  the  cavity  shape  for  the  latter 
problem  have  been  investigated. 

Another  important  and  still  totally  unresolved  problem  is  that  of  select¬ 
ing  the  most  appropriate  constitutive  relation  for  the  material  of  the  body. 
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The  constitutive  relation  employed  should  adequately  model  the  material 
response  over  the  range  of  deformations  expected  to  occur  in  the  problem. 
However,  the  computed  values  of  the  deformation  fields  generally  depend 
strongly  upon  the  constitutive  assumptions  made.  A  way  out  of  this  dilemma  is 
to  choose  a  constitutive  relation,  solve  the  problem,  check  if  the  constitu¬ 
tive  assumptions  are  valid  over  the  range  of  computed  deformations,  and,  if 
necessary,  resolve  the  problem  with  the  modified  constitutive  relation. 

In  the  last  few  years,  many  new  theories  [8-10]  of  large  deformation 
elasto-plasticity  have  been  proposed.  These  theories  make  different  kinematic 
assumptions  thus  necessitating  the  hypothesizing  of  constitutive  relations  for 
variables  which  may  not  be  simply  related  with  each  other.  In  an  attempt  to 
determine  which  of  these  theories  is  the  most  appropriate  for  the  analysis  of 
penetration  problems;  we  find  the  histories  of  the  effective  stress,  second 
invariant  of  the  strain- rate  tensor,  the  temperature  and  the  spin  at  a  visco¬ 
plastic  target  particle  being  penetrated  by  a  long  rigid  cylindrical  penetra- 
tor 

An  important  failure  mechanism  in  penetration  mechanics  is  the  develop¬ 
ment  of  adiabatic  shear  bands  at  or  near  the  target/penetrator  interface.  In 
order  to  increase  our  understanding  of  the  factors  that  enhance/retard  the 
initiation  and  growth  of  adiabatic  shear  bands,  we  have  analyzed  the  dynamic 
deformations  of  a  viscoplastic  block  undergoing  simple  shearing  deformations. 
The  material  inhomogeneity  or  defect  has  been  modelled  by  assuming  an  initial 
non-uniform  temperature  distribution  within  the  block. 

Results  for  the  aforementioned  problems  are  summarized  below. 

II.  SUMMARY  OF  THE  MOST  IMPORTANT  RESULTS 

From  the  results  completed  under  this  contract,  the  following  salient 
conclusions  can  be  drawn. 
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1)  The  peak  hydrostatic  pressure  at  a  material  particle  is  nearly  12 
times  its  yield  stress  in  a  quasistatic  simple  compression  test. 

2)  The  plastic  spin  at  a  material  particle  is  of  the  same  order  of  magni¬ 
tude  as  the  plastic  strain-rate  at  that  particle. 

3)  Severest  target  deformations  occur  near  the  periphery  of  the  blunt 
nose.  The  peak  strain-rates  there  and  in  the  target  material  ahead 
of  the  penetrator  nose  periphery  are  approximately  an  order  of  magni¬ 
tude  higher  than  that  at  other  points.  Also,  for  all  nose  shapes, 
strain-rates  at  target  particles  abutting  the  penetrator  nose  are  very 
high  and  drop  off  rapidly  within  a  short  distance.  The  exact  thick¬ 
ness  of  this  layer  has  not  yet  been  ascertained. 

4)  For  a  thermoviscoplastic  rod  striking  a  rigid  cavity,  peak  strain- 
rates  occur  within  the  rod  material  adjacent  to  the  cavity  wall,  and 
near  the  bottom-most  point  of  the  free  surface. 

5)  The  maximum  temperature  computed  in  the  two  problems  equalled  nearly 
one-half  of  the  presumed  melting  temperature  for  the  deformable  mate¬ 
rial  . 

6)  The  presence  of  multiple  initial  imperfections  delays  considerably  the 
initiation  of  shear  bands  as  compared  to  the  case  when  there  is  only 
one  initial  perturbation. 

7)  The  consideration  of  dipolar  effects  delays  the  initiation  of  shear 
bands,  and  inertia  forces  start  playing  a  noticeable  role  at  an 
applied  overall  strain-rate  of  5,000  s '  ^ . 

III.  BRIEF  REVIEW  OF  THE  COMPLETED  WORK 

In  their  analysis  of  the  steady  state  axisymmetric  deformations  of  a 
rigid  perfectly/plastic  target  Batra  and  Wright  [7]  found  that  the  nondimen- 
sional  resisting  force  F  acting  on  the  hemispherical -nosed  rigid  cylindrical 
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penetrator  is  given  by 

F  -  3.903  +  0.0773  a,  l  <  a  <  5,  (1) 

where  a  -  pv2/a0  is  a  nondimens ional  number,  p  is  the  mass  density  of  the 
o 

target  material,  a0  its  yield  stress,  and  vQ  is  the  striking  speed  of  the 
penetrator.  Wright  [11]  recently  used  this  result  to  estimate  one  of  the 
terms  in  his  re-examination  of  Tate's  theory.  Wright  [2]  and  Pidsley  [12] 
have  pointed  out  that  if  the  equation  of  motion  for  steady  flow  is  integrated 
along  the  central  streamline,  there  is  a  contribution  from  transverse  gra¬ 
dients  of  shear  stress,  unlike  the  case  for  a  perfect  fluid.  Figure  1,  taken 
from  Batra  and  Wright's  paper,  shows  the  contributions  from  various  terms  and 
clearly  establishes  that  the  transverse  gradients  of  shear  stress  make  a 
significant  contribution  to  the  total  and  that  the  deviatoric  component  stays 
constant  at  approximately  0.75  (compared  to  the  theoretical  exact  value  of 
0.667)  out  of  a  total  of  8.5. 

Batra  [13]  studied  the  effect  of  the  nose  shape  and  various  material 
parameters  on  the  steady-state  penetration  of  viscoplastic  targets  by  rigid 
cylindrical  penetrators.  The  constitutive  relation  for  the  viscoplastic 
material  was  taken  to  be 

a  -  -p  1  +  <r0(l  +  bI)mD/y3I,  21  -  trD2  ,  (2) 


H  i 


»>v  * 


••  *•  .. 

wo*  rv*  now 

Figure  1.  Contribution  of  Various  Terms  in  the  Bernoulli 
Equation. 
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where  D  is  the  stretching  tensor,  aQ  the  yield  stress  in  simple  compression, 
a  t!  _  Cauchy  stress  tensor,  p  the  hydrostatic  pressure  not  determined  by  the 
deformation  history  of  the  target,  and  b  and  m  are  material  parameters. 

Figure  2  shows  that  the  change  in  nose  shape  from  hemispherical  to  ellipsoidal 
reverses  the  curvature  of  the  normal  stress  versus  the  angular  position  curve 
from  convex  downwards  to  concave  downwards.  Also,  as  is  evident  from  Figure 
3,  the  axial  resisting  force  depends  strongly  upon  the  nose  shape  and  rather 
weakly  upon  m  and  a. 

When  the  effects  of  work-hardening  of  the  target  material  are  also  con¬ 
sidered,  and  the  constitutive  relation  (2)  is  changed  to 

a  -  -p  1  +  aj  (1  +  bl)m(l  +  — )n  D/,/31  (3) 

where  0O  and  n  are  material  parameters,  and  is  an  internal  variable 

whose  growth  depends  upon  the  plastic  working,  Batra  [14]  noticed  that  the 
results  are  very  sensitive  to  the  finite  element  mesh  used.  This  is  due  to 
the  absence  of  any  diffusive  term  in  the  constitutive  relation  for  0. 
Subsequently,  Batra  [13]  accounted  for  material  softening  caused  by  the  tem¬ 
perature  rise  due  to  plastic  working  and  modified  the  constitutive  relation  to 

t!> 

a  -  -  p  1  +  ct0  '1  +  bl )m( 1  +  — )n(l  -  7 9)  D//3I  (4) 

*o 


Figure  2.  Distribution  of  the  normal  pressure  on  the  pene- 
trator  nose  for  three  different  nose  shapes. 


_  Blunt  nose  (rn/r0  -  0.2);  hemis¬ 
pherical  nose  (rn/rQ  -  1.0);  .  ellipsoidal 

nose  (rn/rQ  -  2.0). 
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(ALFA) 


Figure  3.  Dependence  of  the  axial  force  upon  various  fac¬ 
tors  . 

where  y  is  the  thermal  softening  coefficient,  and  8  the  temperature  rise. 

For  values  of  material  parameters  that  correspond  to  a  typical  steel,  Batra 
[15]  computed  the  strain- rate  distribution  depicted  in  Figure  4.  This  agrees 
qualitatively  with  the  experimental  observation  of  Wingrove  [16]  that  a  shear 
band  forms  and  thus  very  high  rates  occur  in  the  target  material  near  the 
periphery  of  a  flat-ended  projectile  (cf.  Fig.  5).  Figure  6,  taken  from  Lin 
and  Batra' s  [17]  paper,  shows  the  stress  and  strain-rate  histories  at  three 
material  particles  that  were  near  the  centroidal  axis  in  the  undeformed  target 
being  penetrated  by  a  blunt-nosed  penetrator.  The  time  is  clocked  from  the 


Figure  4.  Distribution  of  the  strain-rate  invariant  I  near 
the  periphery  of  a  blunt-nosed  penetrator. 


Figure  5.  Shear  band  produced  below  a  flat-ended  projectil 
showing  the  extent  of  shear  in  the  band.  Projec 
tile  velocity  140  m/s ,  flat-nosed  penetrator. 
(Wingrove  [16]) 


Figure  6.  Histories  of  stress  and  strain-rate  at  three 

material  particles  that  were  initially  near  the 
centroidal  axis. 

instant  these  particles  are  on  the  boundary  of  the  spatial  domain  ana¬ 
lyzed.  Clearly  the  effective  stress  at  each  of  these  particles  attains  a 
maximum  and  then  decreases  slowly  due  to  the  softening  caused  by  the  tempera¬ 
ture  rise.  This  drop  in  the  stress  is  more  pronounced  for  a  blunt  nosed 
penetrator  as  compared  to  that  for  a  hemispherical  nosed  one.  Since  a  shear 
band  usually  initiates  at  a  point  after  the  stress  there  has  reached  the  peak 
value,  there  is  more  likelihood  of  a  band  being  formed  for  a  blunt-nosed 
penetrator  as  compared  to  that  for  a  hemispherical  one. 

Batra  and  Gobinath  [18]  showed  that  the  compressibility  of  the  target 
material  reduces  the  axial  resisting  force  experienced  by  the  penetrator  by 
f  nearly  10  percent,  and  the  peak  values  of  the  strain- rate  invariant  I  and  the 

l  normal  sttess,  both  of  which  occur  at  the  target  particle  situated  on  the 
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penetrator  nose- tip,  by  8 . 8  percent  and  3.2  percent,  respectively. 

Batra  and  Lin  [19]  have  analyzed  the  steady  state  axisymmetric  deforma¬ 
tions  of  a  rod  striking  a  rigid  cavity.  This  problem  is  more  challenging  than 
the  one  involving  a  deformable  target  and  rigid  penetrator  because  of  the 
presence  in  it  of  a  free  surface  whose  shape  is  not  known  a  priori .  When  the 
rod  material  is  assumed  to  exhibit  strain  and  strain- rate  hardening  and  ther¬ 
mal  softening,  the  axial  force  experienced  by  the  rod  was  found  to  depend 
strongly  upon  a,  the  peak  value  of  the  nondimens ional  strain- rate  invariant 
was  found  to  depend  weakly  upon  a,  and  the  maximum  temperature  equal  to  nearly 
one -third  the  presumed  melting  temperature  of  the  rod  was  found  to  occur.  All 
of  the  material  particles  abutting  the  cavity  wall  had  the  same  temperature, 
as  shown  in  Figure  7..  Figure  8  depicts  the  history  of  the  effective  stress  at 
four  material  particles,  two  of  which  initially  are  located  near  the  centroidal 


Figure  7.  Temperature  distribution  in  a  viscoplastic  rod 
upset  at  the  bottom  of  a  hemispherical  cavity. 
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Figure  8.  Stress  histories  at  four  material  particles,  two 
of  these  particles  are  initially  near  the  cen- 
troidal  axis  and  the  other  two  are  near  the  free 
surface . 

axis  and  the  other  two  near  the  free  surface  of  the  undeformed  rod.  The 
stress  history  experienced  by  the  material  particles  initially  near  the  cen- 
troidal  axis  suggests  that  the  likelihood  of  a  shear  band  being  formed  near 
the  cavity  wall  is  quite  high. 

In  each  of  the  two  aforementioned  problems ,  the  peak  hydrostatic  pressure 
at  a  material  point  was  found  to  be  nearly  twelve  times  the  flow  stress  of  the 
material  and  the  peak  plastic  spin  to  be  10^-10^  sec'^-.  We  note  that  the 
constitutive  relation  employed  neither  accounted  for  the  dependence  of  the 
flow  stress  upon  the  hydrostatic  pressure  nor  accorded  any  importance  to  the 
plastic  spin.  Also,  no  fracture  or  failure  criterion  was  included  in  the 
model . 

In  problems  aimed  at  developing  an  understanding  of  factors  that 
enhance/retard  the  initiation  and  growth  of  adiabatic  shear  bands  in  armor 
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materials,  we  found  that  the  consideration  of  dipolar  effects  [20]  delays 
significantly  the  initiation  of  adiabatic  shear  bands.  Shear  bands  initiating 
at  neighboring  defects  that  would  grow  independently  in  nonpolar  materials 
seem  to  merge  together  and  develop  as  a  single  band  in  dipolar  materials. 
Inertia  forces  start  playing  a  significant  role  [21]  when  the  applied  strain- 
rate  in  a  simple  shearing  problem  exceeds  5,000  sec"^.  During  the  time  the 
deformation  is  localizing,  the  shear  stress  in  the  region  of  localization 
collapses  and,  if  the  rate  of  collapse  of  the  stress  is  dramatic,  an  unloading 
elastic  wave  emanates  [22]  out  of  this  region. 
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STEADY-STATE  PENETRATION  OF  VISCOPLASTIC  TARGETS 

R  C.  BATRA 

Department  of  Engineering  Mechanics.  University  of  Missouri  Rolla.  Rolla.  M0fe540l-0249.  USA 

Abstract  The  problem  of  steady-state  penetration  by  a  semi-infinite,  rigid  cylindrical  penetrator 
with  an  ellipsoidal  nose  into  an  infinite,  rigid  viscoplastic  target  has  been  studied  The  target 
material  is  assumed  to  obey  a  generalized  form  of  Von  Mises  yield  criterion  to  account  for  the 
strain-rale  dependence  C  ontact  between  target  and  penetrator  is  assumed  to  be  smooth.  Computed 
results  show  that  the  deformation  field  adjacent  to  the  nose  of  the  penetrator  is  significantly 
different  in  the  ellipsoidal  case  from  what  it  is  when  the  nose  is  hemispherical  Results  presented 
graphically  include  the  dependence  of  the  axial  resisting  force  on  penetrator  speed,  the  ratio  of  the 
major  to  minor  axes  of  its  ellipsoidal  nose,  and  the  strain-rate  hardening  parameter  of  the  target 
Also  depicted  are  the  normal  pressure  over  the  penetrator  nose  and  the  velocity  held  in  different 
parts  ol  the  target 


INTRODUCTION 

In  simple  theories  of  penetration,  the  material  properties  of  target  and  penetrator  are  often 
represented  only  by  constant  characteristic  stresses,  as  for  example  in  Tate  [I]  Although 
this  approach  leads  to  results  that  are  qualitatively  correct,  it  can  be  difficult  to  use 
quantitatively  Some  of  the  problems  have  to  do  with  actual  deformations  in  target  and 
penetrator  including  lateral  motion,  and  others  are  associated  with  the  fact  that  the  plastic 
flow  stress  is  determined  only  bv  the  deviatonc  components  of  stress  whereas  the  spherical 
or  pressure  component,  which  may  be  quite  large  ahead  of  the  penetrator  and  contributes 
significantly  to  the  retardation  of  the  penetrator.  is  unrelated  to  flow  stress  le.g.  see  Wright 
C-]T  In  developing  an  engineering  model  for  penetration  and  perforation.  Ravid  and 
Bodner  [3]  assumed  simple  kinematics  for  the  flow  around  the  penetrator  and  then 
adjusted  the  unknown  parameters  by  utilizing  an  upper  bound  theorem  of  plasticity 
modified  to  include  dynamic  effects. 

In  [4],  Batra  and  Wright  have  presented  a  detailed  numerical  solution  to  the  following 
idealized  penetration  problem  It  was  assumed  that 

(II  the  rod  is  semi-infinite  in  length  and  that  the  target  is  infinite  with  a  semi-infinite  hole. 

(2)  the  rate  of  penetration  and  all  flow  fields  are  steady  as  seen  from  the  nose  of  the 
penetrator. 

(3)  no  shear  stress  can  be  transmitted  across  the  target  penetrator  interface. 

|4>  the  deforming  material  was  taken  to  be  rigid  perfectly  plastic. 

They  studied  the  problem  of  the  deforming  target  and  a  rigid  penetrator  having  a  circular 
cylindrical  body  and  a  hemispherical  nose 

Batra  and  Wright's  calculations  revealed  that  strain  rates  in  the  target  material  that  is 
ahead  of  the  penetrator  arc  of  the  order  of  10'  sec  1  Since  many  materials  used  in  such 
applications  have  strain-rate  sensitive  properties,  we  extend  herein  the  previous  work  to 
viscoplastic  materials  Also,  the  penetrator  nose  is  taken  to  be  ellipsoidal  As  in  the  previous 
work  [4],  the  objective  here  is  to  study  the  idealized  penetration  problem  in  detail  and 
possibly  shed  some  light  on  the  aforementioned  factors.  The  problem  studied  herein 
simulates  approximately  the  following  situation  the  penetrator  is  in  the  intermediate  stages 
of  penetration  so  that  the  active  target  penetrator  interface  is  at  ieast  one  or  two  penetrator 
diameters  away  from  either  target  face,  and  the  remaining  penetrator  is  still  much  longer 
than  several  diameters  and  is  still  traveling  at  a  speed  close  to  its  striking  velocity  Tor 
this  case  the  first  two  assumptions  stated  in  the  second  paragraph  above  arc  quite 
reasonable  and  are  also  made  in  this  work.  It  should  be  emphasized  that  we  have  not 
incorporated  any  fracture  or  failure  criterion  in  our  work  Thus  the  material  is  presumed 
to  undergo  unlimited  plastic  deformations 
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FORMULATION  OF  THE  PROBLEM 

We  presume  that  the  deformations  of  the  target  appear  to  be  independent  of  time  to  an 
observer  situated  on  the  penetrator  nose  and  moving  with  it  at  a  uniform  velocity  r0e.  e 
being  a  unit  vector  along  the  direction  of  motion  of  the  rigid  penetrator.  We  use  a 
cylindrical  co-ordinate  system  attached  to  the  center  of  the  penetrator  nose,  with  ;-axis 
pointing  into  the  target. 

Equations  governing  the  target  deformations  are: 


div  v  =  0. 

ID 

p(v  grad)v  =  div  <r. 

(2) 

Here  v  is  the  velocity  of  a  target  particle  as  seen  by  an  observer  situated  on  the  penetrator. 
P  is  the  mass  density  and  a  is  the  Cauchy  stress  tensor.  Equation  (1)  implies  that  the 
deformations  of  the  target  are  isochoric.  and  eqn  (2)  expresses  the  balance  of  linear 
momentum.  We  neglect  the  elastic  deformations  of  the  target  and  assume  that  its  material 
obeys  the  following  constitutive  relation  for  a 

a  =  —pi  +  2p(/)D.  1  ^  0. 

(3) 

2 /!</)  =  (Tod  +  hir  v  3  /, 

|4| 

D  =  [gradv  •+■  (gradv)r]  2. 

(5) 

=  UrD* 

(6) 

In  these  equations,  p  is  the  hydrostatic  pressure  that  is  not  determined  by  the  deformation 
history  of  the  target.  D  is  the  stretching  tensor.  rr0  is  the  yield  stress  in  simple  tension  or 
compression,  parameters  h  and  m  describe  the  strain-rale  hardening  of  the  material,  and 
trlD^)  equals  the  sum  of  the  diagonal  terms  of  the  square  matrix  D;  Equation  (3)  can 
also  be  viewed  as  a  constitutive  relation  for  an  incompressible  viscous  fluid  with  viscosity 
coelficient  equal  to  n0(  1  +  hl)m  (2X  3  /).  Implicit  in  eqn  (3)  is  the  assumption  that  the  Von 
Mises  yield  surface  is  given  by 


trls;)  =  y  ctrI  I  +  hl)Zm.  (7) 

s  =  a  +  pi.  (8) 

The  tensor  s  is  the  deviatoric  stress  tensor. 

Ravid  and  Bodner  [3]  have  used  a  constitutive  relation  similar  to  eqn  (3)  and  assumed 
that 


2p\l)  =  <t0(1  +  Clogl0(2/  v  3)1  v  3/.  (9) 

where  C  is  a  material  constant  and  is  taken  to  be  zero  for  strain  rates  lower  than  unity 
The  constitutive  relation  for  a  Bingham  soiid  [5]  results  by  taking  m  =  I  and  interpreting 
h  properly  Zienkicwicz  et  ul.  [6]  look 

2pU)  =  [n„  +  (2/  3)1'"]  x  3/  1 10) 

and  asserted  that  it  corresponds  to  Perzvna  s  viscoplastic  model.  In  eqn  (10)  y  and  n  are 
temperature-dependent  material  constants.  Our  choice  [eqns(3.4)]  was  motivated  by  the 
desire  to  generalize  the  power  law  model  used  by  Burns  [7]  and  Shawki  ct  ul.  [8]  so  that 
it  is  also  valid  for  quasistatic  tests.  This  generalized  constitutive  model  that  also  includes 
thermal  softening  and  strain-hardening  has  been  used  by  Wright  and  Batra  [9],  and  Batra 
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[10]  to  study  adiabatic  shear  bands.  For  the  simple  shearing  stress  state,  a  curve  fit  to 
Costin  et  al.'s  [11]  experimental  data  gives  b  =  104  and  m  =  0.025  for  a  hard  steel. 

Equation  (1)  and  equations  obtained  by  substitution  of  eqns  ( 3>— ( 5)  into  eqn  (2)  are  the 
field  equations  which  together  with  suitable  boundary  conditions  are  to  be  solved  for  p 
and  v.  Before  stating  the  boundary  conditions,  we  non-dimensionalize  the  variables  as 
follows. 

6  =  ota0,  s  =  s-<r0.  P  =  P't r0. 

f  -  r  r0.  :  =  :r0.  v  =  vr0.  F  =  bi0.r0.  (11) 

Here  (r,  -)  denote  the  co-ordinates  of  a  point  with  respect  to  the  cylindrical  co-ordinate 
system  chosen,  and  r0  is  the  radius  of  the  cylindrical  body  of  the  penetrator  The  field 
equations  written  in  terms  of  the  non-dimensional  variables  are: 

div  v  =  0  (12) 

-gradp  +  div  |^[(gradv  +  (gradv)T][  =  x(v  grad)v.  (13) 

where 

#i(/)  =  (l  +  Wr  2V  3/.  (14) 

and 

x  =  pii  a o  .  (15) 

is  a  non-dimensional  number.  In  writing  eqns  ( 12)  ( 14).  we  have  dropped  the  superimposed 
bars  and  have  used  grad  and  div  to  denote  the  gradient  and  divergence  operators  in  non- 
dimensional  co-ordinates. 

For  the  boundary  condition  on  the  penetrator-target  interface,  we  assume  that 

t  (<rnl  =  0.  (16) 

v  n  =  0.  (17) 

where  n  and  t  are.  respectively,  a  unit  normal  and  a  unit  tangent  vector  on  the  interface. 
The  boundary  condition  (16)  represents  smooth  contact  between  the  penetrator  and  target. 
This  appears  reaonable  since  a  thin  layer  of  material  at  the  interface  either  melts  or 
is  severely  degraded  by  adiabatic  shear.  The  boundary  condition  (17)  represents  no 
interpenetration  of  the  target  material  into  the  penetrator  and  vice  versa.  At  points  far 
away  from  the  penetrator  we  require  that 

|v  +  e|  -*0  as  (r;  +  r:)'  :  —  r . .  (IK) 

l<rnl  — ►  0  as  r  j  rn.  (19) 

where  e  is  a  unit  vector  along  the  positive  r-axis.  as  before.  That  is.  the  target  material 
ahead  of  the  penetrator  nose  and  far  away  from  it  appears  to  approach  the  penetrator 
with  a  uniform  velocity  and  the  one  behind  the  nose  but  very  far  from  it  is  virtually 
traction  free. 

Note  that  the  field  eqns  (13)  are  nonlinear  in  v.  A  solution  of  eqns  (12)  and  (13)  under 
the  boundary  conditions  (16)  (19),  if  there  exists  one.  will  depend  upon  (he  rate  at  which 
quantities  in  1 18)  and  ( 19)  decay  to  zero  Even  for  the  prescribed  rale  of  decay,  the  solution 
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may  not  be  unique.  We  will  gloss  over  these  rather  difficult  questions  and  seek  an 
approximate  solution  of  these  equations  numerically.  The  hope  is  that  the  numerical 
solution  is  meaningful  for  the  physical  problem  at  hand. 

FINITE  ELEMENT  FORMULATION  OF  THE  PROBLEM 
We  replace  the  infinite  target  region  by  the  bounded  region  R  shown  in  Fig.  I  and  the 
boundary  conditions  (18)  and  (19)  by 

=  0.  rr  =  0  on  tne  surface  AB. 

<rrz  =  0.  v,  =0  on  the  axis  of  symmetry  DE, 

r,  =  0.  iz  =  -1.0  on  the  bounding  surface  EFA.  (20) 

That  the  region  considered  is  adequate  is  justified  by  the  computed  results  presented  below 
which  show  that  noticeable  deformations  of  the  target  occur  only  in  the  region  surrounding 
the  penetrator  and  at  target  particles  whose  distance  form  the  penetrator  is  at  most  2r0. 

Referring  the  reader  to  [12]  for  details,  we  simply  note  that  a  weak  formulation  of  the 
problem  defined  by  eqns  (12).  (13).  (16).  (17)  and  (20)  is  that  equations 


Fig  1  The  region  modeled 
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hold  for  every  smooth  functions  i/t  and  ^  defined  on  R  such  that  4>,  =  0  on  AB  and  DE; 
4>  =  0  on  EFA  and  ^  n  =  0  on  the  target -penetrator  interface  BCD.  Since  eqns  (21)  and 
(22)  are  nonlinear  in  v.  the  following  iterative  technique  has  been  employed: 


i//(div  v‘)dr  =  0, 


(23) 


p‘ldiv^)di  - 


^1‘  ' ) ;  D‘ :  [grad  4>  +  (grad  0)T] 


J  dr  =  i  j  [(v' 
Jr 


"  1  grad)v']  •  <p  dr. 
r  (24) 


where  i  is  the  iteration  number.  The  solution  was  taken  to  have  converged,  if  at  each  nodal 
poim. 

|v‘  —  v‘ '  ll  $0.01  |v* ~ 1 U 


where 


i  1 2  2.2 

(v(-  =  r,  +  v‘. 

COMPUTATION  AND  DISCUSSION  OF  RESULTS 

The  finite  element  code  developed  earlier  [4]  to  solve  the  problem  when  the  target 
material  is  modeled  as  rigid-perfectly  plastic  and  the  penetrator  nose  is  hemispherical  has 
been  modified  to  solve  the  present  problem.  It  employs  six-noded  isoparametric  triangular 
elements  with  rr  and  rz  approximated  by  quadratic  functions  over  an  element  and  p  by  a 
linear  function  defined  in  terms  of  its  values  at  the  vertices  of  the  triangular  element.  The 
accuracy  of  the  code  was  established  by  solving  a  simple  problem  for  an  incompressible 
Navter  Stokes  fluid. 

The  computational  procedure  was  started  by  taking  v°  =  0.  n(l0)  -  104  It  took  19 
iterations  for  the  solution  to  converge  when  x  =  2  and  the  nose  is  hemispherical.  The 
number  of  iterations  required  to  obtain  the  solution  decreased  with  an  increase  in  the  ratio 
r„  r„  when  x  =  2.  However,  for  rn/r0  =  I.  the  number  of  iterations  increased  with  x.  but 
for  rn  r„  =  2  it  decreased  with  an  increase  in  x.  Numerical  experiments  with  different  grids 
were  conducted  first  by  increasing  the  number  of  elements  used  and  then  by  varying  their 
size  but  keeping  the  number  of  elements  used  constant.  The  grid  with  8  rows  of  elements 
in  the  axial  direction.  8  uniformly  spaced  rows  of  elements  in  the  circumferential  direction 
and  8  rows  behind  the  plane  :  =  0.  was  found  to  be  optimum  in  the  sense  that  the  change 
in  the  normal  stress  at  the  penetrator  nose  tip  was  less  than  0.1%  with  further  refinements 
of  the  grid.  The  grid  used  had  a  pattern  similar  to  that  employed  in  [4],  The  dividing 
surfaces  between  elements  intersected  the  r-axis  at  points  distant  0.475.  0.969.  1.483.  2.017. 
2.573.  3.151.  3.752.  4.377  from  the  nose  tip.  The  :-co-ordinate  of  the  horizontal  surfaces 
between  rows  of  elements  behind  the  r  =  0  plane  were  -0.156.  —0.335.  —0.540.  —0.773. 
-  1.041.  -  1.347.  -  1.697.  -2.09.  For  x  =  2.0.  rn,r0  =  2.0.  m  =  0.0.  and  h  =  0.0.  the  total 
axial  force,  defined  below  by  eqn  (29),  obtained  by  using  4.  6.  7.  12  and  13  quadrature 
points  was  computed  to  be  2.6117.  2.6479.  2.6524.  2.6714  and  2.6720  respectively.  The 
results  presented  below  have  been  obtained  by  using  4  quadrature  points. 

Figure  2  depicts  the  velocity  of  the  target  particles  for  x  =  4  and  r„  r0  =  2.0  as  seen  by 
an  observer  sitting  on  the  penetrator  nose.  It  is  apparent  that  significant  deformations  of 
the  target  occur  at  points  within  2r0  of  the  penetrator  boundaries.  The  velocity  fields  for 
other  cases  studied  have  a  similar  pattern.  At  target  points  that  lie  to  the  rear  of  the  center 
of  the  penetrator  nose,  the  flow  quickly  becomes  parallel  to  the  axis  of  the  penetrator.  In 
Fig.  3  is  shown  the  non-dimensional  normal  stress  on  the  penetrator  nose  for  different 
choices  of  various  parameters.  Values  of  parameters  for  various  curves  in  this  and 
subsequent  figures,  unless  specified  otherwise,  are  identified  in  Table  1  The  change  in  nose 
shape  from  hemispherical  to  ellipsoidal  reverses  the  curvature  of  the  normal  stress  vs  0 
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Table  1.  Legend  for  figures 
Curve  type  Parameter  values 

-  -  2.0.  b  -  10*.  m  =  0.03.  i  •=  2.0 

rv 

-  =  2.0.  b  =  I04.  m  =  0.03.  i  *  2.0 

ro 

-  =  20.  h  =  0.  m  -  0.  I  «  2.0 

r0 

-  =  2.0.  b  =  0.  m  =  0.  i  =  4  0 
r„ 

-  =  I  0.  b  =  0.  m  =  0.  J  -  2  0 
r„ 

-  =  0.2.  />  =  0.  m  =  0.  j  =  2.0 
rl( 


curve  from  convex  downwards  to  concave  downwards.  This  change  was  observed  even  for 
r„  r0  =  1.1.  For  x  =  2,  the  normal  stress  at  the  nose  tip  was  essentially  unchanged  for 
0  2  $  r„  r0  $  2.4.  For  rn,r0  =  0.2,  the  nose  shape  is  essentially  fiat  and,  as  expected,  the 
normal  stress  on  it  stays  uniform  and  drops  off  rapidly  near  the  periphery.  The  increase 
in  x.  m  or  h  increases  the  normal  stress  on  the  penetrator  nose. 

The  non-dimensional  axial  resisting  force  expenenced  by  the  penetrator  as  various 
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Fig.  4.  Dependence  of  axial  force  on  various  factors. 

parameters  are  changed  is  plotted  in  Fig.  4.  Equations  obtained  by  fitting  curves  to  the 
computed  data  are: 


F  =  2.58(  1 

+  0.00738a  -  0.000569a2) 

1  x  $  5. 

^  =  2.0,  b  =  m  = 
fo 

^  0. 

(25) 

h  =  2.61  1(1 

+  6.761m  +  20.465m2), 

0  $  m  <  0.03, 

-  =  2.0.  b  =  104. 

r  o 

w 

II 

1  -J 

p 

(26) 

F  =  4.456<  1 

+  7.28m  +  23.478m2). 

0^ 

m  ^  0.03, 

© 

II 

-c 

o' 

II 

c|  o 

x  =  2.0. 

(27) 

r 

/  x2 

/ 

'  V  (r  Y~ 

F  =  1 3.041 

1.0  -  1.3086  —  +  0.9504 

(  —  ) 

-  0.3432 

^  +  0.0477  - 

ro 

Vo/ 

\ 

To/  Vo)  J 

0.2  $  -  ^  2.4.  6  =  m  =  0.  x  =  4  (28) 

fo 
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F  =  10.857 


1.0 


.0455  ~  +  0.6294(  -A  -  0.1975(  -  j  +  0.0246( 

ro  V  ro/  \r  o  7  \  r  o/ 


0.2  ^  ^  2.4.  h  =  m  =  0.  *  = 

r,, 


where 


-=r- 


<rn)  sin  Ocos  </> 


sin;  0  +  (  —  |  cos;  0 

\r()J 


dO. 


(29) 


<t>  is  the  angle  which  the  unit  normal  to  the  nose  makes  with  the  penetrator  axis  as  shown 
in  Fig.  I.  The  corresponding  axiai  force  in  physical  units  is  given  by  Finr^a,,).  It  is  clear 
from  these  results  that  changing  the  nose  shape  from  hemispherical  to  ellipsoidal  with 
rn  r„  =  2.0  reduces  the  axial  resisting  force  to  one-half  of  its  value.  For  a  fixed  rn  r„  and 
m,  the  axial  force  depends  rather  weakly  upon  x.  This  weak  dependence  of  F  upon  x  is 
perhaps  one  reason  why  the  choice  of  constant  target  resistance  in  simple  theory  of  Tate 
[I]  gives  good  qualitative  results.  Even  though  the  coefficient  of  m 2  in  eqns  (25)  and  (26) 
is  nearly  3  times  that  of  m.  the  dependence  of  the  axial  force  F  upon  the  strain-rate 
hardening  coefficient  m  is  essentially  linear.  It  is  so  because  the  values  of  m  for  typical 
steels  arc  much  smaller  than  one.  This  linear  dependence  of  F  upon  m  becomes  transparent 
from  the  plots  of  Fig.  4. 

Computed  results  indicate  that  for  I  $  x  $  5.  rn  r„  =  2.0.  h  =  m  =  0.  the  hydrostatic 
pressure  at  the  nose  tip  increases  essentially  linearly  from  5.6  to  8.2  and  the  normal 
compressive  stress  from  6.22  to  8.84  It  is  apparent  that  the  hydrostatic  pressure  contributes 
significantly  to  the  values  of  er„  and  <r„.  Figure  5  shows  that  the  principal  stress  component 
-<r„  along  the  axis  in  front  of  the  penetrator  falls  off  rapidly  with  the  distance.  The 
stresses  at  points  situated  more  than  3 rn  from  the  nose  tip  cannot  be  accurately  determined 
since  the  velocity  gradients  there  are  extremely  small  The  values  of  rr„  at  points  on  the 
centerline  which  are  at  a  distance  of  3r„  and  more  from  the  nose  tip  equal  essentially  the 
hydrostatic  pressure. 

In  Fig.  5  is  also  plotted  the  variation  of  the  strain-rate  measure  /  along  the  axial  line. 
These  nondimensiona!  values  need  to  be  multiplied  by  r0,ro.  which  typically  equals  I  O', 
to  arrive  at  the  corresponding  dimensional  effective  strain-rate  measure  /.  Thus  strain- 
rates  of  the  order  of  2  x  10' sec  '  occur  at  points  in  the  vicinity  of  the  penetrator  nose. 
On  the  axial  line,  peak  values  of  /  are  higher  near  the  nose  tip  for  r„.r0  =  2.0  as  compared 
with  those  for  lower  values  of  r„  rn.  However,  near  the  periphery  of  the  penetrator  nose 
and  its  cylindrical  body,  the  strain-rate  measure  /  for  the  blunt  nose  increases  by  an  order 
of  magnitude  whereas  that  for  the  ellipsoidal  and  hemispherical  nose,  it  drops  by  a  factor 
of  10.  This  is  shown  in  Fig.  6  wherein  is  also  plotted  the  nondimensional  velocity,  tangential 
to  the  penetrator  nose,  for  different  values  of  various  parameters.  Again  the  nose  shape  is 
the  major  influencing  factor.  However.  rz  at  points  on  the  centerline  ahead  of  the  penetrator 
nose  is  not  affected  that  much  by  the  nose  shape  except  for  r„  r„  =  0.2.  The  decay  rate  of 

|r,  +  e|  for  I  $  —  ^  2.4  is  nearly  constant  but  is  appreciably  less  for  rn.r„  =  0.2.  The  nose 

ro 

shape  does  change  rather  noticeably  the  relative  r-velocily  of  points  on  the  lines  r  =  0  and 
r  =  2.09.  This  follows  from  the  results  depicted  in  Figs  7  and  8  which  also  show  that  more 
of  the  target  material  at  the  sides  of  the  penetrator  deforms  for  a  penetrator  with  an 
ellipsoidal  nose,  even  though  it  is  not  true  ahead  of  the  penetrator.  as  noted  earlier.  Two 
of  the  curves  in  these  figures  essentially  overlap. 

We  note  that  the  target  material  within  some  distance  inside  the. boundary  EFA  does 
not  deform  at  all  and  the  deformations  are  essentially  independent  of  r  near  the  boundary 
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Fig  b  Variation  of  tangential  velocity  and  strain-rate  measure  on  the  penelrator  nose 
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Fig  i  Variation  of  relative  r-velocity  with  radial  co-ordinate  r  on  the  line  :  =  0 
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AB  It  is  therefore  reasonable  to  conclude  that  the  target  region  chosen  for  computation 
is  sulhcient  to  obtain  a  good  description  of  the  deformations  in  the  vicinity  of  the  penetrator 
nose. 


CONCLUSIONS 

For  the  range  of  values  of  r„,rn.  x  and  m  studied,  the  ratio  r„.r0  of  the  radius  of  the 
nose  to  that  of  the  cylindrical  body  of  the  penetrator  has  the  most  effect  on  the  axial 
resisting  force  experienced  by  the  penetrator.  The  retardation  force  depends  rather  weakly 
on  x  and  m.  In  all  cases,  significant  target  deformations  occur  only  within  2 r0  of  the 
penetrator.  Whereas  for  rn,r0  =  1.0.  the  target  material  adjacent  to  the  sides  of  the 
penetrator  appears  to  extrude  rearwards  in  a  uniform  block  that  is  separated  from  the 
bulk  of  the  target  material  by  a  narrow  region  with  a  sharp  velocity  gradient,  such  is  not 
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the  case  for  rn  r0  =  2.0.  In  the  latter  case  there  is  no  noticeable  region  of  steep  velocity 
gradients. 

For  x  =  5  and  rn  r0  =  2.0  the  target  material  seemed  to  separate  from  the  sides  of  the 
penetrator  nose.  This  provides  a  limiting  value  for  the  validity  of  our  calculations  since 
flow  separation  has  not  been  built  into  our  model.  Maximum  normal  stress  and  the 
maximum  hydrostatic  pressure  occur  at  the  nose  tip.  and  fall  off  rapidly  away  from  that 
point.  These  peak  values  depend  significantly  upon  %. 

We  hope  that  results  presented  here  will  prove  useful  in  devising  or  checking  the  results 
from  simpler  engineering  theories  of  penetration. 
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STEADY  STATE  DEFORMATIONS  OF  A  RIGID 
PERFECTLY  PLASTIC  ROD  STRIKING  A  RIGID  CAVITY 

R.  C.  BATRA  and  PEI-RONG  LIN 

Deparment  of  Engineering  Mechanics.  University  of  Missouri-Rolla,  Rolla.  MO  65401-0249.  U  S  A. 

Abstract — The  axisvmmetric  steady  state  deformations  of  an  infinite  cylindrical  rod  made  of  a 
ngid/perfectly  plastic  matenal  and  striking  a  known  cavity  in  a  rigid  target  are  analyzed  by  the  finite 
element  method.  The  contact  between  the  deforming  rod  and  the  target  surface  is  assumed  to  be 
smooth  It  is  found  that  the  axial  force  experienced  by  the  rod  depends  strongly  upon  the  square  of  its 
speed.  Results  computed  and  presented  graphically  include  the  velocity  field  in  the  deforming  region, 
the  dependence  of  the  shape  of  the  upset  head  of  the  striker  upon  its  speed,  and  the  distribution  of 
normal  tractions  upon  the  cavity  wall. 


INTRODUCTION 

Penetration  of  metal  targets  by  projectiles  is  influenced  by  such  variables  as  material 
properties,  impact  velocity,  projectile  shape,  target  support  position,  and  relative  dimensions 
of  the  target  and  the  projectile.  Recently,  emphasis  has  been  placed  on  kinetic  energy 
penetrators,  which  for  terminal  ballistic  purposes  may  be  considered  as  long  metal  rods 
traveling  at  high  speeds.  Wright  [1],  in  his  survey  article  on  long  rod  penetrators.  ellucidated 
vividly  some  of  the  problems  with  the  existing  penetration  models.  In  another  extensive  review 
article,  Backman  and  Goldsmith  [2]  discussed  superbly  the  work  done  in  penetration  mechanics 
until  1977.  Jonas  and  Zukas  [3]  reviewed  various  analytical  methods  for  the  study  of  kinetic 
energy  projectile-armor  interaction  at  ordance  velocities  and  placed  particular  emphasis  on 
three-dimensional  numerical  simulation  of  perforation.  Anderson  and  Bodner  [4]  have  recently 
reviewed  the  status  of  the  ballistic  impact  modeling.  A  penetration  model  that  is  not  too 
difficult  to  use  has  been  proposed  by  Ravid  and  Bodner  [5j.  They  studied  the  penetration 
problem  by  presuming  a  kinematically  admissible  flow  field  in  the  target  and  found  the 
unknown  parameters  by  utilizing  an  upper  bound  theorem  of  plasticity  modified  to  include 
dynamics  effects. 

In  an  attempt  to  shed  some  light  on  questions  raised  by  Wright  [1],  Batra  and  Wright  [6] 
recently  studied  an  idealized  penetration  problem  that  simulates  the  following  situation. 
Suppose  that  the  penetrator  is  in  the  intermediate  stages  of  penetration  so  that  the  active 
target/penetrator  interface  is  at  least  one  or  two  penetrator  diameters  away  from  either  target 
face,  and  the  remaining  penetrator  is  much  longer  than  several  diameters  and  is  still  traveling  at 
a  speed  close  to  its  striking  velocity.  This  situation  has  been  idealized  as  follows.  It  is  assumed 
that  the  rod  is  semi-infinite  in  length,  the  target  is  infinite  with  a  semi-infinite  hole,  the  rate  of 
penetration  and  all  flow  fields  are  steady  as  seen  from  the  nose  of  the  penetrator,  and  that  no 
shear  stress  can  be  transmitted  across  the  target/penetrator  interface.  This  last  assumption  is 
justified  on  the  grounds  that  a  thin  layer  of  material  at  the  interface  either  melts  or  is  severely 
degraded  by  adiabatic  shear.  These  idealizations  make  it  possible  to  decompose  the  penetration 
problem  into  two  parts  in  which  either  a  rigid  rod  penetrates  a  deformable  target  or  a 
deformable  rod  is  upset  at  the  bottom  of  a  hole  in  a  rigid  target.  Of  course,  in  the  combined 
case  the  contour  of  the  hole  is  unknown,  but  if  it  can  be  chosen  so  that  normal  tractions  match 
in  the  two  cases  along  the  entire  boundary  between  penetrator  and  target,  then  the  complete 
solution  is  known  irrespective  of  the  relative  motion  at  the  boundary.  Even  without  matching 
the  normal  tractions,  it  would  seem  that  valuable  qualitative  information  about  the  flow  field 
and  distribution  of  stresses  can  be  gained  if  the  chosen  contour  is  reasonably  close  to  those  that 
are  found  in  experiments. 

Whereas  Batra  and  Wright  [6|  studied  the  problem  of  the  deforming  target  and  a  rigid 
penetrator,  we  analyze  herein  the  companion  problem  of  a  deformable,  semi-infinite  and 
cylindrical  penetrator  striking  a  known  semi-infinite  cavity  in  an  infinite  and  rigid  target.  Only 
the  axisymmetric  and  steady  state  problem  in  which  the  penetrator  material  is  rigid/perfectlv 
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plastic  has  been  studied.  This  problem  is  more  challenging  than  the  one  studied  earlier  by 
Batra  and  Wright  (6j  because  of  the  presence  in  it  of  free  surfaces  whose  shapes  are  not  known 
a  priori.  We  hope  that  the  kinematic  and  stress  fields  found  in  this  study  would  help  in  devising 
and/or  checking  results  from  simpler  engineering  theories  of  penetration. 


FORMULATION  OF  THE  PROBLEM 

We  describe  the  deformations  of  the  cylindrical  rod  upset  at  the  bottom  of  a  semi-infinite 
cavity  in  an  infinite  rigid  target  with  respect  to  a  cylindrical  coordinate  system  with  origin  at  the 
center  of  the  hole  and  z-axis  pointing  into  the  rod.  Equations  governing  the  steady  state 
axisymmetric  deformations  of  the  rod  are 

divv  =  0,  (1) 

divo  =  pv,  (2.1) 

=  p(v  •  grad)v.  (2.2) 

Here  v  is  the  velocity  of  a  rod  particle,  p  is  the  mass  density  and  o  is  the  Cauchy  stress  tensor. 
Equation  (1)  implies  that  the  deformations  of  the  rod  are  isochoric.  and  eqn  (2)  expresses  the 
balance  of  linear  momentum.  The  operators  grad  and  div  signify  the  gradient  and  divergence 
operators  on  fields  defined  in  the  present  configuration.  We  neglect  the  elastic  deformations  of 
the  rod  and  assume  that  it  is  made  of  a  homogeneous  and  isotropic  material  that  obeys  the 
Von-Mises  yield  criterion  and  the  associated  flow  rule.  Thus  we  take  the  following  constitutive 
relation  for  o.  e.g.  see  Prager  and  Hodge  [7J. 


,,=  -pl*WiD ■ 

D  =  (grad  v  +  (grad  v)r)/2. 


(3) 

(4) 

(5) 


In  eqn  (3)  p  is  the  hydrostatic  pressure  which  cannot  be  determined  from  a  knowledge  of  the 
deformation  because  of  the  assumption  of  material  incompressibility.  1  is  the  unit  matrix,  o„  is 
the  flow  stress  of  the  material  of  the  rod  in  simple  compression  or  tension  and  /  is  the  second 
invariant  of  the  strain-rate  tensor  D.  Equation  (1)  and  the  equation  obtained  by  substituting 
eqn  (3)  into  (2)  are  the  field  equations  to  be  solved  for  v  and  p  subject  to  a  suitable  set  of 
boundary  conditions.  In  terms  of  the  non-dimensional  variables 


v  =  v/t/„.  d  =  o/o„,  p=p/o0,  r  =  r/r0,  z  =  z/r„ 


these  field  equations  are 


where 


div  v  =  0, 

-gradp  +  div((grad  v  +  (grad  v)r)/2V3  /)  =  a(v  ■  grad)v 


a  =  pv;Jo„ 


(6) 

(7) 

(8) 

(9) 


is  a  non-dimensional  number.  In  eqn  (6).  v„  is  the  speed  of  the  rod  and  rn  its  radius.  In  eqns  (7) 
and  (8)  and  hereafter  the  superimposed  bars  over  the  non-dimensional  variables  have  been 
dropped.  The  operators  grad  and  div  in  eqn  (8)  imply  the  gradient  and  divergence  operators  in 
terms  of  the  non-dimensional  coordinates.  We  note  that  there  is  only  one  non-dimensional 
number  a  that  governs  the  steady  state  deformations  of  the  rod. 

For  the  boundary  conditions  on  the  rod/cavity  interface,  we  assume  that 

t  •  (on)  =  0.  (10) 

v-n  =  0.  (H) 

where  n  and  t  are.  respectively,  a  unit  normal  and  unit  tangent  vector  on  the  interface.  The 
boundary  condition  (10)  represents  smooth  contact  between  the  rod  and  target.  This  appears 
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reasonable  since  a  thin  layer  of  material  at  the  interface  eithei  melts  or  is  severely  degraded  by 
adiabatic  shear.  The  boundary  conditon  (11)  represents  no  interpenetration  of  the  rod  material 
into  the  target  and  vice  versa.  On  the  free  surface  of  the  rod, 

cm  =  0.  (12) 

v  •  n  =  0.  (13) 

where  n  is  a  unit  outward  normal  to  the  surface.  The  boundary  condition  (13)  implies  that  the 
velocity  of  particles  on  the  surface  is  tangent  to  the  surface.  Out  of  the  boundary  conditions 
(12)  and  (13)  only  the  former  is  needed  to  specify  the  boundary-value  problem  completely 
provided  the  shape  of  the  free  surface  is  known.  Since  such  is  not  the  case,  we  use  eqn  (13)  to 
check  whether  the  presumed  free  surface  is  correct  or  not.  On  the  rod  cross-section  far  from 
the  cavity  bottom 

|v  +  e.|—»0  as  z  — (14) 

and  on  the  deformed  rod  material  at  the  cavity  outlet 

|on|-*0  as  \rz  +  z2\u2— *  *  (15) 

Equation  (14)  implies  that  the  end  of  the  rod  far  from  the  bottom  of  the  cavity  is  moving  with  a 
uniform  velocity  in  the  negative  z-direction  and  eqn  (15)  is  equivalent  to  the  statement  that  the 
surface  of  the  deformed  rod  near  the  cavity  outlet  is  traction  free.  In  order  to  state  the  problem 
more  precisely  one  should  specify  the  rates  at  which  quantities  indicated  in  eqns  (14)  and  (15) 
decay.  But  at  the  present  time  there  is  little  hope  of  establishing  an  existence  and  uniqueness 
theory  for  the  problem  and.  therefore,  specification  of  rates  of  decay  in  eqns  (14)  and  (15)  is 
not  required.  Here  we  seek  an  approximate  solution  of  the  problem  numerically. 


FINITE  ELEMENT  FORMULATION  OF  THE  PROBLEM 
To  solve  the  problem  numerically  we  usually  consider  a  finite  region.  Thus  the  region  of  the 
deformable  rod  analyzed  is  shown  in  Fig.  1  wherein  is  also  indicated  its  spatial  discretization. 
We  note  that  since  the  shape  of  the  free  surface  is  not  known  a  priori  we  estimate  one  and  will 
subsequently  check  whether  or  not  it  is  an  appropriate  one.  If  not,  we  will  modify  the  same  and 
keep  on  iterating  till  a  prespecified  criterion  is  met.  This  is  elaborated  upon  below.  The 
boundary  conditions  ( 10)— ( 12),  (14)  and  (15)  are  replaced  by  the  following: 


orz  =0,  vr  =  0  on  the  axis  of  symmetry  AB.  (16.1) 

v  •  n  =  0  and  t  •  (on)  =  0  on  the  cavity  surface  BC.  (16.2) 

on  =  0  on  the  free  surface  DEF.  (16.3) 

v.  =  -1.0,  ur  =  0  on  AF,  (16.4) 

v  =  u,n  and  t  •  (on)  =  0  on  the  outlet  surface  CD.  ( 16.5) 


As  before  n  is  a  unit  outward  normal  to  a  surface  and  t  is  a  unit  tangent  to  the  surface.  The 
boundary  condition  (16.5)  implies  that  the  rod  particles  at  the  exit  surface  CD  are  traveling 
normal  to  the  surface  with  a  uniform  speed  vr  and  the  tangential  traction  on  it  is  zero.  The 
value  of  v ,  is  computed  by  equating  the  amount  of  material  flowing  out  through  CD  to  that 
flowing  in  through  AF.  We  now  obtain  a  variational  statement  for  the  boundary  value  problem 
defined  by  eqns  (7),  (8)  and  (16). 

Let  ^  be  a  smooth  vector  valued  function  defined  on  the  region  R,  shown  enclosed  by 
ABCDEF  in  Fig.  1,  and  <pr  —  0  on  AB,  -  n  =  0  on  BC,  tj>  =  0  on  AF  and  9  •  n  =  0  on  CD. 
Also  let  v>  be  a  bounded  scalar  valued  function  defined  on  R.  Multiplying  both  sides  of  eqn  (7) 
with  tg,  taking  the  scalar  product  of  both  sides  of  eqn  (8)  with  <|>,  integrating  the  resulting 
equations  on  the  domain  R,  using  the  divergence  theorem,  the  traction  boundary  conditions  in 
(16)  and  the  above  stated  side  conditions  on  we  arrive  at  the  following. 

|  v/(divv)dK  =  0,  (17.1) 

Jr 
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Fig  1  The  region  studied  and  its  finite  element  discretization. 

I  p(div  dV'  -  J  j^/D:  (grad  0  +  (grad  0)r)dV'  =  aJ  {(v  •  grad)v}  •  0 dV.  (17.2) 

Here  a  single  dot  implies  the  scalar  product  between  two  vectors  and  the  symbol  :  stands  for 
the  scalar  product  between  two  second  order  symmetric  tensors.  Thus  a  weak  formulation  of 
the  boundary-value  problem  given  by  eqns  (7),  (8)  and  (16)  is  to  find  p  and  v  defined  on  R  such 
that  eqns  (17)  holds  for  every  and  0.  y  and  grad  0  are  square  integrable  over  R,  0  satisfies 
the  homogeneous  essential  boundary  conditions  stated  above  and  v  satisfies  the  essential 
boundary  conditions  outlined  in  eqns  (16). 

We  refer  the  reader  to  Becker  el  al.  [8]  for  details  of  obtaining  a  finite  element  solution  of 
eqns  (17).  Suffice  it  to  say  that  eqns  (17)  are  solved  by  using  the  following  iterative  scheme. 

[  V(divv/n)dV  =  0,  (18.1) 

Jr 

j  p'"(div0)dV/  -  J  -^-^rrD'":(grad0  +  (grad0)7")dV  =  a'|^{(vw",-grad)vm}-0dV/  (18.2) 

Here  m  is  the  interation  number.  Note  that  the  hydrostatic  pressure  p  appears  linearly  and  its 
previous  values  are  not  needed.  The  initial  velocity  field  is  taken  to  be  zero  for  small  values  of 
a.  If  the  solution  is  known  for  some  a,  it  is  taken  as  the  initial  solution  when  solving  the 
problem  for  the  next  higher  value  of  a.  The  iteration  process  in  (18)  is  stopped  if.  at  each  nodal 
point. 

IK -v'"-1!!  <0.01  ||v'"-,||  (19) 
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where  ||v||2  =  (it;  +  v2.).  When  the  velocity  field  corresponding  to  (a  -  1)  was  taiten  as  the 
initial  estimate  for  the  solution  to  be  computed  for  an  assigned  value  of  a,  it  took  nearly  10 
iterations  for  the  convergence  criterion  (19)  to  be  satisfied. 


TREATMENT  OF  THE  FREE  SURFACE 
In  the  preceding  finite  element  formulation  of  the  problem  the  shape  of  the  free  surface  DEF 
was  presumed  to  be  given.  Since  it  is  not  known  a  priori  we  presume  one.  solve  the  problem 
and  then  see  if  the  condition  (13)  is  satisfied  on  it  or  not.  Improving  upon  the  presumed  shape 
of  the  free  surface  DEF  if  (13)  is  not  satisfied  required  a  considerable  effort.  We  note  that 
problems  in  which  the  shape  of  the  free  surface  is  to  be  determined  as  a  part  of  the  solution  of 
the  problem  have  been  solved  by  Zienkiewicz  et  al.,  [9],  Oden  and  Lin  [10).  and  Batra  et  al. 
(11).  For  us  the  following  technique  proved  quite  effective. 

Starting  from  the  point  F  we  found  the  parabolic  curve  passing  through  F  and  the  next  two 
nodes  downstream  from  it.  We  then  found  the  tangent  to  the  curve  and  a  unit  normal  n  at  the 
middle  node.  We  evaluated  the  magnitude  of  the  error 

e=v-n/||v||  (20) 

at  the  middle  node  and  repeated  this  procedure  for  all  of  the  nodes.  Note  that  the  node  F  is 
fixed  and  the  last  node  D  on  the  free  surface  was  assumed  to  be  on  the  straight  line  passing 
through  the  two  nodes  immediately  preceding  to  it  on  the  free  surface  DEF.  That  is.  the 
curvature  of  DE  near  the  point  D  was  presumed  to  be  zero.  Points  where  the  magnitude  of  the 
error  e  was  not  smaller  than  a  preassigned  small  number  were  moved  along  n  or  -n  according 
as  e  was  positive  or  negative  and  the  distance  moved  was  proportional  to  the  magnitude  of  e.  A 
check  was  made  to  ensure  that  two  different  nodes  on  the  presumed  free  surface  are  not 
mapped  into  the  same  location  when  modifying  the  shape  of  the  free  surface.  The  procedure 
was  repeated  until  the  error  e  at  each  node  on  the  free  surface  DEF  was  less  than  0.1.  It  was 
found  computationally  very  efficient  to  switch  to  the  following  method  for  adjusting  the  free 
surface  subsequently.  Let  H  and  K  be  two  successive  nodes  on  the  free  surface  downstream 
from  G.  N  be  the  normal  vector  to  the  previously  assumed  free  surface  at  K.  T  be  a  vector 
parallel  to  v  computed  at  K  and  we  wish  to  find  the  new  location  K’  of  K.  It  is  at  the  point  of 
intersection  K'  of  N  with  the  circular  arc  that  passes  through  G  and  H.  and  has  a  tangent  vector 
at  K'  parallel  to  T.  This  eliminates  the  iikelihood  of  two  nodes  ending  up  at  the  same  location 
during  a  refinement  of  the  shape  of  the  free  surface.  With  this  technique  the  magnitude  of  e  at 
each  node  point  on  the  surface  DEF  could  be  reduced  to  less  than  5%  and  the  average  of  |e| 
for  all  nodes  on  DEF  to  less  *han  1.5%  in  at  most  six  iterations. 

COMPUTATION  AND  DISCUSSION  OF  RESULTS 
The  finite  element  program  developed  earlier  [6)  to  solve  the  companion  problem  of  a  rigid 
rod  penetrating  into  a  deformable  target  was  modified  to  study  the  present  problem.  It  employs 
6-noded  triangular  elements  and  within  each  element  the  velocity  (pressure)  field  is 
approximated  by  a  quadratic  (linear)  function  defined  in  terms  of  its  values  at  6  (3  corner) 
nodes.  Thus  both  the  velocity  and  pressure  fields  are  continuous  across  interelement 
boundaries.  We  note  the  element  satisfies  the  Babu§ka-Brezzi  [12]  condition.  The  velocity 
boundary  condition  in  (16.2)  is  accounted  for  by  using  the  method  of  Lagrange  multipliers.  The 
sample  problem  used  to  establish  the  validity  of  the  code  has  been  discussed  in  [6]. 

In  Fig.  2  is  plotted  the  computed  velocity  field  for  or  =  5.1.  The  details  of  the  velocity  field 
within  the  vicinity  of  the  point  where  the  curvature  of  the  free  surface  changes  sharply  are  also 
shown.  As  is  rather  obvious  from  the  plotted  results  the  Lagrange  multipier  method  is  quite 
effective  in  satisfying  the  essential  boundary  condition  (16.2)  on  the  cavity  surface.  We  note 
that  on  the  computed  free  surface  the  velocity  of  points  is  along  the  tangent  to  the  surface.  All 
of  the  results  presented  herein  are  for  a  fixed  shape  of  the  cavity  given  by 

c  =  0.04r4. 

This  shape  of  the  cavity  was  chosen  after  several  trials;  the  selection  criterion  being  that  the 
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Fig.  2.  The  computed  velocity  field  for  a  =  5. 1. 


part  of  the  surface  over  which  the  deformed  rod  particles  exhibit  a  tendency  to  separate  away 
from  the  surface  BC  is  as  small  as  possible.  Ideally  the  contact  should  be  maintained  over  all  of 
BC  but  such  a  goal  proved  essentially  impossible  to  attain  especially  for  different  striking 
speeds  of  the  rod. 

The  velocity  field  for  other  values  of  a  was  quite  similar  to  that  shown  in  Fig.  2.  As  the  speed 
of  the  rod  increased  the  free  surface  moved  away  from  the  fixed  cavity  so  that  the  thickness  of 
the  region  at  the  outlet  increased.  The  shapes  of  the  free  surface  for  different  values  of  a  are 
shown  in  Fig.  3.  The  radius  of  curvature  near  the  bottom  of  the  surface  decreases  sharply  as  the 
speed  of  the  rod  is  increased.  The  plot  of  the  thickness  at  the  outlet  versus  a  in  Fig.  3  shows 


a  =  I  8 
a  =  22 
a  =26 


(  a  )  - a  =3  l 


Fig  3  a.  Shapes  of  the  free  surface  for  different  values  of  a.  b.  Thickness  of  the  outlet  region 

versus  a. 
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Fig.  4.  Variation  of  the  second  invariant  I  of  the  stretching  tensor  D  in  the  deforming  region. 

that  the  thickness  increases  linearly  with  a.  Since  the  speed  of  the  rod  at  the  inlet  is  always 
taken  to  be  1.0,  the  speed  of  the  particles  at  the  outlet  decreases  so  as  to  satisfy  the  balance  of 
mass. 

Figure  4  depicts  the  variation  of  the  second  invariant  /  of  the  strain-rate  tensor  D  in  the 
deforming  region.  As  one  would  expect  severe  deformations  take  place  in  the  region  near  the 
bottom  of  the  cavity.  Larger  values  of  I  occur  near  the  bottom  of  the  free  surface  where  the 
flow  reverses.  For  typical  values  of  the  radius  of  the  penetrator  these  non-dimensional  values  of 
l  are  to  be  multiplied  by  10s  indicating  thereby  that  peak  strain-rates  of  the  order  of  K)5— 10°  s~ 1 
occur  in  the  vicinity  of  the  point  on  the  free  surface  where  the  flow  reverses.  We  note  that  the 


Fig.  5.  Principal  stresses  in  the  deforming  region. 
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non-dimensional  values  of  /  and  their  locations  in  the  deforming  region  did  not  vary  much  with 
a.  In  Figs  5  and  6  are  plotted,  respectively,  the  principal  stresses  and  the  hydrostatic  pressure 
in  the  deforming  region  for  two  different  values  of  a.  The  lines  are  oriented  along  the  axes  of 
principal  stresses  and  their  lengths  are  proportional  to  the  magnitudes  of  principal  stresses  at 
that  point.  The  arrows  at  the  ends  of  a  line  indicate  that  the  corresponding  principal  stress  is 
tensile  otherwise  it  is  compressive.  Whereas  for  a  =  1.8  the  material  particles  whose  velocity  in 
the  z -direction  is  opposite  to  that  of  the  incoming  rod  experience  tensile  stresses  such  is  not  the 
case  at  the  higher  value  of  a.  This  transition  seems  to  take  place  around  a  =  3.6.  For  a  =  1.8 
one  of  the  principal  stresses  at  points  in  the  region  between  the  free  surface  and  the  bottom  of 
the  cavity  is  compressive;  the  other  one  is  essentially  zero.  At  the  higher  value  of  a  both 
principal  stresses  are  compressive.  This  is  possibly  due  to  the  increase  in  the  hydrostatic 
pressure  (cf.  Fig.  6)  with  the  speed  of  the  rod.  Since  the  strain  rates  at  points  near  the  inlet  and 
the  outlet  surface  are  negligibly  small,  and  the  values  of  stresses  at  these  points,  as  given  by 


Fig.  h  Distribution  of  the  hydrostatic  pressure  p  in  the  deforming  region 
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Arc  length 

Fig.  7.  Variation  of  the  normal  and  axial  tractions  on  the  cavity  wall. 

eqn  (3),  are  computed  by  taking  the  ratio  of  two  small  numbers,  the  magnitudes  of  stresses  at 
these  points  may  not  be  very  accurate.  This  is  one  possible  explanation  for  the  small  tensile 
stresses  at  points  near  the  inlet  and  outlet  surfaces.  At  a  few  isolated  points  in  the  vicinity  of 
the  bounding  surfaces  near  the  inlet  and  outlet  regions  the  computed  values  of  the  hydrostatic 
pressure  were  relatively  small  negative  numbers. 

There  is  a  tendency  for  the  material  particles  to  leave  the  cavity  wall  near  the  outlet  region. 
This  is  indicated  by  the  negative  values  of  the  computed  normal  traction  at  these  points.  The 
normal  tractions  and  the  axial  tractions  at  different  points  on  the  cavity  wall  are  plotted  in  Fig. 
7.  The  arc  length  along  the  cavity  surface  is  measured  from  the  point  of  intersection  of  the 
cavity  and  the  axis  of  the  rod.  For  the  same  cavity  surface  the  point  where  separation  tends  to 
occur  moves  away  from  the  centroidal  axis  of  the  rod  as  the  speed  of  the  rod  is  increased.  We 
note  that  in  our  computations  material  particles  were  not  allowed  to  leave  the  cavity  surface. 
Since  the  magnitude  of  the  normal  tractions  at  these  points  is  rather  small  and  the  distance  of 
these  points  from  the  center  of  the  rod  is  of  the  order  of  2r0,  the  computed  results  especially  in 
the  severely  deformed  region  near  the  bottom  of  the  cavity  are  quite  meaningful.  It  is  now 
obvious  that  to  keep,  say,  the  horizontal  distance  of  the  point  on  the  cavity  wall  where 
separation  tends  to  occur  from  the  centroidal  axis  of  the  rod  constant,  the  shape  of  the  cavity 
should  be  adjusted  as  the  rod  speed  is  varied.  The  axial  traction  plotted  in  Fig.  7  indicates  that 
the  contribution  to  the  total  axial  force  from  points  where  flow  separation  would  have  occurred 
is  very  small.  This  is  due  to  the  fact  that  the  angle  between  the  cavity  surface  at  these  points 
and  the  axis  of  the  rod  is  very  small  so  that  the  normal  to  the  cavity  surface  at  these  points  is 
essentially  perpendicular  to  the  rod  axis. 

How  the  total  axial  force  acting  on  the  cavity  surface  and  experienced  by  the  rod  depends 
upon  a  is  shown  in  Fig.  8.  Unlike  the  case  of  a  rigid  rod  penetrating  into  a  rigid/perfectly 
plastic  target  where  the  axial  force  acting  on  the  penetrator  depended  weakly  on  the  penetrator 
speed,  here  the  dependence  of  the  axial  force  upon  a  is  rather  strong.  We  note  that  whereas 
herein  the  calculations  have  been  performed  with  one  cavity  surface,  in  the  actual  penetration 
problem  the  surface  separating  the  deforming  penetrator  and  target  regions  will  alter  with  the 
speed  of  the  penetrator  when  all  other  parameters  are  kept  fixed.  In  the  approximate  theory  of 
Tate  [13]  the  axial  force  acting  at  the  penetrator/target  interface  is  presumed  to  be  constant. 
We  have  not  investigated  the  dependence  of  the  axial  force  upon  the  shape  of  the  assumed 
cavity.  Since  there  is  very  little  experimental  data  available  in  the  open  literature  it  is  hard  to 
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Fig.  8.  Dependence  of  the  total  axial  force  acting  on  the  cavity  wall  upon  a. 


judge  whether  or  not  the  cavity  surface  envisaged  here  is  close  to  those  observed 
experimentally. 

We  note  that  strain-rates  and  therefore  the  rod  deformations  are  negligibly  small  near  the 
inlet  and  outlet  regions.  This  ensures  that  the  assumed  finite  region  considered  for  solving  the 
problem  numerically  is  sufficient.  The  computed  solution  seemed  to  be  stable  in  the  sense  that 
either  superimposed  small  perturbations  died  out  to  zero  or  following  slightly  different  paths 
gave  the  same  solution. 


CONCLUSIONS 

For  the  fixed  cavity  the  thickness  at  the  outlet  increases  linearly  with  a,  and  the  total  axial 
force  acting  on  the  cavity  wall  depends  strongly  upon  a.  Most  severe  deformations  occur  in  the 
region  near  the  cavity  bottom  and  the  point  where  the  curvature  of  the  free  surface  changes. 
Peak  strain-rates  in  the  range  105-10<’s-1  invariably  occurred  at  or  near  the  bottom  E  (Fig.  1) 
of  the  free  surface.  The  hydrostatic  pressure  increases  considerably  with  a.  Whereas  at  lower 
speeds  of  the  rod  tensile  stresses  developed  at  points  on  the  exist  side  of  the  flow  these  were 
overcome  by  the  increase  in  the  value  of  the  hydrostatic  pressure  at  higher  speeds  of  the  rod. 
The  point  on  the  cavity  wall  where  the  flow  had  a  tendency  to  separate  from  the  wall  moved 
away  from  the  axis  of  the  rod  as  the  striking  speed  increased. 
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Abstract — Steady  state  axisymmetnc  deformations  of  a  homogeneous,  isotropic,  compressible  and 
rigid  perfectly  plastic  target  being  penetrated  by  a  rigid  cylindrical  penetrator  with  a  hemispherical 
nose  are  studied  by  the  finite  element  method.  The  steady  state  is  reached  with  respect  to  an  observer 
situated  on  the  penetrator  nose  and  moving  with  it.  Tillotson's  equation,  restricted  to  mechanical 
deformations,  is  used  to  express  the  pressure  as  a  function  of  the  mass  density.  Contact  between  the 
penetrator  and  the  target  is  assumed  to  be  smooth.  The  effect  of  compressibility  of  the  material  is 
delineated  by  comparing  results  for  compressible  and  incompressible  materials.  Also  studied  is  the 
effect  of  the  penetrator  speed  on  target  deformations 


INTRODUCTION 

In  an  attempt  to  shed  some  light  on  the  validity  of  the  approximations  made  in  simple  theories 
of  penetration  due  to  Alekseevskii  [1]  and  Tate  [2j,  and  to  understand  better  the  significance  of 
various  kinematic  variables.  Batra  and  Wright  [3]  studied  in  detail  a  steady  state  penetration 
problem  that  simulates  the  following  situation.  Suppose  that  the  penetrator  is  in  the 
intermediate  stages  of  penetration  so  that  the  active  target/penetrator  interface  is  at  least  one 
or  two  penetrator  diameters  away  from  either  target  face,  and  the  remaining  penetrator  is 
much  longer  than  several  diameters  and  is  still  travelling  at  a  uniform  speed.  This  stage  of 
penetration  can  be  idealized  as  one  in  which  deformations  of  the  target  appear  to  be  steady  to 
an  observer  situated  on  the  penetrator  nose.  Wright  and  Batra  [3]  presumed  that  the  target  is 
made  of  a  rigid/perfectly  plastic  material,  and  is  being  penetrated  by  a  long  cylindrical  rigid  rod 
with  a  hemispherical  nose.  Subsequently  Batra  [4}  showed  that  the  axial  resisting  force 
experienced  by  the  rigid  penetrator  is  considerably  reduced  if  its  nose  shape  is  ellipsoidal  rather 
than  hemispherical  and  also  investigated  the  effect  of  the  dependence  of  the  flow  stress  upon 
the  strain-rate.  He  [5]  has  extended  this  work  to  the  case  when  the  target  material  is  thermally 
softening  but  strain  and  strain-rate  hardening. 

Pidsley  [6]  recently  studied  a  complete  penetration  problem  in  which  both  target  and 
penetrator  materials  were  assumed  to  be  compressible  rigid/perfectly  plastic.  Whereas  he  gave 
a  detailed  numerical  solution  for  one  set  of  material  and  geometric  parameters,  we  study  the 
steady  state  penetration  problem  similar  to  the  one  analyzed  by  Batra  and  Wright  [3]  and 
investigate  the  effect  of  the  speed  of  the  penetrator  and  the  compressibility  of  the  target 
material.  We  hope  that  the  kinematic  and  stress  fields  found  in  this  study  will  be  useful  in 
identifying  key  variables  to  be  included  in  simpler  engineering  theories  of  target  penetration. 
Since  the  continuity  equation  has  no  diffusive  term,  it  is  more  challenging  to  solve  the  problem 
numerically  for  compressible  materials.  Pidsley  used  the  finite  difference  code  HELP  developed 
by  Hageman  and  Walsh  [7]  to  solve  the  problem.  However,  we  use  the  finite  element  method 
and  have  developed  the  requisite  code. 

We  note  that  no  failure  or  fracture  criterion  is  included  in  our  study.  Thus  the  target  material 
is  assumed  to  undergo  unlimited  amount  of  deformations. 


FORMULATION  OF  THE  PROBLEM 

We  describe  the  axisymmetric  deformations  of  the  target  with  respect  to  a  set  of  cylindrical 
coordinate  axes  with  origin  at  the  center  of  the  hemispherical  nose  of  the  rigid  penetrator  and 
z-axis  pointing  into  the  target.  Equations  governing  the  steady-state  deformations  of  the  target 
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are 


div(pv)  =  0.  (1.1) 

div  a  =  p(v  •  grad)v.  (1.2) 

Here  p  is  the  current  mass  density  of  a  target  material  particle,  v  is  the  velocity  relative  to  an 
observer  situated  on  the  penetrator  nose  tip  and  moving  with  it.  o  is  the  Cauchy  stress  tensor, 
and  the  operators  grad  and  div  signify  the  gradient  and  divergence  operators  on  fields-defined 
in  the  present  configuration.  Equations  (1.1)  and  (1.2)  express,  respectively,  the  balance  of 
mass  and  the  balance  of  linear  momentum  with  zero  body  forces. 

We  presume  that  the  target  material  is  compressible  and  obeys  the  Von-Mises  yield  criterion 
and  the  associated  flow  rule.  That  is+ 

a=  -p(p)l  +  y^D. 

D  =  D  -  ^  (tr  D)l. 

D  =  (grad  v  +  (grad  v)r)/2. 

/2  =  ^tr(D;). 

p(p)  =  a(—  -  l)  +  b(  —  ~  l)  . 

Vp,.  '  Vo  ' 

In  these  equations  o„  is  the  flow  stress  of  the  target  material  in  simple  compression,  p,>  is  the 
mass  density  of  the  undeformed  target  material.  A  and  B  are  material  constants.  D  is  the 
strain-rate  tensor  and  D  its  deviatoric  part.  Equation  (2.5)  is  obtained  from  the  Tillotson 
equation  by  neglecting  the  parts  that  apply  to  hot  expanded  metal  states  and  the  change  in 
temperature  of  a  material  point.  The  constant  A  is  related  to  the  bulk  modulus  of  the  material. 
Equation  (2.1)  may  be  regarded  as  a  constitutive  relation  of  a  compressible  Non-Newtonian 
fluid  with  shear  viscosity  given  by  oJ2V?>l.  Equation  (2)  when  substituted  into  (1.2)  gives 


(2.1) 

(2.2) 

(2.3) 

(2.4) 

(2.5) 


-gradp  +  (div((grad  v  +  (grad \)r)/2 /) 

-  grad((div  v)/3/))  =  p(v  •  grad)v  (3) 

which  along  with  (1.1)  are  the  field  equations  for  p  and  v. 

We  now  non-dimensionalize  the  variables  as  follows: 

f  -  r/r„,  i  =  c/r„,  v  =  v/i>„.  d  =  o/o„.  p=p/a„.  p  =  p/p„, 

(4) 

A  =  A/a„.  B  =  8/o„. 

Here  r„  is  the  radius  of  the  cylindrical  part  of  the  penetrator.  v„  its  speed  in  the  c-direction  and 
the  pair  (r.  r)  denotes  the  cylindrical  co-ordinates  of  a  target  particle.  Hereafter  we  use  only 
the  non-dimensional  variables  and  drop  the  superimposed  bars.  Equations  (1.1)  and  (3)  in 
terms  of  the  non-dimensional  variables  are 


(grad  p)  •  v  +  p(div  v)  =  0. 


-  grad  p  +  — -  (div((grad  v  +  (grad  v)'  )/2 1)  -  grad((div  v)/3/))  =  a(v  •  gradlv. 
V  3 


V3 

where 

is  a  non-dimensional  number. 


n  =  f)„vr,/a„ 


(5.1) 

(5.2) 

(6) 


*  This  constilulivc  relation  is  discussed  briefly  in  the  note  at  the  end  ol  the  paper 
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For  the  boundary  conditions  on  the  target/penetrator  interface  we  take 

v-n  =  0.  (7.1) 

t  •  (on)  -  0,  (7.2) 

where  n  and  t  are.  respectively,  unit  normal  and  tangent  vectors  at  a  point  on  the  interface. 
The  boundary  conditions  (7)  ensure  that  target  particles  do  not  penetrate  into  the  rigid 
penetrator  and  the  contact  surface  is  frictionless.  At  points  far  away  from  the  penetrator  nose, 
we  assume  that 

|v  +  e|— *0  as  (r:  +  r2)1 2— » *.  (8.1) 

|  on  |  — *  0  as  z  — * -*,  r>l.  (8.2) 

In  eqn  (8.1)  e  is  a  unit  vector  along  the  z-axis.  This  boundary  condition  states  that  target 
particles  far  from  the  penetrator  but  not  on  the  back  surface  appear  to  move  as  a  rigid  body 
and  those  on  the  back  surface  are  traction  free.  In  order  to  state  the  problem  precisely  one 
needs  to  specify  the  rates  at  which  quantities  in  eqn  (8)  decay  to  zero.  We  now  assume  that  the 
nonlinear  and  coupled  equations  (5)  under  the  boundary  conditions  (7)  and  (8)  have  a  solution 
and  find  an  approximation  to  that  solution  by  the  finite  element  method. 


FINITE  ELEMENT  FORMULATION  OF  THE  PROBLEM 

For  a  numerical  solution  of  the  problem  we  first  recall  that  the  target  deformations  are 
assumed  to  be  axisvmmetnc  and  therefore  consider  the  finite  region  shown  in  Fig.  1  and 

l 
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impose  the  following  boundary  conditions: 

a..  =0,  v,  =  0  on  the  bottom  surface  AB,  (9.1) 

r -011  =  0.  v  •  n  =  0  on  the  target/penetrator  interface  BCD,  (9.2) 

ar:  =  0,  vr  =  0  on  the  axis  of  symmetry  DE,  (9.3) 

t'r  =  0,  v.  =  -1.0  on  the  boundary  surface  EFA.  (9.4) 


The  adequacy  of  the  studied  region  is  verified  by  ensuring  that  the  computed  results  especially 
in  the  vicinity  of  the  target/penetrator  interface  are  unaffected  by  the  location  of  the  boundary 
EFA. 

Referring  the  reader  to  Becker  et  al.  [8]  for  details,  we  note  that  a  weak  formulation  of  the 
problem  defined  on  the  region  R  (shown  enclosed  by  ABCDEFA  in  Fig.  1)  by  eqns  (5)  and 
boundary  conditions  (9)  is  that  equations 

I  y;((grad  p)  •  v  +  p(div  v))  du  =  0,  (10.1) 

Jr 

|  p(div  <]>)  du  -  j  j  D :  (grad  $  +  (grad  0)r)  dt>  =  a  j  ((v  •  grad)v)  •  $  du,  (10.2) 

hold  for  arbitrary  smooth  functions  ip  and  $  defined  on  R  such  that  <pr  =  0  on  AB,  =  0  on 
EFA  and  <|>  •  n  =  0  on  the  target/penetrator  interface  BCD.  In  these  equations  A :  B  =  tr(ABr) 
for  linear  transformations  A  and  B.  Since  these  equations  are  nonlinear  in  v  and  p,  the 
following  iterative  technique  has  been  employed.  At  the  ith  iteration,  equations 


y>((grad  p')  •  V  '  +  p'  '(div  v1))  dv  =  0, 

Jr 

^  :  (grad  <|i  +  (grad  $)'r)  du  +  a  J*  ((v'~'  •  grad)v')  •  0  dv 

=  |p(p'"‘)(d iv«j»)du  +  ^  [(tr  D'_1)div  <p]  du 


(10.3) 


(10.4) 


are  solved  for  u'  and  p'.  The  iterative  process  is  stopped  when  at  each  nodal  point 

**  llu'-u'-'ll  +  lp'-p'-'lstdlu'-'ll  +  lp'-'l]  (10.5) 

where  ||v||2  =  v;  +  u:,  and  e  is  a  preassigned  small  number. 

The  lack  of  a  diffusive  term  in  the  balance  of  mass  eqn  (5.1)  necessitates  that  the  test 
functions  y  and  the  trial  solutions  p  be  chosen  from  different  functional  spaces.  This  is  usually 
referred  to  as  Petrov-Galerkin  formulation  [9],  We  use  9-noded  rectangular  elements  for  v  and 
the  four-noded  rectangular  elements  for  p  and  employ  the  basis  functions  given  by  Heinrich  et 
al.  [  10]  to  generate  the  test  functions  y>.  These  basis  functions  involve  four  constants  and  the 
rate  of  convergence  of  the  solution  depends  rather  strongly  upon  the  values  of  these  constants. 
The  test  functions  and  the  trial  solutions  v  are  taken  from  the  same  space  of  functions.  Thus 
we  have  used  the  Petrov-Galerkin  formulation  for  the  continuity  equation  and  the  Galerkin 
formulation  for  the  balance  of  linear  momentum. 


COMPUTATION  AND  DISCUSSION  OF  RESULTS 

A  computer  code  based  on  eqns  (10.3)  and  (10.4),  and  employing  9-noded  rectangular 
elements  has  been  developed.  The  two  components  of  the  velocity  are  taken  as  unknowns  at 
each  node  and  the  mass  density  is  assumed  to  be  unknown  only  at  the  four  corner  nodes.  The 
accuracy  of  the  computer  code  was  established  by  solving  a  hypothetical  problem  involving  the 
flow  of  a  compressible  Navier-Stokes  fluid  in  a  circular  pipe  and  achieving  a  favorable 
comparison  between  the  computed  and  analytical  results.  The  sample  problem  studied  and  the 
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comparison  between  the  computed  and  analytical  results  is  given  in  the  Appendix.  Even 
though  the  formulation  of  the  problem  does  not  require  that  the  mass  density  be  prescribed  on 
any  part  of  the  boundary,  the  numerical  solution  of  the  problem  necessitated  that  it  be  assigned 
values  on  a  part  of  the  boundary.  In  the  results  presented  below  the  non-dimensional  mass 
density  was  set  equal  to  1.0  on  the  part  EFA  of  the  boundary.  Also  the  boundary  condition 
v  •  n  =  0  on  the  target/penetrator  interface  BCD  was  satisfied  by  using  the  method  of  Lagrange 
multipliers.  The  finite  element  discretization  of  the  region  studied  is  depicted  in  Fig.  1.  Note 
that  only  one  non-dimensional  number  a  governs  the  deformations  of  the  target  material.  Thus 
all  of  the  results  below  are  expressed  in  terms  of  a.  However  the  compressibility  of  the 
material  is  governed  by  the  values  of  A  and  B  in  eqn  (2.5)  which  we  took  as  A  =  144.231. 
B  =  125.0.  These  values,  taken  from  Pidsley’s  [6j  paper,  are  for  aluminum.  For  a  prescribed 
value  of  a,  the  problem  was  first  solved  by  presuming  that  the  target  material  is 
incompressible.  This  solution  for  v,  and  v.  and  p  =  1.0  everywhere  was  taken  as  the  initial 
estimate  of  the  solution  for  the  compressible  target  material.  We  recall  that  the  basis  functions 
for  y;  involve  four  constants.  The  number  of  iterations  required  to  obtain  the  converged 
solution  depended  rather  strongly  upon  the  values  selected  for  these  constants.  In  general, 
however,  the  number  of  iterations  needed  for  the  solution  to  converge  increased  with  the 
increase  in  the  value  of  a.  The  parameter  e  in  the  convergence  criterion  (10.5)  was  assigned  the 
value  0.04. 

In  Fig.  2  are  plotted  the  variation  of  the  normal  stress,  tangential  velocity  and  the  second 
invariant  /  of  the  strain-rate  tensor  on  the  penetrator  nose  for  or  =  5.36  and  for  compressible 
and  incompressible  target  materials.  The  compressibility  of  the  target  material  affects  most  the 
values  of  /  near  the  penetrator  nose-tip.  Over  most  of  the  nose  surface  the  values  of  /  and  the 
normal  stress  are  lower  for  the  compressible  target  material  as  compared  to  that  when  the 
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Fig  2.  The  variation  of  the  normal  stress,  tangential  velocity  and  the  second  strain-rate  invariant  /  on 
the  penetrator  nose  for  a  =  5.36. - incompressible  material. . compressible  material 
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target  material  is  modelled  as  incompressible.  Whereas  the  peak  value  of  /  at  the  penetrator 
nose-tip  drops  by  8.78%  when  the  target  material  is  presumed  to  be  compressible,  that  of  the 
normal  stress  decreases  by  only  3.19%.  In  the  immediate  vicinity  of  the  nose  tip  the  tangential 
velocities  for  compressible  and  incompressible  targets  essentially  coincide  with  each  other. 
However  for  the  angular  position  8  greater  than  12°,  the  tangential  velocity  for  the 
compressible  target  is  lower  than  that  for  the  incompressible  target. 

Figure  3  depicts  the  effect  of  the  compressibility  of  the  target  material  on  the  variation  of  the 
normal  stress,  strain-rate  invariant  /  and  the  relative  2-velocity  along  the  axial  line.  As 
evidenced  by  the  values  ot  /,  most  of  the  target  deformations  occur  within  a  distance  of  1.5 
times  the  radius  of  the  penetrator  for  both  compressible  and  incompressible  target  materials.  In 
each  case  the  values  of  /  drop  off  to  zero  at  points  on  the  axial  line  whose  distance  from  the 
nose  surface  is  greater  than  twice  the  radius  of  the  penetrator.  Typically  these  non-dimensional 
values  of  /  need  to  be  multiplied  by  105  to  get  their  dimensional  counterparts.  Thus  peak 
strain-rates  of  the  order  of  105s-1  occur  for  both  compressible  and  incompressible  target 
materials.  The  values  of  absolute  z-velocity  decay  to  zero  a  little  bit  slowly  for  compressible 
targets  as  compared  to  that  for  incompressible  targets.  The  rate  of  decay  of  the  absolute  values 
of  a.,  along  the  axial  line  in  the  two  cases  is  nearly  the  same.  The  values  of  | a,. |  are  less  when 
the  target  material  is  modeled  as  compressible  as  compared  to  that  when  it  is  taken  to  be 
incompressible.  The  difference  between  the  two  values  is  primarily  due  to  the  difference  in  the 
values  of  the  hydrostatic  pressure  in  the  two  cases.  For  0  s  z  ^  2,  this  difference  in  the  values  of 
p  stayed  nearly  constant  and  equalled  0.48. 

Figure  4  shows  the  variation  of  v.  with  r  at  2  =0  and  z  =  -5.0  and  also  the  dependence  of 
the  axial  resisting  force  F  experienced  by  the  penetrator  upon  or.  For  both  compressible  and 
incompressible  target  materials  the  axial  resisting  force  F  depends  linearly  upon  o  and  the  two 
lines  are  nearly  parallel  to  each  other.  In  each  case  the  dependence  of  F  upon  a  is  rather  weak 
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Fie  3  Variations  of  (—<>..(.  ; -velocity  and  the  second  strain-rate  invariant  /  on  the  axial  line  tor 
a  =  - incompressible  material.  compressible  material 
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Fig  4  Dependence  of  the  axial  force  upon  a.  Variations  ol  i\  with  r  on  planes  c  =  0  and  ;  =  -5.0. 
a  =  5..Vi.  —  incompressible  material.  compressible  material. 


and  explains  why  the  choice  of  constant  target  resistance  in  the  simple  theory  of  Tate  [2]  gives 
good  qualitative  results.  Note  that  the  dimensional  values  of  F  equal  (^rr,a„)  times  its 
non-dimensional  values.  The  plots  of  t».  versus  r  indicate  the  deformations  of  the  target  spread 
farther  to  the  side  of  the  penetrator  than  ahead  of  it.  This  is  true  both  for  compressible  and 
incompressible  target  materials.  Whereas  target  points  on  the  axial  line  and  ahead  of  the 
penetrator  situated  at  a  distance  of  three  times  the  penetrator  radius  have  zero  absolute 
z-velocity,  that  on  the  sides  of  the  penetrator  and  situated  at  the  same  distance  from  the 
penetrator  surface  have  non-zero  :: -velocity.  The  variation  of  u_.  with  r  on  the  surface  z  =  -5.0 
indicates  that  the  material  in  the  region  1  <  r  s  3.5  extrudes  as  a  rigid  block. 

The  value  of  the  compression  c,  defined  as  ( 1  -  p"1).  at  the  nose  tip  was  found  to  depend 
linearly  upon  a,  and  for  4  s  a  s  8  the  two  are  related  by 

c  =  (4.75  +  0.575o')/l(X). 

After  having  delineated  the  difference  in  the  solutions  for  compressible  and  incompressible 
target  materials  we  next  studied  the  effect  of  the  penetrator  speed  upon  the  deformations  of  the 
compressible  target.  We  recall  that  the  corresponding  results  for  the  incompressible  target 
material  have  been  given  by  Batra  and  Wright  [3 J.  Results  presented  in  Fig.  5  indicate  that  the 
speed  of  the  penetrator  has  virtually  no  effect  on  the  tangential  velocity  of  the  target  particles 
abutting  the  penetrator  hemispherical  nose.  However,  the  normal  stress  does  increase  with  a 
on  most  of  the  nose  surface  except  near  its  periphery  where  the  normal  stress  decreases  with  a. 
The  normal  stress  at  H  =  65°  seems  to  be  affected  little  by  the  values  of  a.  Such  a  behavior  was 
also  observed  at  ti  =  45°  for  incompressible  target  materials  [3|.  The  peak  value  of  the  normal 
stress,  which  occurs  at  the  penetrator  nose  tip,  increases  from  7.25  to  9.5  when  »  is  increased 
from  2.35  to  6. 12.  Whereas  in  the  previous  work  [3],  for  a  =  6. 13  the  target  particles  seemed  to 
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Fig.  5.  Distribution  of  the  tangential  velocity,  normal  stress  and  the  second  strain-rate  invariant  on 
the  hemispherical  nose  of  the  penetrator  for  different  values  of  a  and  compressible  target 

material. - a  =  2.63; - cr  =  3.44; - a  =  4.35; - .-  a  =  5.36; - - 

a  =6.49 

separate  away  from  the  penetrator  nose  surface  near  its  periphery,  no  such  behavior  was 
observed  in  the  present  case.  The  finite  element  mesh  used  herein  is  finer  than  the  one 
employed  earlier.  Also  the  mass  density  is  taken  as  the  independent  variable  here  whereas  in 
[3]  the  hydrostatic  pressure  p  was  taken  as  an  independent  variable  because  of  the  assumption 
of  incompressibility.  The  values  of  the  second  strain-rate  invariant  in  the  vicinity  of  the 
penetrator  nose  tip  are  affected  most  by  a  but  those  at  points  for  which  20°  s  8  s  70°  seem  to 
be  affected  less.  The  increase  in  the  value  of  the  normal  stress  on  the  penetrator  nose  with  a 
occurs  primarily  because  the  hydrostatic  pressure  p  increases  with  a. 

In  Fig.  6  is  plotted  the  effect  of  a  upon  the  variation  of  v.,  (-a..)  and  /  at  target  particles  on 
the  axial  line  ahead  of  the  penetrator.  The  penetrator  speed  does  not  affect  to  any  noticeable 
degree  the  values  of  /  and  v..  However,  the  values  of  at  points  on  the  axial  line  increase 

with  <*;  this  increase  is  mainly  due  to  the  higher  values  of  p.  For  all  values  of  or  considered  here 
the  target  particles  within  a  distance  of  almost  three  times  the  penetrator  radius  undergo 
deformations;  those  outside  this  region  hardly  deform. 

On  the  axial  line  nearly  uniaxial  strain  conditions  prevail.  Thus  the  magnitude  of  the 
deviatoric  stress  defined  as  p  +  a..,  should  equal  2/3  the  flow  stress.  For  points  on  the  axial 
line  for  which  0.0  <  :  s  1.5,  and  for  <*  =  5.36.  5;.  was  computed  to  be  0.771  and  0.741  for 
compressible  and  incompressible  materials,  respectively. 

We  note  that  the  region  studied  is  adequate  since  the  target  particles  within  a  distance  of  one 
penetrator  radius  of  the  boundary  EFA  hardly  deform  and  the  deformation  of  particles  situated 
within  one  penetrator  radius  of  the  boundary  surface  AB  are  independent  of  r. 

We  are  not  aware  of  any  experimental  work  available  in  the  open  literature  with  which  we 
could  compare  our  computations.  Nevertheless,  the  kinematic  and  stress  fields  found  in  this 
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Distance  along  the  axial  line 

Fig  6  Variation  of  the  e-velocity.  (-a.-)  and  the  second  strain-rate  invariant  /  along  the  axis  ahead 
of  the  penetrator  for  different  values  of  o  and  compressible  target  material.  Sec  Fig.  5  for 

explanation. 


study  should  help  in  devising  or  checking  the  results  from  simpler  theories  of  target 
deformations. 


CONCLUSIONS 

For  the  same  speed  of  the  rigid  cylindrical  penetrator  with  a  hemispherical  nose,  the  peak 
values  of  the  strain-rate  invariant  and  the  normal  stress,  both  of  which  occur  at  the  target 
particle  situated  at  the  penetrator  nose-tip.  are  lowered  bv  8.7 X  and  3. 19%.  respectively,  for 
compressible  target  materials  as  compared  to  that  for  incompressible  materials.  The  axial 
resisting  force  experienced  by  the  penetrator  is  nearly  10%  less  when  the  target  material  is 
compressible  as  compared  to  that  when  it  is  incompressible. 

Both  for  compressible  and  incompressible  target  materials,  significant  deformations  occur 
only  at  points  that  are  less  than  three  penetrator  radii  from  the  penetrator.  and  the  target 
seems  to  deform  farther  to  the  side  than  ahead  of  the  penetrator.  In  each  case  the  target 
material  adjacent  to  the  sides  of  the  penetrator  appears  to  extrude  rearwards  as  a  uniform 
block  that  is  separated  from  the  stationary  target  by  a  deforming  region. 

Note  added  in  proof 

Recent  unpublished  work  by  Batra  and  Liu  involving  dynamic  plane  strain  simple  compression  ol  viscoplastic  solids 
defined  by  the  constitutive  relation  (2.1)  has  shown  that,  for  short  times  after  the  application  ol  the  load,  material  lines 
in  a  direction  perpendicular  to  the  axis  of  the  compressive  load  are  also  shortened  The  constitutive  relation 

a  =  -p{p)\  +  -D  (Nl) 

•  V  3  / 
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with  /  defined  hv  eqn  (2.4)  rules  out  this  anomalous  behavior  The  constitutive  relation  (Nl)  can  be  rewritten  as 

--Im'>-$"d>'*7&“-  IN21 

Thus,  for  a  given  D.  eqns  (2.1)  and  (Nl)  give  identical  values  of  the  deviatoric  stress 

For  the  steady  state  problem  studied  herein,  constitutive  relations  (2.1)  and  (Nl)  give  essentially  the  same  results 
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University  of  Missouri-Rolla 


REFERENCES 


[1|  V  P  ALEKSEEVSKII.  Combust.  Expls  Shock  Waves  2.  63  ( 1466) 

|2)  A.  TATE.  J.  Mech  Pins.  Solids  IS.  387  (1467) 

[3|  R  C.  BATRA  and  T  W  WRIGHT.  Ini  J  Engng  Sa.  24.  41  ( I4K6) 

[4]  R  C.  BATRA.  Ini.  J  Engng  Sa.  25,  1131  (1987).' 

[5 j  R.  C  BATRA.  Comput.  Mech.  3.  1  (1988) 

|6)  P  H  PIDSLEY,  J  Mech.  Phys.  Solids  32.  315  (1984). 

[7|  L  J  HAGEMAN  and  J.  M.  WALSH.  HELP.  A  Multi-Material  Eulerian  Program  for  Compressible  Fluid  and 
Elastic -Plastic  Flows  in  Two  Space  Dimensions  and  Time  Systems  Science  and  Software  Report  3SR-350  (1971). 

[ 8 1  E  BECKER.  G.  CAREY  and  J  T  ODEN.  Finite  Elements — An  Introduction.  Vol  I  Prentice-Hall. 
Englewood  Cliffs.  N.J.  (1981) 

[9]  T  J.  R  HUGHES.  The  Finite  Element  Method.  Linear  Statu  and  Dynamic  Finite  Element  Analysis 
Prentice-Hall.  Englewood  Cliffs.  N.J.  (1987). 

|l(l|  J.  C.  HEINRICH.  P  S.  HUYAKORN  and  O  C.  ZIENKIEW1CZ.  Ini  J.  Sum  Meth  Engng  II.  131  (1977). 

( Received  29  October  1987) 


APPENDIX 


In  order  to  assess  the  validity  of  the  developed  finite  element  code  and  establish  the  reliability  of  the  computed 
results  we  solved  the  following  problem  for  compressible  Navier-Stokes  fluids  We  replaced  eqn  (2.1)  by 

a  =  -p{p)l  +  D.  p(p)  =  p/l(X)  (All 

and  the  balance  of  linear  momentum  ( 1.2)  by 

div  o  =  p(v  ■  grad  v)  +  b  (A2) 

where  h  is  the  body  force  per  unit  mass.  The  fields 

t =  r;/2.  v.=  -r;l 2.  p  =  U:  -  l/r:)/4.  (A3) 


with 


r_'_  r'r  r  I) 

2(Hl*  8  8  12’ 


+  ;/4  r 


( A4) 


he  A  I  I  he  limit  element  grid  lor  the  sample  problem 
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Tabic  A I  Comparison  of  analytical  and  numerical  solution 


Point 

Analytical  values 

Computed  values 

Difference 

r. 

-t\ 

P 

-  f. 

P 

r. 

r 

p 

P 

2.53125 

1  401,25 

II  34124 

2.5307"> 

1  4335X 

ll  342X52 

0  0|9 

1  944 

0  472 

O 

3  44531 

1  M0f>3 

0  35435 

3  45XX 

1  hXhl3 

I)  359417 

0.391 

2  774 

1  432 

R 

3  125 

1  X750 

0.52250 

3.11 4VK 

1  UI7M 

1)  52810 

o  321 

2.274 

i  o7: 

s 

2.53125 

1  X2XI3 

0  h  1077 

2  50X15 

1  X53I4 

II  1,1X034 

o  mi 3 

l  3h 8 

1  189 

1 

3  7X125 

2.57X13 

0  K45K5 

3  7643X 

2  5988  1 

0  838507 

o  314 

0X03 

ti.Hfth 

satisfy  the  balance  of  mass  and  the  balance  of  linear  momentum  Values  of  r,  and  r.  as  given  bv  the  presumed  fields 
were  assigned  on  the  boundary  faces  BC.  CD  and  DA  of  the  domain  shown  in  Fig  Al.  those  of  v,  and  o.  were 
prescribed  on  AB  and  values  ol  p  were  prescribed  on  AD  and  DC  Note  that  b,  and  b.  would  appear  on  the  right-hand 
side  of  eqn  (10  4)  and  their  values,  as  given  by  eqn  (A4).  were  inpul  into  the  computer  code  In  Tahlc  Al  are 
compared  the  values  ol  i',.  r.  and  t>  as  computed  from  expressions  (A3)  and  the  computer  code  with  i  in  eqn  <  10.5)  sel 
equal  to  0.005  ll  look  10  iterations  for  the  solution  to  converge.  The  maximum  error  in  «>,.  r.  or  p  at  any  ol  the  node 
points  was  found  to  be  less  than  3.25 
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SUMMARY 

Steady  state  penetration  of  viscoplastic  targets  by  rigid 
cylindrical  penetrators  with  an  ellipsoidal  nose  has  been  stu¬ 
died  by  the  finite  element  method.  The  target  material  is 
assumed  to  obey  a  generalized  form  of  Von-Mises  yield  crite¬ 
rion  to  account  for  strain-rate  and  strain  hardening  depen¬ 
dence.  The  ratio  of  the  major  to  minor  axes  of  the  ellip¬ 
soidal  nose  is  varied  so  chat  Che  nose  shape  corresponds  to 
essentially  a  flat  surface,  a  hemisphere  or  an  ellipsoid. 
Frictional  forces  at  the  contact  surface  are  assumed  to  be 
negligible.  Computed  results  show  that  the  nose  shape  effects 
significantly  the  deformations  of  the  target  material  in  the 
vicinity  of  the  penetrator/target  interface,  and  the  axial 
resisting  force  >xperience-l  by  the  penetrator  is  considerably 
higher  for  the  blunt  nose  as  compared  to  that  for  the  ellip¬ 
soidal  nose.  The  axial  force  depends  rather  weakly  upon  the 
speed  of  the  penetrator  for  every  nose  shape. 

INTRODUCTION 

In  an  attempt  to  shed  some  light  on  the  approximations 
made  in  simpler  theories  of  Tat  [1]  and  Alekseevski  [2], 

Batra  and  Wright  [3]  gave  a  detailed  numerical  solution, 
obtained  by  Che  finite  element  method,  to  the  following  ideal¬ 
ized  axisymmetric  penetration  problem.  They  assumed  chat  the 
penetrator  is  semi-infinite  in  length,  Che  target  is  infinite 
with  a  semi-infinite  hole,  the  rate  of  penetration  and  all 
flow  fields  are  steady  as  seen  by  an  observer  situated  on  the 
penetrator  nose,  no  shear  stress  can  be  transmitted  across  the 
target/penetrator  interface,  the  penetrator  nose  is  hemispher¬ 
ical,  and  Che  target  material  is  rigid/perfectly  plastic. 

This  problem  simulates  the  intermediate  stages  of  the  penetra¬ 
tion  process  wherein  the  active  target/penetrator  interface  is 
at  least  one  or  two  penetrator  diameters  away  from  either  tar- 


*Work  supported  by  ARO  Contract  DAAG  29-85-K-0238 
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Errata  "Effect  of  Nose  Shape  and  Strain-Hardening  on  Steady  State  Penetration 
of  Viscoplastic  Targets"  by  R.C.  Batra 
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geC  face,  Che  remaining  penecracor  is  scill  much  longer  chan 
several  diamecers  and  is  craveling  ac  a  speed  close  co  ics 
scriking  velocicy,  and  Che  CargeC  macerial  is  much  sofcer  chan 
Che  macerial  of  Che  penecracor.  Subsequencly ,  Bacra  [4]  stu¬ 
died  che  problem  when  Che  penecracor  nose  can  be  ellipsoidal 
and  che  cargec  macerial  viscoplascic.  Herein  we  excend  chis 
work  further  and  invescigate  also  che  influence  of  strain- 
hardening  of  Che  cargec  macerial. 

FORMULATION  OF  THE  PROBLEM 


Ue  use  cylindrical  coordinate  axes  attached  to  che  pene¬ 
cracor  nose  with  che  positive  2-axis  pointing  into  che  target. 
Equacions  governing  che  steady  state  deformations  of  che  tar¬ 
get  are 


div  v  ■  0, 

(1) 

div  o  *  p(ygrad)  v  . 

(2) 

Here  v  is  che  velocity  of  a  cargec  particle  as  seen  by  an 
observer  sicuaced  on  che  penecracor  note  and  moving  wich  it,  p 
is  che  current  mass  density,  o  is  che  Cauchy  stress  censor, 
and  div  and  grad  denote,  respectively,  Che  divergence  and 
gradient  operators.  Equation  (1)  expresses  che  balance  of 
mass  and  equation  (2)  Che  balance  of  linear  momentum  in  che 
absence  of  body  forces.  In  equation  (I)  it  has  been  assumed 
Chat  che  deformations  of  che  cargec  are  isochoric. 

We  neglecc  che  elastic  deformations  of  che  target  and 
assume  Che  following  constitutive  relations  for  ics  macerial. 


o  -  -p  1  +  2u(I»4>)  D,  if  D  i*  0  , 

tr(s2)  <  \  o|  (1  +  f“)2n  ,  if  D  -  0  , 
J  o  ~ 

s  -  o  +  p  1  , 

0 

2u  (I,*)  (1  +  b  I)®  (1  +  i-)"  , 

/3  I  *0 

(1  +  )n  "  Cr  * 

2I2  -  tr (D2) 

2D  ■  grad  v  +  (grad  v)^  . 


(K 

(4) 

(5) 

(6) 

(7) 

(3) 

(9) 


Here  p  is  Che  hydrostatic  pressure  that  is  not  determined  by 
che  deformation  history,  D  is  che  stretching  tensor,  oQ  '-s  the 
yield  scress  in  simple  cension  or  compression,  parameters  b 
and  m  describe  che  scrain-race  hardening  of  che  macerial,  a 
superimposed  doc  stands  for  macerial  cime  derivaclve,  and 
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parameters  4>0  and  n  characterize  the  strain  hardening  of  the 
material.  Equations  (3).  (5)  and  (6)  imply  that 

(|tr  s2)i  -  ~  o0  (1  +  bl)B  (1  +  4»/^0)n  (10) 

which  may  be  viewed  as  generalized  Von-Kises  yield  criterion 
when  the  flow  stress  (given  by  the  right-hand  side  of  (10))  at 
a  material  particle  depends  upon  its  strain-rate  and  work¬ 
hardening  parameter  <J>.  Equation  (7)  describes  the  evolution 
of  the  work-hardening  parameter  'll. 

A  constitutive  relation  similar  to  equation  (3)  has  been 
used  by  Zienklewlcz  et  al.  [5]  who  took 

2u  -  [0O  +  <21//3-r>  1/tl3//3  I  ,  (ID 

where  y  and  n  are  functions  of  the  temperature.  In  their 
study  of  the  penetration  and  perforation  problem,  Ravid  and 
Bodner  [6]  postulated  chat 

2u(I)  -  o0  (1  +  c  log10  (2I//3))//3  I  (12) 

where  c  is  a  material  constant  and  is  sec  equal  co  zero  for 
strain  rates  lower  than  unity.  We  note  chat  where  as  the 
effect  of  work  hardening  is  Included  in  equation  (6),  it  was 
not  considered  by  Zienkiewicz  et  al.  and  Ravid  and  Bodner. 

Equation  (1)  and  equations  obtained  by  substituting  from 
equations  (3),  (6),  (9)  and  (10)  into  (2)  are  the  field  equa¬ 
tions  which  together  with  appropriate  boundary  conditions  are 
co  be  solved  for  p  and  v.  Before  ::tating  the  boundary  condi¬ 
tions,  we  introduce  non-dimensional  variables  as  follows: 

a  -  a/oQ  ,  p  -  p/o0  ,  s  -  s/o0  ,  v  -  v/v0  ,  r  -  r/rQ  , 

z  -  z/r0  ,  b  -  b  v0/r0  ,  a  -  p  v^/oQ  ,  i|>  -  <l>  .  (13) 

H«re  v0  is  the  speed  of  the  penecrator  and  rQ  the  radius  of 
its  cylindrical  body. 

The  governing  equations  in  terms  of  non-dimensional  var¬ 
iables  are 


div  v  ■  0  , 

(14) 

div  o  -  a(v*grad)v  , 

(15) 

tr (oD) 

<»+♦/*)"*  (Y*8rad) 

(16) 

•'or  e 

:  ■  "P  i  +  jf-J  (1  +  bD®  (l  +  !|»/<{i0)n  D  ,  (17) 


-53- 


grad  and  div  now  sCand  for  Che  gradient  and  divergence  opera¬ 
tors  in  non-dimensional  variables  and  we  have  dropped  the  sup¬ 
erimposed  bars.  Henceforth,  unless  stated  otherwise,  we  will 
use  non-dimensional  variables  only.  We  note  that  besides  the 
ratio  of  the  major  to  minor  axes  of  the  ellipsoidal  nose  of 
the  penetrator,  there  are  five  non-dimensional  numbers  a,  b, 
m,  4»0  and  n.  All  of  these  characterize  the  target  material 
and  a  is  also  proportional  Co  the  square  of  the  penetrator 
speed. 


For  the  boundary  conditions  we  cake 


t* (an)  -  0  , 

(18) 

vn  ■  0  , 

on  Che  target/penetrator  interface,  and 

(19) 

|v  +  e|  ♦  0  as  (r2  +  z2)i  ♦  -  ,  z  >  —  , 

(20) 

| on |  +  0  as  z  ♦  -  •,  r  ^  r0  » 

(21) 

at  points  far  away  from  the  penetrator.  In  equations  (18)  - 
(21),  n  and  t  are,  respectively,  a  unit  normal  and  a  unit 
tangent  vector  on  the  bounding  surface,  and  e  is  a  unit  vector 
along  the  positive  z-axis.  The  boundary  condition  (18) 
asserts  chat  the  tangential  traction  on  the  target/penetrator 
interface  is  zero.  This  is  justified  since  a  thin  layer  of 
material  at  the  interface  either  melts  or  is  severely  degraded 
by  adiabatic  shear.  Equation  (19)  implies  no  interpenetration 
of  the  target  material  into  the  penetrator  and  vice  versa  and 
equation  (20)  states  that  the  target  material  ahead  of  the 
penetrator  back  surface  and  far  away  from  the  penetrator 
appears  to  be  moving  with  a  uniform  velocity.  That  target 
particles  behind  the  nose  and  on  the  bounding  face  perpendi¬ 
cular  to  the  axis  of  Che  penetrator  are  traction  free  is 
embodied  in  equation  (21). 

Note  that  the  governing  equations  (14)  -  (17)  are  coupled 
and  nonlinear  in  v  and  i p.  Their  solution,  if  there  is  one, 
under  the  boundary  conditions  (18)  -  (21)  will  depend,  in 
general,  upon  the  rate  at  which  quantities  in  (20)  and  (21) 
decay  to  zero.  For  a  prescribed  rate  of  decay  of  quantities 
in  equations  (20)  and  (21),  the  solutions  may  be  unique.  Also 
equation  (16)  giving  the  evolution  of  the  strain-hardening 
parameter  does  not  have  any  diffusive  term  in  it.  We  will 
gloss  over  the  rather  difficult  questions  of  existence  and 
uniqueness  and  seek  an  approximate  solution  by  the  finite 
element  method.  We  hope  Chat  this  solution  is  meaningful  for 
Che  physical  problem  at  hand. 

The  first  step  in  find  a  numerical  solution  of  the 
problem  is  to  approximate  the  infinite  target  region  by  a 
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finite  region  R  shown  in  Fig.  1.  Whether  this  region  is  suf¬ 
ficient  or  not  can  be  determined  by  solving  the  problem  for 


Fig.  1.  The  region  to  be  studied. 

two  domains,  one  larger  and  containing  the  other,  and  ascer¬ 
taining  that  the  two  sets  of  solutions  essentially  agree  with 
each  other  at  target  points  in  the  vicinity  of  the  penetrator. 
The  region  shown  in  Fig.  1  was  found  to  be  adequate  and  has 
been  used  in  results  presented  below.  The  boundary  conditions 
(18)  and  (19)  are  taken  to  hold  on  the  target/penetrator 
interface  BCD.  On  other  bounding  surfaces  following  boundary 
conditions  instead  of  equations  (20)  a.id  (21)  me  appl  ed 

vr  »  0  ,  orz  ■  0  ,  on  the  axis  of  symmetry  DE  , 

vz  «  -1 .0  ,  vr  •  0  ,  on  the  bounding  surface  EFA  ,  (22) 

vr  *  0  ,  ozz  •  0  ,  on  the  bounding  surface  BA  . 

COMPUTATION  AND  DISCUSSION  OF  RESULTS 


In  order  to  obtain  an  approximate  solution  of  the  problem, 
**  first  cast  equations  (14)  -  (17),  (18),  (19)  and  (22)  into 
ia  equivalent  weak  form  by  using  the  Galerkin  method.  The 
loader  is  referred  to  the  book  by  Becker  ec  al.  (7]  and  Che 
'aper  by  Zienkiewicz  et  al.  (5]  for  details  and  to  references 
for  a  weak  form  of  somewhat  similar  equations.  The 
■asence  of  a  diffusive  term  in  equation  (16)  necessitates  the 
se  of  either  a  superfine  mesh  or  a  reasonable  fine  mesh  with 
artificial  diffusive  term  added  to  the  equation.  Brooks 
r'd  Hughes  [8]  have  discussed  in  detail  the  justification  for 
•■’eluding  such  a  term  and  have  given  equivalent  ways  of 
thieving  the  same  objective.  We  added  a  diffusive  term  and 
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conducted  numerical  experiments  with  its  coefficient  6.  The 
two  sets  of  values  of  p,  vr  and  v2  obtained  with  6  *  10“^  and 
10~7  at  nodes  near  the  penetrator/ target  interface  differed  by 
less  than  one  percent  at  each  node.  The  results  presented 
below  are  for  6  ■  10“^. 


A  finite  element  code  employing  six  noded  isoparametric 
triangular  elements  has  been  developed.  The  incompressibility 
constraint  (14)  and  the  boundary  condition  (19)  on  the  tar- 
get/penetrator  interface  are  satisfied  by  using  the  method  of 
Lagrange  multipliers.  On  an  element,  a  quadratic  variation  of 
the  velocity  field  and  a  linear  variation  of  the  pressure 
field  defined  in  terms  of  its  values  at  the  vertices  are 
employed.  Six  quadrature  points  are  used  to  evaluate  various 
integrals  numerically,  the  governing  nonlinear  equations  are 
solved  iteratively.  The  iterative  process  is  stopped  when,  at 
each  nodal  point, 

||vi  -  yi-1||  +  1 1  <  e(||vi"1||  +  |'l>i-1|)  ,  (23) 

where  j|v|p  *  vj  +  v|,  and  e  is  a  preassigned  small  number. 
Values  of  p*  are  not  included  in  equation  (23)  since  p  appears 
linearly  in  the  governing  equations.  The  validity  of  the 
computer  code  has  been  established  by  solving  a  hypothetical 
problem  that  simulates  the  flow  of  an  incompressible  Navier- 
Scokes  fluid  in  a  circular  pipe.  The  results  for  this  sample 
problem  are  given  in  reference  3. 

Values  of  various  material  parameters  b,  m,  t|/0,  n  and  o0 
found  from  the  experimental  ~ata  in  torsion  of  Costin  et  al. 
[9]  are  assumed  to  be  valid  under  a  more  general  state  of 
stress  studied  herein.  These  and  ocher  parameters  have  been 
assigned  the  following  values  in  the  numerical  results  pre¬ 
sented  below. 


n  •  .09  ,  i |iQ  -  .017  ,  b  »  10^  sec-^ 


.025 


p  -  7800  kg/m3  ,  o0  -  180  x  106  Pa  ,  e  -  0.02  . 


(24) 


All  of  the  results  stated  below  are  in  terms  of  non-dimen¬ 
sional  variables. 


Figure  2  depicts  Che  effect  of  the  nose  shape  on  some  of 
the  variables.  For  a  relatively  blunt  nose  for  which  Che  ratio 
of  the  radius  of  the  nose  tip  to  chat  of  the  cylindrical  por¬ 
tion  of  Che  penetrator  is  taken  as  0.2,  the  normal  pressure  on 
the  nose  surface  stays  essentially  constant  and  drops  off 
sharply  to  zero  near  its  periphery.  The  abcissa  is  the  angu¬ 
lar  position  9  in  degrees.  For  rn/r0  <  1.0,  the  normal 
pressure  versus  9  curve  is  concave  downwards,  it  is  concave 
upwards  for  .rn/r0  <  1.0  and  for  values  of  between  0  and  70 
degrees.  In  every  case  the  peak  values  of  the  normal  pressure 
occurs  for  9-0*  and  the  least  value  at  9  ■  90*  where  the 
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Figure  2.  Effect  of  the  ratio  of  the  nose  radius  to  the 

radius  of  the  penetrator  on  the  normal  stress  dis¬ 
tribution  at  the  penetrator  nose  and  on  the  vari¬ 
ation  of  -ozz  and  I  at  points  of  the  target  along 
the  axis  ahead  of  the  penetrator. 

penetrator  nose  meets  its  cylindrical  part.  On  the  axial  line 
(r  «  0),  the  magnitude  of  the  principal  stress  ozz  decreases 
slowly  with  the  distance  from  the  nose  tip  for  the  blunt  nose 
as  compared  to  that  for  other  values  of  rn/rQ.  The  rate  of 
decay  of  (-ozz)  with  the  distance  from  the  nose  tip  increases 
notizesbly  with  the  increase  in  the  ratio  rn/r0.  Said 
differently,  the  deformations  spread  to  a  larger  distance  for 
a  blunt  nosed  penetrator  as  compared  to  those  for  a  penetrator 
with  a  sharp  nose.  This  is  also  confirmed  by  the  plot  of  the 
strain-rate  measure  I  on  the  axial  line.  Note  that  the  peak 
values  of  I  Increases  markedly  with  the  increase  in  the  ratio 
tn/rQ.  However,  with  an  increase  in  the  ratio  rn/r0,  I  drops 
■’ft  sharply  with  the  distance  from  the  nose  tip.  In  all  five 
cases  studied,  I  essentially  becomes  zero  at  points  distance 
-r0  from  the  nose  tip.  Thus  at  these  points  nonzero  values 
r  3zz  are  due  to  the  contribution  from  the  hydrostatic  pres¬ 
sure  p.  The  dimensional  values  of  I  typically  equal  10^  times 
::i*ir  non-dimensional  values.  Thus  strain-rates  of  the  order 
105  -  106  sec”l  occur  at  points  near  the  target/penetrator 
•ncerface . 


The  non-dimensional  axial  resisting  force  F  experienced  b% 
: ne  rigid  penetrator  is  given  by 


*/2 


L 

(n*on)  sin0  cos$  [sl"^9  +  0^)2  cos0|$  d0  , 
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where  angles  6  and  $  are  defined  in  Fig.  1.  The  corresponding 
axial  force  in  physical  unics  is  given  by  F(irro0o).  A  curve 
fit  by  the  method  of  least  squares  to  the  calculated  values  of 
F  for  different  values  of  rn/r0  gives 

r  r 

F  =  19.536  (1  -  .647  —  +  .  1409(— )2)  , 

o  o 

a  *  4.0  ,  b  =*  10^  sec-^  ,  m  =  .025  ,  n  =  .09  . 

We  note  that  an  increase  in  the  values  of  rjj/r0  from  0.2  to 
2.0  results  in  more  than  three-fold  decrease  in  the  value  of 
F. 


The  influence  of  the  values  of  strain-hardening  exponent  n 
upon  various  quantities  is  shown  in  Fig.  3  for  an  ellipsoidal 


Figure  3.  Effect  of  strain-hardening  parameter  on  the  normal 
stress  distribution  at  the  penetrator  nose  and  on 
the  variation  of  -ozz  and  I  at  points  of  the  target 
along  the  axis  ahead  of  the  penetrator. 

nosed  penetrator.  The  normal  pressure  on  the  penetrator  nose 
increases  with  n  and  so  does  the  magnitude  of  the  principal 
stress  ozz  on  the  axial  line.  However,  the  values  of 
strain-rate  measure  I  at  points  on  the  axial  line  seem  to  be 
affected  a  little  by  the  value  of  n.  A  curve  fit  to  the  com¬ 
puted  values  of  the  axial  resisting  force  F  gives 

F  =  4.525(1  +  l . 793n)  , 

a  »  4  ,  rn/r0  ■  2  ,  b  ■  10^  sec-^  ,  ra  *  0.25  , 
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implying  thereby  that  F  depends  linearly  upon  n. 

Results  plotted  in  Fig.  4  illustrate  the  effect  of  the 
speed  upon  various  for  a  blunt  nose  penetrator.  Recall  that 


Figure  4.  Effect  of  speed  of  the  penetrator  on  the  normal 

stress  distribution  at  the  penetracor  nose  and  on 
the  variation  of  -ozz  and  I  at  points  of  the  targec 
along  the  axis  ahead  of  the  penetrator. 

c.ie  non-dimensional  number  a  is  proportio *nl  to  chr  square  of 
ciio  speed.  The  normal  pressure  on  the  penetrator  nose 
increases  with  a.  Note  that  for  lower  speeds,  the  normal 
pressure  near  Che  extremity  of  Che  penetrator  nose  increases 
■'’-':'>re  dropping  off  sharply  to  zero.  For  higher  speeds  this 
Phenomenon  is  absent  and  Che  normal  pressure  is  uniformly 
-i^tributed  on  most  of  the  nose  surface  and  drops  off  rapidly 
eir  the  edges.  At  a  point  on  the  axial  line,  the  magnitude 
che  principal  stress  oz2  increases  with  the  speed.  At 
->ir  speeds,  the  peak  value  of  I  at  points  on  che  axial  line 
not  occur  at  the  nose  tip  but  at  a  poinc  whose  distance 
i  che  nose  tip  is  nearly  equal  to  che  radius  of  che  pene- 
r'Cor.  At  higher  speeds,  che  strain-rate  I  becomes  essen- 
•  >Uy  constant  over  a  distance  equal  to  half  che  radius  of 
-  penetracor  and  then  gradually  drops  off  to  zero.  In 
-her  case  I  approaches  zero  at  points  on  che  axial  line 
’**■’  distance  from  che  nose  tip  is  2.5  r0  or  more.  The 
--vndence  of  che  axial  resisting  force  F  upon  a  is  given  by 

'  1^.641(1  +■  .04211a)  ,  rn/r0  *  .2  ,  b  »  10^  sec” l, 

31  -  -025  .  n  -  .09. 
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Thus  F  depends  upon  a  rather  weakly.  This  provides  a  justi¬ 
fication  for  the  assumption  of  the  constant  target  resistance 
in  the  simple  theory  of  Tate  [ 1 ] . 

In  Fig.  S  is  plotted  the  effect  of  speed  on  target  defor¬ 
mations  when  the  penetrator  nose  is  ellipsoidal  rather  than 


WT!0»  •  (0*C*?S 


Figure  5.  Effect  of  the  penetrator  speed  on  the  normal  stress 
distribution  at  the  penetrator  nose  and  on  the  var¬ 
iation  of  -ozz  and  I  at  points  of  the  target  along 
the  axis  ahead  of  the  penetrator. 

hemispherical  nosed  penetrator  in  Ref.  3,  the  value  of  the 
normal  pressure  on  the  penetrator  nose  at  6  equal  to  approxi¬ 
mately  forty-five  degrees  is  unaffected  by  the  penetrator 
speed.  At  points  on  the  penetrator  nose  for  which  6  is  less 
than  45°,  the  normal  pressure  increases  with  a  and  on  the 
remainder  of  the  surface,  it  decreases  with  a.  However,  the 
relative  decrease  is  smaller  than  the  percentage  increase  at 
points  situated  equal  angular  distance  from  the  6  »  43°  line. 
At  a  point  on  the  axial  line  ahead  of  the  penetrator,  whereas 
the  magnitude  of  azz  increases  with  that  of  I  essentially 
remains  unchanged  except  at  a  point  on  the  nose  tip  where  it 
increases  slightly  with  a.  The  rate  of  decay  of  (-ozz)  and  I 
with  Che  distance  from  Che  nose  tip  is  essentially  unchanged 
as  Che  penetrator  speed  increases.  Values  of  I  essentially 
approach  zero  at  points  on  the  axial  line  whose  distance  from 
the  penetrator  nose  tip  is  at  least  1.5  r0. 

A  curve  fit  Co  the  computed  values  of  the  axial  resistance 
force  F  at  various  values  of  a  gives 

F  -  4.955(1  +  .0151a)  , 
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implying,  as  before,  chat  F  depends  weakly  upon  a. 

Figures  6  and  7  compare,  respectively,  the  variation  of  I 
and  in  the  target  region  considered  for  two  different  nose 
shapes  of  the  pene tracer.  When  the  penetrator  nose  is  ellip¬ 
soidal,  peak  values  of  I  occur  at  Che  nose  tip  and  I  decreases 
as  one  moves  away  from  it  in  any  direction.  However,  for  a 
blunt  nosed  penetrator  peak  values  of  I  occur  at  points  on  the 
targec/penetrator  interface  where  Che  nose  meets  the  cylin¬ 
drical  portion  of  the  penetrator.  The  peak  value  of  I  at  any 


figure  6.  Comparison  of  the  variation  of  the  strain-rate 

measure  I  in  the  targec  region  for  an  ellipsoidal 
nose  and  a  blunt  nose. 

point  in  the  target  is  nearly  three  times  larger  when  the 
Penetrator  has  a  blunt  nose  as  compared  to  that  when  it  has  a 
•parp  ellipsoidal  nose.  The  strain-hardening  parameter  exhi- 
'-Cs  a  totally  different  behavior.  For  a  blunc-nosed  pene- 
■  ■".icor  che  peak  value  of  1 1>  occurs  at  che  nose  cip  and  it  drops 
'■t  rather  sharply  as  one  moves  away  from  it.  It  stays  prac- 
:  -ca Lly  unchanged  on  most  of  the  penetrator/targec  interface, 
or  an  ellipsoidal  nosed  penetracor  the  peak  value  of  ^  does 
ot  occur  at  che  nose  cip  but  at  a  point  on  che  penetra- 
•r'carget  interface  which  is  on  the  penetracor  nose  and  for 
-Pich  9  =  l”.  From  this  point  on  the  value  of  stays  prac- 
‘--cally  che  same  on  the  remainder  of  che  cargec/penecrator 
■>ter£ace.  Again  Che  peak  value  of  at  any  target  point  for 
pfunc  nosed  penecracor  is  roughly  two  and  a  half  times  as 
:r^e  as  chat  for  an  ellipsoidal  nosed  penetrator.  A  possible 
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Figure  7.  Comparison  of  the  variation  of  the  strain-hardening 
parameter  <//  in  the  target  region  for  an  ellipsoidal 
nose  and  a  blunt  nose. 

reason  for  the  difference  in  the  variations  of  and  I  in  the 
two  cases  is  that  the  governing  equation  for  is  convection 
dominated  whereas  that  for  I  does  have  a  natural  diffusive 
term  in  it.  As  depicted  in  Fig.  6  of  Ref.  3,  the  tangential 
velocity  on  the  penetrator  nose  is  quite  different  in  the  two 
cases.  Whereas  for  a  blunt  nosed  penetrator  it  increases  from 
zero  at  first  slowly  and  then  quite  rapidly  when  9  is  close  to 
80°,  for  the  ellipsoidal  nosed  penetrator  it  increases  rapidly 
for  0  _<  0  <  18®  and  gradually  afterwards. 

Finally  we  note  that  an  extension  of  this  work  to  the  the 
thermo-mechanical  problem  when  the  target  material  can  also 
thermally  soften  will  appear  in  Ref.  10. 

CONCLUSIONS 

During  the  intermediate  stages  of  penetration  of  a  target 
all  of  whose  dimensions  are  very  large  as  compared  to  the 
penetrator  radius  and  whose  material  is  much  softer  as  com¬ 
pared  to  that  of  the  penetrator,  the  shape  of  the  penetrator 
nose  has  most  influence  upon  target  deformations.  For  a  blunt 
nosed  penetrator  the  deformations  spread  to  a  larger  distance 
away  from  the  target/penetrator  interface  as  compared  to  chat 
for  an  ellipsoidal  nosed  penetrator.  In  every  case  Che  axial 
resisting  force  experienced  by  the  penetrator  depends  weakly 
upon  the  penetrator  speed.  Also  the  deformation  related  field 
variables  depend  smoothly  upon  the  strain  hardening  exponenc. 
The  peak  values  of  strain-rate  measure  I  and  strain-hardening 
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parameter  i>  are  higher  for  a  blunt  nosed  penetrator  as  com¬ 
pared  to  those  for  an  ellipsoidal  nosed  penetrator.  Points 
where  these  peak  values  occur  are  situated  quite  differently 
in  the  two  cases. 
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Abstract  Steady  state  thermomechanical  deformations  of  a  target  hit  by  a  rigid  cylindrical  penetrator  with  an  ellipsoidal 
nose  are  studied.  The  material  of  the  target  is  assumed  to  be  thermally  softening  but  strain  and  strain-rate  hardening.  Results 
computed  and  presented  graphically  include  the  pressure  distribution  on  the  nose  of  the  penetrator,  dependence  of  the  axial 
resisting  force  upon  the  speed  of  the  penetrator,  and  the  variation  of  field  quantities  such  as  the  temperature  and  strain-rate 
in  the  target.  Computed  results  show  that  the  ratio  of  the  major  to  minor  axes  of  the  ellipsoidal  nose  has  a  significant  effect 
on  the  deformations  of  the  target  particles  in  the  vicinity  of  the  penetrator  nose. 


1  Introduction 

In  an  attempt  to  shed  some  light  on  the  validity  of  the  approximations  made  in  simple  theories  of 
penetration  due  to  Alekseevskii  (1966)  and  Tate  (1967,  1969),  Batra  and  Wright  (1986)  studied  in 
detail  the  penetration  problem  that  simulates  the  following  situation.  Suppose  that  the  penetrator 
is  in  the  intermediate  stages  of  penetration  so  that  the  active  target/penetrator  interface  is  at  least 
one  or  two  penetrator  diameters  away  from  either  target  face,  and  the  remaining  penetrator  is  much 
longer  than  several  diameters  and  is  still  traveling  at  a  uniform  speed.  Thus  steady  state  deformations 
of  the  target,  presumed  to  be  made  of  a  rigid/perfectly  plastic  material,  and  being  penetrated  by  a 
long  cylindrical  rigid  rod  with  a  hemispherical  nose  were  analysed.  Subsequently  Batra  (1987) 
showed  that  the  axial  resisting  force  experienced  by  the  rigid  penetrator  is  considerably  reduced  if 
its  nose  shape  is  ellipsoidal  rather  than  hemispherical  and  also  investigated  the  effect  of  the  depend¬ 
ence  of  the  flow  stress  upon  the  strain-rate.  Herein  we  study  the  steady  state  penetration  problem 
when  the  target  material  is  thermally  softening  but  strain  and  strain-rate  hardening. 

Pidsley  (1984)  has  recently  given  a  detailed  numerical  solution  of  the  penetration  problem  in 
which  both  materials  are  considered  to  be  deformable  and  rigid/perfectiy  plastic.  We  refer  the  reader 
to  his  paper  for  more  references  on  this  subject.  Even  though  we  study  a  somewhat  simpler  situation, 
our  material  model  is  more  general  in  that  we  account  for  the  effect  of  strain  and  strain-rate 
hardening  and  thermal  softening.  We  note  that  the  peak  strains  and  strain-rates  encountered  during 
steady  state  deformations  of  the  target  are  of  the  order  of  10  and  105  sec-1  respectively.  Also  the 
temperature  at  target  points  may  rise  to  as  much  as  half  of  the  melting  point  of  the  target  material. 
We  study  the  effect  of  these  competing  factors  as  well  as  of  the  penetrator  speed  and  the  shape  of 
its  nose  on  the  deformations  of  the  target. 


2  Formulation  of  the  problem 

Since  the  axisymmetric  deformations  of  the  target  appear  to  be  independent  of  time  to  an  observer 
situated  on  the  penetrator  nose  and  moving  with  it,  we  choose  a  cylindrical  co-ordinate  system 
attached  to  the  nose  tip  with  the  positive  z-axis  pointing  into  the  target  material.  With  respect  to 
these  axes  translating  with  a  uniform  velocity  v0e,e  being  a  unit  vector  along  the  penetrator  axis 
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and  in  the  direction  of  its  motion,  equations  governing  the  thermomechanical  deformations  of  the 


target  are 

dive  =  0,  (l) 

dive  =  grad)®,  (2) 

—  divf  +  tr(<rD)  =  <?(»  •  grad)  U,  (3) 

tr(oD)  =  (o  •  grad)x,  (4) 

D  =  (grad  c  +  (grad  v)^fl .  (5) 


Equations  (1)  through  (4)  written  in  Eulerian  description  of  motion  express  respectively  the  balance 
of  mass,  linear  momentum,  internal  energy  and  the  evolution  of  the  work  hardening  parameter  x. 
In  Eq.  (4)  we  have  neglected  the  elastic  deformations  of  the  target  and  in  Eq.  (3)  assumed  that  all 
of  the  plastic  working,  rather  than  a  part  of  it  (e.g.  Farren  and  Taylor  192S)  is  dissipated  in  the  form 
of  heat.  The  operators  grad  and  div  signify  the  gradient  and  divergence  operators  on  fields  defined 
in  the  present  configuration.  Furthermore,  <r  is  the  Cauchy  stress  tensor,  q  is  the  mass  density  of  the 
target  material,  v  is  the  velocity  of  the  target  particle  relative  to  the  penetrator,  q  is  the  heat  flux  per 
unit  present  area,  D  is  the  strain-rate  tensor,  and  U  the  specific  internal  energy.  Equations  (1) 
through  (4)  are  to  be  supplemented  by  constitutive  relations  and  boundary  conditions. 

We  assume  the  following  constitutive  relations  for  the  target  material. 


q  =  --t  grad0  , 

(6) 

U=c6, 

(7) 

a  =  —  p\  +  2p{I,9,  ip)D,  ifD#0, 

(8) 

tr(s2)  <  ^(1  -a0)2(  1  +  — )  ,  ifD  =  0, 

3  V  Vo/ 

(9) 

S  =  <F  +/7  1, 

(10) 

2p(I,0,y>)  =  j^(l  +  bini  -  off)  (l  +  ^)", 

(ID 

x(ip)  =  <r0ip(\  +ip/\pQy, 

(12) 

/2  =  }tr(D2). 

(13) 

Equation  (6)  is  Fourier’s  law  of  heat  conduction,  k  is  the  thermal  conductivity,  9  is  the  change  in 
the  temperature  of  a  material  particle  from  that  in  the  underformed  configuration,  c  is  the  specific 
heat  which  is  assumed  to  be  constant,  p  is  the  hydrostatic  pressure  not  determined  by  the  deformation 
history,  and  <r0  >s  the  yield  stress  in  a  simple  tension  or  compression  test.  The  material  parameters 
b  and  m  describe  the  strain-rate  sensitivity  of  the  material,  the  material  parameter  a  describes  its 
thermal  softening,  and  %  and  n  characterize  the  strain  hardening  of  the  material.  An  integral  form 
of  Eq.  (12)  with  x  interpreted  as  the  true  stress  and  ip  the  plastic  strain  represents  the  stress-strain 
curve  in  a  quasistatic  reference  test.  Equation  (8)  may  be  interpreted  as  a  constitutive  relation  for 
a  Non-Newtonian  fluid  whose  viscosity  p  depends  upon  the  strain-rate,  temperature  and  a  material 
parameter  ip.  Equation  (8)  implies  that 

G^)*_i75'*(l+*'r(,-*')(,  +  s)'  (,4) 

which  can  be  viewed  as  a  generalized  Von-Mises  yield  criterion  when  the  flow  stress  (given  by  the 
right-and  side  of  (14))  at  a  material  particle  depends  upon  its  strain-rate,  strain  and  temperature.  A 
constitutive  relation  similar  to  Eq.  (8)  has  been  used  by  Zienkiewicz  et  al.  (1981)  who  took 

2p=[a0+(2II\fi7yi']l]/3I,  (15) 
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where  y  and  n  are  functions  of  0.  They  asserted  that  it  represents  Perzyna’s  viscoplastic  model.  For 
a  simple  shearing  deformation,  Litonski  (1959)  proposed  that 

r  =  c(l  —  a0)(l-+  byY'y"  (16) 

where  x  and  y  equal  the  shear  stress  and  shear  strain,  and  c  is  a  material  constant.  Note  that  this 
relation  implies  that  t  is  zero  whenever  y  =  0.  Another  stress-strain  law  proposed  by  Lindholm  and 
Johnson  (1984),  based  on  fitting  curves  to  experimental  data  obtained  from  torsion  tests,  is 

t  =  (A  +  £y")(l  +  C  ln(y/y0))  ~~  (17) 

“m  ~  “0 

where  9m  is  the  melting  point  of  the  material,  90  is  a  reference  temperature,  and  A,  B  and  C  are 
material  constants.  Lin  and  Wagoner  (1986)  recently  reported  that  the  following  curve 

a  =  556 (e  -  0.014)°-2,9(e/0.02)0018 ( 1  -  0.0012(0  -  298))  MPa  ,  (18) 

fitted  well  their  experimental  data  derived  from  a  uniaxial  tension  test  on  Armco  I.  F.  steel.  In  Eq. 
(18),  a  and  e  are  the  axial  stress  and  the  axial  strain  respectively  and  0  is  in  °K.  The  linear  dependence 
of  the  flow  stress  upon  temperature  has  also  been  observed  by  Bell  (1968). 

The  constitutive  relation  (8)  with  p  given  by  Eq.  (1 1)  is  an  attempt  to  generalize  the  one  used  by 
Wright  and  Batra  (1986)  for  simple  shearing  deformations  of  nonpolar  and  dipolar  materials.  They 
used  it  to  study  shear  strain  localization  phenomenon  in  metals  and  derived  it  by  using  arguments 
similar  to  those  employed  by  Green,  Mclnnis  and  Naghdi  (1968).  A  curve  fit  to  the  torsion  test 
data  of  Costin  et  al.  (1979)  for  a  1018  cold  rolled  steel  gives  n  =  0.09,  rp0  —  0.017,  b—  104  sec-1  and 
m  =  0.025. 

Before  stating  the  boundary  conditions  we  non-dimensionalize  the  variables  as  follows. 


5  =  a/o0,  p  =  p/a0,  s  =  s/a0,  r  =  o/u0,  f  =  r/r0,  f  =  z/r0,  6  =  0/0„  , 

„  09) 

B  —  b—,  a  —  ad0,  ci  =  qvI)o0,  $  -  kj(ecv0r0),  90-cr0l(Qc). 

r0 

Substituting  from  Eqs.  (6)  through  (12)  into  the  balance  laws  (1)  through  (4),  rewriting  these  in 
terms  of  non-dimensional  variables,  and  denoting  the  gradient  and  and'divergence  operators  in  non- 
dimensional  coordinates  by  grad  and  div,  we  arrive  at  the  following  set  of  equations. 


div  v  =  0, 

(20) 

div®  =  a(v  ■  grad)». 

(21) 

tr(®D)  +  div  (grad  0)  =  (p  ■  grad)0  , 

(22) 

tr(®D) 

„  —7  -t;  =  (o-grad)v, 

(I  +  rp/rpoY 

(23) 

where 

a  =  —pl  +  -r*~(I  +  6/)"(I  —a0)(l  +  —  j  D, 

J/37  \  WoJ 

(24) 

and  we  have  dropped  the  superimposed  bars. 

We  assume  smooth  contact  at  the  target/penetrator  interface.  Thus  the  boundary  conditions  on 
this  surface  are 


/  •  (an)  =  0 , 
v  n  —  0, 
q  ■  n=>  h(9  —  9a), 


(25.1) 

(25.2) 

(25.3) 


where  h  is  the  heat  transfer  coefficient  between  the  penetrator  and  the  target,  0O  is  an  average 
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temperature  of  the  penetrator,  and  n  and  t  are,  respectively,  a  unit  normal  and  a  unit  tangent  vector 
to  the  surface.  At  points  far  away  from  the  penetrator 

|®  +  *|-»0,  0 -*  0,  tp  -*  0  as  (r2  +  z2)1  -*  oo,  z>— oo,  (26.1) 

|o«|-»0,  |f*a|  -»  0,  tp  — ♦  0  as  z  -*  — oo,  r^r0.  (26.2) 

The  boundary  condition  (26.1)  states  that  target  particles  at  a  large  distance  from  the  penetrator 
appear  to  be  moving  at  a  uniform  velocity  with  respect  to  it  and  experience  no  change  in  their 
temperature.  Equation  (26.2)  implies  that  far  to  the  rear  the  traction  and  heat  flux  fields  vanish. 

Note  that  the  governing  Eqs.  (20)— (23)  are  coupled  and  are  nonlinear  in  c,  0  and  tp.  Their 
solution,  if  there  is  one,  may  not  be  unique  and  will  depend,  in  general,  upon  the  rates  at  which 
quantities  in  (26)  decay  to  zero.  Since  we  are  unable  to  solve  these  equations  analytically  and  prove 
any  uniqueness  theorem,  we  will  seek  a  numerical  solution  of  these  equations  which  we  hope  will 
be  physically  meaningful. 

3  Finite  element  formulation  of  the  problem 

The  numerical  solution  of  the  problem  necessitates  the  consideration  of  a  finite  region.  Since  the 
target  deformations  are  axisymmetric,  only  the  target  region  R  shown  in  Fig.  1  is  studied.  The 
adequacy  of  the  finite  domain  considered  will  be  verified  by  solving  the  problem  for  two  separate 
regions,  one  larger  and  containing  the  other,  and  ensuring  that  the  two  sets  of  computed  values  of 
various  field  quantities  are  close  to  each  other.  The  boundary  conditions  (26)  are  replaced  by  the 
following. 

on  the  surface (27.1) 


on  the  axis  of  symmetry  DE,  (27.2) 

vr  —  0,  vz=-1.0,  0  =  0,  ip  -0  on  the  boundary  surface  EFA.  (27.3) 

Referring  the  reader  to  Becker  et  al.  (1981)  and  Zienkiewicz  et  al.  (1981)  for  details,  we  simply  note 
that  a  weak  formulation  of  the  problem  defined  on  the  region  R  by  Eqs.  (20)-(24)  and  boundary 
conditions  (25)  and  (27)  is  that  equations 

J  A  (div  r)  d»  =  0,  (28.1) 

R 


A  _  vv 

<r„  =  0,  vr  —  0,  —  =  0 
or 

30 

arI  =  0,  vr  =  0,  —  =  0 
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J/>(div^)du  -  j/t(/,0,v>)(D:(grad^  +  (grady)r))di;  =  a  J ((»  grad)  p)  <pdr, 


(28.2) 


J  ((*•  grad)  0)  i;  da  + /Jj  grad  0- grad  da +  h  J  drjdA  =  j»f(<r:D)du—  J  hBai\dA , 
r  r  a  ,R  R  i,R 


I ((a  •  grad)y)^dt>  =  }  -  £de, 

*  '  *(l+*Y 

\  Vo/ 


(28.3) 


(28.4) 


hold  for  arbitrary  smooth  functions  2,  <p,  r\  and  {  defined  on  R  such  that  <p,  =  0  on  AB ,  <p  —  0  on 
EFA,  q>  •  n  =  0  on  the  target/penetrator  interface  BCD ,  and  rj  =  q  =  0  on  EFA.  In  these  equations, 
A :  B  =  tr(A  Br)  for  linear  transformations  A  and  B,  and  6,  R  denotes  the  target/penetrator  interface 
BCD.  Since  these  equations  are  nonlinear  in  a,  0  and  tp,  the  following  iterative  technique  has  been 
employed.  At  the  /'*  iteration,  equations 


f  2  diva' da  =  0, 


(29.1) 


JV(div<*)do  =  f /i (/'"*, 9‘~l, ip*"1) (D*: (grad +  (grad<<»)r))da  =  aj(k  1  •  grad) v")  <pdv,  (29.2) 

R  R 

J ((o'-1  •  grad )6‘)rfdv  +  /JjgradB1  •  gradrf  da  +  h  |  9‘rjdA  ~  1 n(o'~l  :D,_l)da  -  j  hOar]dA,  (29.3) 

R  d,R  R  3,/f 


J((r'-1  •grad)v>'Kdp  =  J-7^ — -ff.—  {da,  (29.4) 

*  *(i  +  sii!  Y 

\  Vo  / 

are  solved  for  a*,  &,  tp*  and  pi.  The  iterative  process  is  stopped  when,  at  each  nodal  point, 

Ik  — *'-‘11  +  1 9‘  -  0"'|  +  IV*  -  V'-'I  <  e [Ik* 'll  +  lO'-'l  +  Iv'-'ll  (30) 

where  ||  o||2  =  a2  +  a2,  and  c  is  a  preassigned  small  number.  Values  of  p‘  are  not  included  in  Eq.  (30) 
since  p  appears  linearly  in  Eq.  (29.2). 


4  Computation  and  discussion  of  results 

The  finite  element  code  developed  earlier  (Batra  and  Wright  (1986)]  to  solve  the  problem  when  the 
target  material  is  modeled  as  rigid/perfectly  plastic  and  the  penetrator  nose  is  hemispherical  has 
been  modified  to  solve  the  present  problem.  It  employes  six-noded  triangular  elements  with  a,,  a., 
\p  and  6  approximated  by  quadratic  functions  over  an  element  and  p  by  a  linear  function  defined  in 
terms  of  its  values  at  the  vertices  of  the  triangular  element.  The  validity  of  the  code  was  established 
by  first  modifying  Eqs.  (29)  to  include  arbitrary  source  terms  on  their  right-hand  sides,  and  then 
solving  simple  problems  for  incompressible  Navier-Stokes-Fourier  fluids.  The  source  terms  were 
adjusted  so  that  the  governing  equations  were  satisfied  by  the  presumed  analytic  expressions  for  ar, 
vf,  p,  0  and  ip.  Results  for  a  sample  problem  that  does  not  include  thermal  effects  are  given  in  Batra 
and  Wright  (1986). 

A  major  difference  between  the  problem  studied  herein  and  those  studied  earlier  by  Batra  and 
Wright  (1986)  and  Batra  (1987)  is  that  Eq.  (28.4)  does  not  have  any  diffusive  term  in  it.  This 
necessitates  the  use  of  either  an  ultrafine  mesh  or  a  fine  mesh  with  an  artificial  diffusive  term  included 
in  Eq.  (28.4).  Brooks  and  Hughes  (1982)  have  discussed  in  detail  the  justification  for  including  such 
a  term  and  have  given  equivalent  ways  of  achieving  the  same  objective.  We  added  a  term 

5  j  grad  ip  •  grad  {da  to  the  left-hand  side  of  Eq.  (28.4)  and  computed  results  for  <5  =  IO-6and  10-7. 

R 

The  two  sets  of  values  of  0,  p,  vr  and  vz  differed  by  less  than  one  percent  at  each  node.  The  results 
presented  below  are  for  5  —  10-6.  We  next  ascertained  the  adequacy  of  the  region  considered  by 
increasing  DE  in  Fig.  I  from  3.25 r0  to  4.5 r0.  Again  the  difference  in  the  values  of  6,  p,  vr,  v2  and  tp 
at  points  in  the  vicinity  of  the  penetrator  nose  was  negligible. 
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We  note  that  experimental  data  for  the  range  of  deformations  expected  to  occur  in  the  penetration 
problem  under  study  is  not  available  in  the  open  literature.  Therefore,  values  of  material  parameters 
b,  m,  a,  ip0  and  n  in  Eq.  (14)  found  by  fitting  a  curve  to  the  experimental  data  in  torsion  of  Costin 
et  al.  (1979)  were  assumed  to  be  valid  under  the  more  general  state  of  stress  studied  here.  This 
should  enable  us  to  undertake  the  parametric  study  for  a  reasonable  range  of  values  of  various 
material  parameters.  The  values  of  various  parameters  used  to  compute  numerical  results  are: 

n  =  0.09,  tp0  =  0.017,  b=  \0*scc  ,  m  =  0.025,  a  =  0.000555/°C,  k  =  48W/m°C, 
c  =  473  J/kg  °C,  q  =  7800  kg/m3,  a0  =  1 80  x  106  Pa,  h  =  20  W/m2  °C,  r0  =  2.54  mm, 
e  =  0.02,  0,  =  0. 

However,  the  results  presented  below  are  in  terms  of  non-dimensional  quantities  and  the  variables 
that  are  assigned  values  different  from  those  given  above  are  so  indicated  in  the  figures  along  with 
their  new  values. 

In  Fig.  2  is  plotted  the  pressure  distribution  on  the  penetrator  nose  for  a  relatively  blunt  nose 
(rjr o  =  0.2),  a  hemispherical  nose  (r„/r0  =  1)  and  an  ellipsoidal  nose  (r„/r0  =  2.0).  As  expected  the 
normal  pressure  on  the  blunt  nose  stays  essentially  uniform  over  most  of  its  surface  and  drops  off 
sharply  near  its  extremities.  Note  the  change  in  the  curvature  of  the  pressure  curve  in  going  from 
hemispherical  to  an  ellipsoidal  nose.  The  non-dimensional  axial  resisting  force  decreased  from 
17.091  for  the  blunt  nose  to  8.902  for  the  hemispherical  nose  and  further  to  5.085  for  the  ellipsoidal 
nose.  The  axial  resisting  force  F  is  given  by 

F=  2  \  (n  •  <rn) sin 0 cos <p  £sin20  +  ^  cos20  J^d0, 

where  angles  0  and  <p  are  defined  in  Fig.  1.  The  corresponding  axial  force  in  physical  units  is  given 
by  F(nr\a^).  The  normal  pressure  on  the  hemispherical  and  the  ellipsoidal  nose  surface  for  the 


F1«.Z  Distribution  of  the  normal  pressure  on  the  peneirator 

nose  for  three  different  nose  shapes.  -  Blunt  nose 

(r,/r0  =  0.2);  - hemispherical  nose  (r„/r0  =  1.0); 

. ellipsoidal  nose  (rjr0  =  2.0) 


Ftj.  3.  Variation  of  the  strain-rate  measure  /  and  the  temper¬ 
ature  change  6  on  the  penetrator  nose  for  three  different  nose 
shapes.  - Blunt  nose  (rjr0  =  0.2);  . hemispher¬ 
ical  nose  (r./r0  =  1 .0); - ellipsoidal  nose  (r,/ra  =  2.0); 

a  =  4.0 
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angle  0  greater  than  nearly  75°  is  not  plotted  because  of  the  difficulties  encountered  in  computing 
it  accurately.  The  mesh  in  this  region  was  not  fine  enough  to  yield  reliable  values.  Figure  3 
depicts  the  variation  of  the  strain-rate  invariant  /  and  temperature  change  6  on  the  nose  surface. 
Whereas  the  maximum  value  4.21  of  /  occurs  at  the  penetrator  nose  tip  for  the  ellipsoidal  nose,  it 
assumed  very  high  values  at  the  extremities  for  the  blunt  nose.  For  these  two  nose  shapes  significant 
values  of  /  occur  near  the  nose  tip  and  the  nose  periphery  respectively.  For  the  hemispherical  nose 
shape  I  decreases  almost  linearly  from  its  maximum  value  of  2.16  at  the  nose  tip  to  0.4  at  its 
periphery  (0  —  90°).  The  dimensional  values  of  /  equal  1.1  (10s)  times  the  non-dimensional  values. 
The  values  of  temperature  at  the  nose  tip  do  not  depend  that  much  on  the  nose  shape.  However  the 
temperature  decreases  with  6  for  the  ellipsoidal  and  the  hemispherical  nose,  it  increases  with  0  for 
the  blunt  nose.  Because  of  the  high-strain  rates  near  the  vicinity  of  the  periphery  of  the  blunt  nose, 
there  is  a  lot  of  heat  generated  in  this  narrow  region.  Since  material  particles  near  the  periphery  of 
the  nose  are  moving  relatively  slowly,  not  much  of  the  heat  produced  is  convected  or  transported 
away.  In  Fig.  4  is  plotted  the  variation  of  the  strain  rate  /  and  temperature  change  9  on  the  axial 
line.  For  the  blunt  nose,  the  deformation  has  spread  to  a  larger  distance  as  compared  to  that  for 
the  ellipsoidal  nose.  Accordingly  the  temperature  drops  off  slowly  for  the  blunt  nose  than  it  does 
for  the  other  two  cases.  The  actual  temperatures  in  °C  are  obtained  by  multiplying  their  non- 
dimensional  values  by  48.9.  Thus  temperatures  as  high  as  605  °C  occur  at  and  near  the  nose  tip. 
The  maximum  value  of  strain-rate  /  on  the  axial  line  appears  to  occur  at  a  point  slightly  away  from 
the  nose  tip.  This  initial  rise  is  probably  only  an  artifact  and  the  curves  should  be  extended  smoothly 
to  the  nose-tip. 


Fig.  4.  Variation  of  the  temperature  change  and  strain-rate  measure  I  on  the  axial  line  for  three  different  nose  shapes. 
- Blunt  nose  (r„/r0  =  0.2); . hemispherical  nose  (r«/r0  =  1.0); - ellipsoidal  nose  (r,/r„  =  2.0);  a  ~  4.0 


Fig.  5.  Pressure  and  temperature  distribution  on  the 
ellipsoidal  penetrator  nose  for  different  values  of  a. 

- a  «■  1.0; - a  «  2.0; - a  =  3.0; 

. a  —  4.0; - a  «  5.0 
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Figures  5  and  6  illustrate  the  effect  of  speed  of  the  penetrator  on  various  solution  parameters  at 
or  in  the  vicinity  of  its  ellipsoidal  nose  with  r„/r0  =  2.0.  As  shown  in  Fig.  5,  the  normal  pressure 
near  the  penetrator  n*jse  tip  increases  with  the  speed  but  decreases  near  its  periphery.  Near  9  —  45°, 
the  speed  has  no  effect  on  the  normal  pressure.  Such  a  behavior  was  also  observed  for  a  hemispherical 
nose  and  a  rigid/perfectly  plastic  target  material  by  Batra  and  Wright  (1986).  The  dependence  of 
the  non-dimensional  axial  force  F  upon  the  speed  (through  non-dimensional  variable  a) 
is  given  by 

F  =  5.021  +  0.0732  a,  ellipsoidal  nose  (r„/r0  =  2.0) 

F=  8.71  +  0.2145a,  hemispherical  nose. 

For  rigid/perfectly  plastic  materials,  Batra  and  Wright  (1986)  obtained  F -  3.903  +  0.0773a  for  a 
hemispherical  nosed  penetrator,  and  Batra  (1987)  computed  F=  2.58  4-  0.019a  for  a  penetrator 
with  an  ellipsoidal  nose  having  rjr0  =  2.0.  Thus  the  consideration  of  strain  and  strain-rate  hardening 
and  thermal  softening  effects  more  than  doubles  the  axial  resisting  force.  In  every  case  studied  so 
far,  F  depends  upon  a  weakly.  This  weak  dependence  of  Fupon  a  seems  to  explain  why  the  choice 
of  constant  target  resistance  in  the  simple  theory  of  Tate  (1967,  69)  gives  such  good  qualitative 
results.  On  most  of  the  nose  surface,  the  temperature  decreases  with  a.  This  is  shown  in  Fig.  5. 
Figure  6  depicts  that  most  of  the  target  deformations  are  concentrated  near  the  penetrator  nose. 
The  peak  value  of  I  on  the  axial  line  appears  to  occur  not  at  the  nose  tip  but  slightly  away  from  it. 
The  plots  of  strain-rate  invariant  I  and  the  temperature  change  in  the  target  region,  shown  in  Fig. 
7,  confirm  that  significant  target  deforr  s  occur  in  the  vicinity  of  the  target/penetrator  interface. 


F.g.6.  Variation  of  the  hydrostatic  pressure  and  the  strain-rate  mesure  /  on  the  axial  line  for  different  values  of  o.  Explanations 
see  Fig.  5 


tjj  H 
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z  < 
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Fig.  7.  Temperature  and  strain-rate  distribution  in  the  target  region  for  a  =  3.0. 
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How  different  material  parameters  influence  the  deformations  of  the  target  is  demonstrated  by 
results  presented  in  Figs.  8  through  1 1 .  Figure  8  shows  that  strain-rate  hardening  increases  the 
normal  pressure  more  than  the  work-hardening  does.  The  inclusion  of  thermal  softening  affects 
little,  if  any,  the  normal  pressure  distribution  on  the  penetrator  nose.  Near  the  nose  tip  (Fig.  9)  the 
inclusion  of  work-hardening  and  strain-rate  hardening  decreases  the  value  of  the  strain-rate  invariant 
I  but  increases  the  temperature.  This  is  due  to  the  fact  that  these  hardening  effects  increase  the 
material's  flow  stress  and  cause  more  plastic  working  which  is  converted  into  heat.  The  thermal 
softening  has  a  noticeable  effect  on  the  temperature  distribution  at  the  penetrator  nose.  From  the 
plot  of  the  strain-hardening  parameter  t p  on  the  penetrator  nose  and  on  the  axial  line  in  Fig.  10, 


Fig.  8.  Effect  of  different  material  parameters  on  the  pressure 
distribution  at  the  penetrator  nose  (r,/rQ  =  2.0;  a  =1.0; 

-  full  model;  - no  thermal  softening; 

- strain  hardening  only; . rigid  perfectly  plastic) 


Fig.  9.  Effect  of  different  material  parameters  on  the 
temperature  and  strain-rate  distribution  at  the  penctra- 
tor  nose  (r./r0  =  2.0).  Explanations  see  Fig.  8 


Fig.  10.  Effect  of  different  material  parameters  on  the  strain-hardening  parameter  ip.  Explanations  see  Fig.  8 
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Fig.  11.  Effect  of  strain-hardening,  strain-rate  hardening  and  thermal  softening  on  relative  ;-velocity,  temperature  change 
and  the  hydrostatic  pressure  at  points  on  the  axial  line.  Explanations  see  Fig.  8 


one  may  conclude  that  the  thermal  softening  reduces  its  value  significantly,  especially  at  points  near 
the  nose  surface.  On  the  axial  line  as  well  as  on  the  nose  surface,  strain-rate  hardening  in  turn 
increases  strain-hardening.  We  should  add  that  a  steady-state  penetration  problem  is  being  studied 
and  thus  it  is  tacitly  assumed  that  the  increased  energy  required  for  deforming  the  target  is  available 
whenever  necessary.  Fig.  11  depicts  the  variation  of  the  hydrostatic  pressure,  temperature  and 
relative  r-velocity  on  the  axial  line.  The  hydrostatic  pressure,  the  temperature  and  the  absolute 
z-velocity  of  a  target  particle  on  the  axial  line  increase  with  the  inclusion  of  hardening  effects  but 
change  very  little  by  the  consideration  also  of  thermal  softening.  In  order  to  investigate  further  the 
effect  of  thermal  softening,  we  arbitrarily  doubled  the  value  of  the  thermal  softening  parameter  a. 
This  increased  the  value  of  the  strain-rate  invariant  I  but  changed  very  little  the  value  of  other 
quantities  on  the  penetrator  nose. 

On  the  axial  line  uniaxial  strain  conditions  prevail  approximately.  Thus  the  magnitude  of  the 
deviatoric  stress  s2z  at  a  point  should  equal  2/3  the  effective  flow  stress  a,  defined  as 

cr,  =  <r0(I  4-  b I)m ( 1  -  ad)(  1  +  \plVoT- 

Of  all  the  points  on  the  axial  line,  the  nose  tip  is  the  most  critical  one  since  the  strain-hardening 
parameter  \p  assumes  very  high  values  there.  In  Table  1  below  are  compared  the  values  of  (-s:z) 
and  2/3  o t  at  the  nose  tip.  Whereas  the  error  is  negligible  when  the  target  material  is  rigid/perfectly 
plastic,  it  is  rather  high  for  the  other  three  cases.  A  possible  reason  for  the  high  error  is  that  values 
of  ip  at  the  center  are  sensitive  to  the  value  of  the  artificial  viscosity  <5  eventhough  other  field  variables 
show  negligible  dependence  upon  the  precise  value  of  S  within  a  certain  range.  To  support  this 
reasoning,  we  list  in  Table  2  values  of  the  same  variables  but  with  the  effect  of  strain-hardening 
neglected.  Note  that  these  are  for  a  higher  value  of  the  speed  of  the  penetrator. 

Finally  we  remark  that  results  presented  here  are  valid  only  for  the  constitutive  model  used 
herein. 
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Table  1.  Values  of  sa  and  2/3  at  the  nose-tip  for  different  material  models  (a  =  I) 


Model 

~°a 

~P 

2/3  e. 

%  difference 

Perfectly  plastic 

6.34 

5.64 

0.7 

0.67 

4.5 

Strain-hardening 

7.05 

6.36 

0.69 

1.12 

40.2 

+  Rate  dependence 

11.07 

9.88 

1.19 

2.01 

40.7 

+  Thermal  softening 

11.10 

10.25 

0.85 

1.42 

40.1 

Table  2.  Values  of  sa  and  2/3  e. 

at  the  nose-tip  for  different  material  models  (a  = 

5) 

Model 

~°u 

-P 

2/3  o. 

%  difference 

Perfectly  plastic 

8.97 

8.35 

0.62 

0.67 

7.46 

Rate  dependence 

1136 

11.26 

1.10 

1.17 

5.98 

+  Thermal  softening 

12.28 

11.39 

0.89 

0.96 

7.29 

Fig.  12.  Variation  of  strain-hardening  parameter  y>  in  the  target  region,  (a  =  3.0,  rjr0  -  2.0) 


Conclusion 

The  computed  results  show  that  during  the  steady  state  portion  of  the  penetration  process,  the 
penetrator  nose  shape  has  a  significant  effect  on  the  deformations  of  the  target.  Whereas  the  strain- 
rates  are  higher  for  the  sharper  ellipsoidal  nose,  deformations  spread  to  a  larger  distance  away  from 
the  nose  surface  for  the  blunt  nose.  The  speed  of  the  penetrator  has  a  weak  effect  on  the  axial 
resisting  force  experienced  by  the  penetrator  even  though  the  hydrostatic  pressure  does  increase 
with  the  speed.  The  inclusion  of  thermal  softening  effects  increases  the  strain-rate  in  the  target 
material  but  does  not  alter  the  pressure  distribution  on  the  penetrator  nose. 
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EFFECT  OF  MULTIPLE  INITIAL  IMPERFECTIONS  ON 
THE  INITIATION  AND  GROWTH  OF  ADIABATIC 
SHEAR  BANDS  IN  NONPOLAR  AND  DIPOLAR 

MATERIALS 

Y.  W.  KWON  and  R.  C.  BATRA 

Department  of  Engineering  Mechanics.  University  of  Missoun-Rolla,  Rolla,  MO  65401-0249.  U  S  A. 

Abstract — Simple  shearing  deformations  of  a  viscoplastic  block  made  of  nonpolar  and  dipolar 
materials,  and  placed  in  a  hard  loading  device  are  studied.  Multiple  defects  in  the  block  are  modeled 
by  perturbing  the  uniform  temperature  within  the  block  when  the  material  just  starts  deforming 
plastically  to  that  given  by  a  cosine  function  which  assumes  relative  maximum  values  at  several  points 
in  the  block.  It  is  found  that  for  simple  materials,  the  deformation  localizes  at  points  where  the 
perturbed  temperature  has  relative  minima  when  the  average  applied  strain-rate  y0  is  500  s'1  and  at 
the  locations  of  the  relative  maxima  of  the  perturbed  temperature  when  the  applied  strain-rare  is  more 
than  1000  s'1.  This  transition  occurs  possibly  due  to  different  thermal  lengths  and  the  dme  scales 
associated  with  the  work-hardening  in  the  two  cases.  For  dipolar  materials  the  deformation  localizes 
near  the  boundaries  of  the  block  abutting  the  loading  device  when  y0  =  500  s_I  but  at  the  locations  of 
the  relative  maxima  of  the  perturbed  temperature  when  y0  =  50,000s*1.  For  both  simple  and  dipolar 
materials  the  initiation  of  the  localization  of  the  deformation  is  considerably  delayed  as  compared  to 
the  case  when  the  temperature  perturbation  has  only  one  bump  with  its  center  coinciding  with  the 
center  of  the  block. 


INTRODUCTION 

Recently  there  has  been  considerable  interest  in  the  study  of  the  localization  of  the  shearing 
deformation  in  bodies  being  deformed  at  very  high  strain-rates.  These  narrow  regions  of  severe 
deformation  are  usually  referred  to  as  adiabatic  shear  bands  because  of  the  little  time  available 
for  the  heat  generated  to  diffuse  away  to  colder  regions  of  the  body.  These  shear  bands  are 
believed  to  be  precursors  to  shear  fractures. 

Most  analytical  studies  (e.g.  Recht  [1],  Staker  [2],  Clifton  [3],  Bums  [4],  Bai  [5]  and  Shawtd 
et  al.  [6]),  aimed  at  delineating  factors  that  inhibit  or  enhance  the  initiation  and  growth  of 
adiabatic  shear  bands,  have  involved  analyzing  the  thermomechanical  deformations  of  a  block 
undergoing  simple  shearing  deformations.  Whereas  Recht,  Staker  and  Clifton  derived 
conditions  based  on  the  assumption  that  the  material  point  becomes  unstable  when  the  shear 
stress  at  that  point  attains  a  maximum  value.  Bums,  Bai  and  Shawki  et  al.  studied  the  growth 
of  small  perturbations  superimposed  on  a  finitely  deformed  body.  In  these  latter  analytical 
studies  the  governing  equations  were  linearized  around  the  finitely  deformed  state.  Clifton  et 
al.  [7J,  Wright  and  Batra  [8,  9j,  Wright  and  Walter  [10]  and  Batra  [11]  studied,  numerically,  the 
effect  of  introducing  a  perturbation  on  a  finitely  deformed  body  and  did  not  need  to  linearize 
the  governing  equations.  Subsequently  Batra  [12]  introduced  the  temperature  perturbation 
when  the  body  just  starts  deforming  plastically  and  investigated  the  effect  of  various  material 
parameters  on  the  value  of  the  average  strain  at  which  the  deformation  begins  to  localize. 
These  perturbations  are  supposed  to  simulate  the  flaws  or  inhomogeneities  present  in  the 
material.  Whereas  in  previous  works  one  or  two  flaws/defects  were  presumed  to  be  present, 
here  we  assume  that  the  flaws  are  periodically  distributed  and  represent  these  by  an  initial 
temperature  distribution  given  by  a  cosine  function  which  assumes  relative  maximum  values  at 
several  points  in  the  specimen.  It  is  found  that  the  presence  of  many  flaws/defects  delays 
considerably  the  initiation  of  the  localization  of  the  deformation. 

Wnght  and  Batra  [9]  and  Batra  [11]  have  shown  that  the  inclusion  of  the  strain  gradients  as 
an  independent  kinematic  variable  has  a  stiffening  effect  in  the  sense  that  it  delays  considerably 
the  onset  of  shear  bands.  Batra  [11]  also  studied  the  interaction  among  shear  bands  in  sir  le 
and  dipolar  materials.  His  numerical  calculations  revealed  that  two  bands  that  will  grow 
independently  in  simple  materials  will  ccalesce  in  dipolar  materials.  Herein  it  is  found  that  in 
dipolar  materials  the  deformation  localizes  near  the  boundaries  of  the  specimen  at  an  applied 
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strain-rate  y0  of  500  s_l  but  at  the  locations  of  the  relative  maxima  of  the  initial  temperature 
when  y0  =  50,000  s'1.  In  simple  materials  the  deformation  localized  at  points  between  the 
locations  of  the  relative  maxima  of  the  initial  temperature  at  y„  =  500  s~‘.  However,  at 
y„  =  50,000  s-1  the  points  where  the  deformation  localized  coincided  with  the  locations  of  the 
relative  maxima  of  the  initial  temperature. 


FORMULATION  OF  THE  PROBLEM 

Equations  governing  the  thermomechanical  deformations  of  a  block  occupying  the  region 
-1  s;y  ss  +1  and  made  of  a  viscoplastic  material  undergoing  simple  shearing  deformations  may 


be  written  [9, 11]  as 

t>  =  -  (s  -  /ay).y,  (1) 

d  =  kd  yy  +  A (s2  +  a2),  (2) 

i  =  p(uy-As),  (3) 

0=/fi(v.yy-J  O^,  (4) 

y-A<s--Mr)/(n-^J  .  (5) 


f 

/ 

A  =  max- 

0, 

-1 

/[b(s2  +  cr)1'2] 

LA'-  vv  /  J 

/ 

These  equations  are  written  in  terms  of  nondimensional  variables.  Equations  (1)  and  (2) 
represent,  respectively,  the  balance  of  linear  momentum  and  the  balance  of  internal  energy. 
The  balance  of  mass  simply  gives  that  the  mass  density  p  stays  constant  since  the  simple 
shearing  deformation  is  isochoric.  Here  v  is  the  x- velocity  of  a  material  particle,  Q  its 
temperature  change  from  that  in  the  reference  configuration,  s  is  the  shear  stress  in  the  x 
direction  acting  on  the  plane  y  =  constant,  a  is  the  dipolar  stress  on  this  plane,  k  the  thermal 
diffusivity,  parameters  y>0  and  n  describe  the  strain-hardening  of  the  material,  parameters  b 
and  m  characterize  the  strain-rate  sensitivity  of  the  material,  the  parameter  a  defines  the 
thermal  softening,  the  parameter  /  is  the  material  characteristic  length,  and  p  is  the  shear 
modulus.  The  parameter  may  be  thought  of  as  an  internal  variable  and  eqn  (5)  its  growth 
equation.  It  is  used  to  describe  the  work-hardening  of  the  material.  Note  that  the  numerator  on 
the  right-hand  side  of  eqn  (5)  equals  the  plastic  working.  A  comma  followed  by  y  signifies 
partial  differentii  .ion  with  respect  to  y  and  a  superimposed  dot  stands  for  the  material  time 
differentiation.  Implicit  in  eqn  (2)  are  the  assumptions  that  the  specific  internal  energy  equals 
the  specific  heat  multiplied  by  the  change  in  temperature,  Fourier's  law  of  heat  conduction 
holds  and  that  all  of  the  plastic  working  is  converted  into  heat.  The  material  characteristic 
length  /  equals  zero  for  nonpolar  materials  and  is  positive  for  dipolar  materials. 

In  eqns  ( 1)— (6)  it  has  also  been  assumed  that  the  strain  rate  y  =  v,r  and  its  gradient  d  =  v_n 
have  additive  decompositions  into  elastic  ye,  de  and  plastic  parts  yp,  dp.  That  is 

y  =  re  +  yp,  d  =  dt  +  dp.  (7) 

For  the  plastic  parts  yp  and  dp  we  have  made  the  following  constitutive  assumptions: 

Act 

yp  =  Ar,  dv  =  ~f-  (8) 

A  in  these  equations  equals  zero  whenever  the  deformations  are  elastic  and  is  positive  for 
plastic  deformations.  In  order  to  determine  whether  a  material  point  is  deforming  elastically  or 
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plastically,  we  presume  that  a  scalar  loading  or  yield  function  /  exists  such  that 

f(s,  a,  6,  yp,  dp)  =  k  (9) 

and 

df 

— -  (r,  a,  6,  As,  Aa)  <  0  (10) 

dA 

for  all  admiss.  ,le  values  of  s,  a  and  6.  k  in  eqn  (9)  describes  the  work  hardening  of  the 
material.  The  criterion  for  elastic  and  plastic  deformation  at  a  material  point  is 


f{s,  a,  6,  0,  0) =£  0,  elastic  (11) 

f(s,  a,  6,  0,  0)  >  0,  plastic.  (12) 


When  eqn  (12)  holds,  eqn  (9)  has  a  unique  solution  for  A  because  of  the  requirement  (10). 
Equation  (6)  gives  the  value  of  A  for  the  choices 


f(s,  a,  d,  yp,  dp)  = 


(s2  +  a2)* 

(1  -«*)[!  +  b{Y*  +  d*)lJ 


(13) 


and 


K 


n 


vv 


(14) 


The  details  of  deriving  eqns  (l)-(6)  may  be  found  in  Wright  and  Batra  [9]  and  Green  et  al. 
[13). 

The  governing  equations  (l)-(6)  are  to  be  supplemented  by  initial  conditions  and  boundary 
conditions.  For  the  former  we  take 


v(y,  0)-y,  d(y,  0)  =  0,  o(y,  0)  =  0,  y>{y,  0)  =  0, 

9(y,  0)  =  (1  +  cos  20ny)/20, 

5(y»  0)  =  [1  -  ad{y,  0)],  (15) 

and  for  the  latter 

u(±l,  t)  =  ±1,  0y(±l,r)  =  O,  a(±l,  t)  =  0.  (16) 

Thus  the  material  point  is  initially  presumed  to  lie  on  its  yield  surface  that  corresponds  to 
quasistatic  deformations,  no  work  hardening  and  its  initial  temperature.  The  initial  stress 
distribution  is  nonuniform  and  yp  =  0  initially  at  all  points  in  the  body. 

We  seek  solutions  of  eqns  ( 1 )— (6)  subjected  to  initial  and  boundary  conditions  (15)  and  (16) 
that  exhibit  the  following  symmetries  and  antisymmetries: 


v{-y,t)=*-v(y,t),  9(-y,t)  =  e(y,t), 

o(~y,t)  =  -o(y,t),  s{-y,t)  =  s(y,t),  y(-y,  t)  =  yr(y,  t). 


(17) 


However  we  do  not  assume  a  priori  that  the  solution  is  periodic.  Hence  we  can  reduce  the 
domain  of  study  to  [0, 1]  and  replace  boundary  conditions  (16)  by 

v(l,  t)  =  1,  u(0,  r)  =  0,  0y(l,/)  =  O,  6  y(0,  t)  =  0, 

(18) 

o(l,  r)  =0,  <7(0,  f)  =  0. 

The  nonlinear  coupled  governing  equations  (l)-(5)  under  the  side  conditions  (15)  and  (18)  are 
solved  numerically  by  the  finite  element  method.  We  note  that  there  are  no  existence  and 
uniqueness  theorems  available  for  such  a  system  of  equations.  Also  there  is  no  hope  of  solving 
these  equations  analytically,  therefore  we  solve  them  numerically. 
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COMPUTATION  AND  DISCUSSION  OF  RESULTS 
We  introduce  the  auxiliary  variables 

u=vy,  g  =  9.y,  p  =  o  y  (19) 


and  rewrite  eqns  (l)-(4)  as 


V  =-(S  -lp).y. 

(20) 

9  =  kgy  +  A(s:  +  a2), 

(21) 

s  =  p(u-As), 

(22) 

.  /  A  \ 

°  =  lP{u.y-j0). 

(23) 

Whereas  Batra  [11]  integrated  eqns  (19)-(23)  and  (5)  by  the  Crank-Nicholson-Galerkin  finite 
element  method  to  find  u,  g,  p,  v,  6,  s,  a  and  ip  at  a  node  point,  we  have  integrated  eqns 
( I )— (5)  by  the  same  method  and  computed  u,  d,  s,  a  and  ^  at  a  node  point.  For  a  test 
problem,  the  computed  results  by  the  two  methods  were  identical  but  the  core  storage  needed 
and  the  CPU  time  used  was  considerably  less  with  the  present  method.  Here  we  divided  the 
domain  [0.1]  into  400  uniform  subdomains  or  finite  elements  and  used  At  =  5  x  lO'6  in 
computing  the  results. 

The  values  of  the  nondimensional  parameters  a,  p,  n,  ip0  and  m  do  not  depend  upon  the 
applied  average  strain-rate  y0  but  those  of  p,  k  and  b  vary  as  Yo>  Yo'  and  y0  respectively.  The 
values  of  these  variables,  which  correspond  to  a  typical  hard  steel,  used  to  compute  numerical 
results  are 


a  =  0.4973,  p  =240.3,  n  =  0.09,  ip0  =  0-017,  m  =  0.025 

and 

p  =  3.928  x  10'5,  k  =  3.978  x  10-3,  b  =  5  x  106  when  y0  =  500s*1, 
p  =  3.928  x  10"\  k  =  3.978  x  10"5,  b  =  5  x  10*  when  y0  =  50,000s-'. 

The  presumed  value  of  the  thermal  softening  coefficient  a  is  approx.  7  times  the  value  for  a 
typical  steel.  This  is  done  so  that  the  peak  in  the  shear  stress-shear  strain  curve  occurs  at  a 
lower  value  of  the  average  strain  and  therefore  the  computational  time  for  the  problem  is 
significantly  reduced.  This  choice  of  the  value  of  a  should  not  affect  the  qualitative  nature  of 
results  presented  here. 

For  homogeneous  deformations  of  the  block  Fig.  I  depicts  the  shear  stress-shear  strain  curve 
for  an  applied  strain-rate  of  500  s"1.  The  peak  in  the  curve  occurs  at  an  average  strain  of  0.093. 
At  the  higher  strain-rate  of  50,000  s"1  the  shear  stress  peak*-  *  at  an  average  strain  of  0.085. 
The  temperature  perturbation  was  introduced  at  an  a»-  -  «ain  corresponding  to  point  I  in 
Fig.  1. 

Figure  2  shows  how  the  temperature,  plastic  strain-rate  and  the  shear  stress  s  evolve  in  the 
steel  block  when  it  is  modeled  as  a  nonpolar  material,  the  initial  temperature  distribution  is 
nonuniform,  and  the  block  is  deformed  at  an  average  strain-rate  of  500  s"1.  Figure  3  depicts  the 
corresponding  results  for  an  applied  strain-rate  of  50,000s"1.  Note  that  the  initial  temperature 
has  the  relative  maximum  value  of  0.10  at  y  =0.0,  0.1,  0.2,  0.3,  0.4,  0.5,  0.6,  0.7,  0.8,  0.9  and 
1.0.  In  each  case  the  shear  stress  s  becomes  essentially  uniform  after  a  short  while  and  stays 
uniform  throughout  the  thickness  of  the  block.  It  first  rises  because  of  the  work  hardening  and 
strain-rate  hardening  effects  and  begins  to  decrease  when  these  hardening  effects  are  overcome 
by  the  thermal  softening.  At  the  lower  strain-rate  of  500s"'  there  is  no  noticeable  drop  in  the 
shear  stress  even  when  the  deformation  has  started  to  localize  but  at  the  higher  strain-rate  of 
50,000  s1  the  shear  stress  s  decreases  considerably  after  the  deformation  begins  to  localize. 
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Fig.  1.  Shear  stress-strain  curve  for  simple  shearing  deformations  of  a  steel  block  at  y0  =  500  s'1 


This  rather  noticeable  drop  in  the  shear  stress  with  increasing  average  strain  after  the 
deformation  begins  to  localize  has  been  pointed  out  by  Wright  and  Walter  [10].  Note  that  the 
temperature  rise  in  the  two  cases  at  places  where  the  deformation  localizes  is  nearly  equal  to 
each  other  at  the  instant  of  the  onset  of  the  localization  process.  However,  the  temperature  is 
more  evenly  distributed  at  y0  =  500  s-1  as  compared  to  that  at  y0  =  50,000  s'1  in  the  sense  that 
the  difference  between  the  maximum  and  the  minimum  temperature  stays  small  for  y„  =  500  s'1 
even  after  the  deformation  localizes.  The  rate  of  growth  of  the  plastic  strain-rate  is  more 
gradual  at  y0  =  500  s-1  as  compared  to  that  at  y0  =  50,000  s'1. 

The  striking  differences  between  the  two  cases  are  outlined  below. 

(a)  First  is  the  average  strain  at  the  instant  the  deformation  begins  to  localize.  For 
yo  =  500s-1  the  average  strain  of  0.624  at  which  the  deformation  has  essentially  localized  is 
nearly  7  times  the  strain  at  which  the  peak  stress  occurs  during  homogeneous  deformations  of 
the  block.  However,  for  y0  =  50,000  s'1  the  localization  of  the  deformation  occurs  at  an 
average  strain  of  0.09  which  almost  equals  the  strain  of  0.0S5  at  which  the  shear  stress  attains 
the  maximum  value  during  homogeneous  deformations  of  the  block.  For  comparison  purposes 
we  note  [14]  that  when  a  temperature  perturbation  with  a  single  bump  centered  at  y  =0  is 
introduced,  the  shear  bands  form  at  an  average  strain  of  0.0814  and  0.206  for  y0  =  500  s'1  and 
50,000  s'1,  respectively.  On  the  hypothesis  that  a  temperature  perturbation  simulates  material 
inhomogeneities  or  flaws  in  the  body,  the  present  results  indicate  the  initiation  of  shear  bands  is 
delayed  for  yn  =  500s_l  when  there  are  more  defects  uniformly  distributed  in  the  specimen, 
and  for  y0  =  50,000  s'1  when  there  is  only  one  defect  present  in  the  specimen. 

(b)  Secondly,  the  places  where  the  peak  temperature  rise  and  the  peak  plastic  strain-rate 
occur  are  quite  different  in  the  two  cases.  Whereas  at  yo  =  500s'1  the  peak  values  of  the 
temperature  rise  and  the  plastic  strain-rate  occur  at  places  where  the  initial  temperature 
perturbation  had  relative  minimum  values,  at  y0  =  50,000  s'1  the  peaks  of  the  temperature  rise 
and  the  plastic  strain-rate  coincide  with  the  locations  of  the  maxima  of  the  temperature 
perturbation.  That  the  centers  of  shear  bands  for  y0  =  50,000  s'1  coincide  with  the  locations  of 
the  maxima  of  the  initial  temperature  perturbation  was  also  observed  in  the  case  the  initial 
temperature  perturbation  assumed  peak  values  at  six  equidistant  points  rather  than  eleven.  In 
this  case,  the  deformation  localized  at  an  average  strain  of  0.il22. 

In  order  to  understand  the  reasons  for  the  above  differences  we  note  that,  for  nonpolar 
materials  (i.e.  /  =  0),  imbedded  in  the  governing  equations  (l)-(6)  are  three  length  scales 
[9, 12]  namely,  the  geometric  length,  a  thermal  length  and  a  viscous  length,  and  seven 
non-dimensional  parameters,  namely,  the  mass  density  p,  the  elastic  modulus  p,  the  rate  of 
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Fig.  2.  Evolution  of  the  temperature-,  plastic  strain-rate  and  the  shearing  stress  in  nonpolar  materials 

at  Yo  =  500  s'1. 


ther  nal  softening  a,  two  for  work  hardening  n  and  V/„.  and  two  for  rate  hardening  b  and  m. 
Whereas  the  thermal  length  varies  as  Yo'n,  the  viscous  length  does  not  depend  upon  y0  but  is  a 
function  of  the  material  parameters.  The  non-dimensional  thermal  length  ( =  thermal 
length/height  of  the  specimen)  is  reduced  from  0.6307  to  0.006307  when  y0  is  increased 
from  500  to  50,000 s“l.  That  the  decrease  in  the  thermal  length  is  among  the  factors  responsible 
for  the  aforementioned  differences  was  confirmed  by  running  the  yo=500s_l  case  with  the 
thermal  conductivity  arbitrarily  set  equal  to  zero.  In  this  case  the  thermal  length  is  zero  and  the 
peak  strain-rates  occurred  at  the  locations  of  the  maxima  of  the  initial  perturbation  as 
expected.  To  seek  an  answer  to  the  question  “For  what  value  of  y„  does  the  above  referred 
transition  occur?”,  we  conducted  numerical  experiments  for  different  values  of  y0  between  500 
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Fig.  3.  Evolution  of  the  temperature,  plastic  strain-rate  and  the  shear  stress  in  nonpolar  materials  at 

y0  =  50,000  s-'. 

and  50.000  s-1  and  found  y0=  1000  s-1  to  be  the  critical  strain-rate  above  which  the 
deformation  localized  at  points  of  maxima  of  the  initial  temperature  perturbation  and  below 
which  it  localized  at  other  points. 

It  appears  that  it  is  the  higher  value  of  the  thermal  length  at  yo  =  500s-1  that  delays  the 
initiation  of  the  localization  process  too.  A  close  examination  of  the  computed  results  indicated 
that  initially  the  material  points  at  the  locations  of  the  temperature  maxima  experienced  peak 
plastic  strain-rates  and  higher  values  of  the  work  hardening  parameter  rfi  as  compared  to  their 
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neighbors.  Since  the  shear  stress  became  essentially  uniform  shortly  after  the  temperature 
perturbation  was  introduced,  the  increased  value  of  ip  resulted  in  reduced  values  of  A,  yp  and 
the  heat  generated  due  to  plastic  working  at  these  points.  This  coupled  with  the  heat 
conduction  tended  to  make  the  temperature  virtually  uniform  throughout  the  specimen.  Of 
course  the  temperature  and  work-hardening  parameter  did  not  become  exactly  uniform 
because  if  they  did,  the  specimen  will  subsequently  undergo  homogeneous  deformations  only. 
To  lend  credence  to  this  hypothesis  we  increased  artificially  the  thermal  conductivity  to  100 
times  its  previous  value  but  kept  yo  =  50Os-1.  In  this  case  no  shear  band  formed  and  the 
temperature  held  within  the  specimen  became  uniform  and  stayed  so  during  the  ensuing 
deformations.  To  understand  the  role  played  by  the  work-hardening,  we  conducted  a  numerical 
experiment  with  n=0  and  yo  =  50Os-1.  This  resulted  in  neglecting  the  effect  of  the 
work-hardening  on  the  deformations  of  the  block.  In  this  case  the  deformation  localized  at 
points  where  the  initial  temperature  perturbation  assumed  peak  values. 

Another  reason  for  the  localization  not  to  occur  at  the  locations  of  the  peaks  in  the  initial 
temperature  could  be  improper  choices  of  the  mesh  and  At.  We  reduced  At  to  one-half  of  its 
value,  but  used  the  same  mesh  and  recomputed  results  for  y0  =  500  s-1.  For  smaller  time  step 
the  deformation  localized  at  the  same  points  as  it  did  earlier  for  the  large  time  step. 

Since  the  distance  between  the  locations  of  the  points  where  the  plastic  strain-rate  eventually 
peaked  and  points  of  maxima  of  the  initial  temperature  perturbation  equals  approximately  the 
length  of  20  elements,  it  is  unlikely  that  the  finite-element  mesh  caused  the  shift. 

When  y0  is  increased  from  500  to  50,000  s-1,  the .  non-dimensional  mass  density  increases 
from  3.93  x  10-5  to  3.93  x  10-1.  In  [14]  the  initial  temperature  perturbation  had  only  one  bump 
centered  at  y  =  0.  It  was  found  that  the  shear  stress  distribution  within  the  specimen  was 
uniform  for  y0  =  500 s-1  but  non-uniform  with  the  lowest  value  at  y  =  0  for  y0  =  50,000 s"‘.  No 
such  non-uniformity  in  the  shear  stress  distribution  has  been  observed  for  the  periodic  initial 
temperature  perturbation  introduced  herein. 

In  Figs  4  and  5  are  plotted  the  evolution  of  the  temperature,  plastic  strain-rate  and  the  shear 
stress  s  when  the  material  of  the  block  is  modeled  as  a  dipolar  material  with  /  =  0.01  and  the 
block  is  deformed  at  an  average  strain-rate  y0  of  500  and  50,000  s-1  respectively.  Whereas  for 
yo=500s-1  only  one  shear  band  with  center  at  y  =  1.0  forms,  at  y0  =  50,000s-1  the 
deformation  localizes  around  the  location  of  the  peaks  of  the  initial  temperature.  In  each  case 
the  initiation  of  the  localization  of  the  deformation  is  delayed  as  compared  to  that  for  nonpolar 
materials  [14];  the  delay  is  more  pronounced  at  the  higher  applied  strain-rate.  Whereas  for 
nonpolar  materials  the  average  strain  at  which  the  deformation  localizes  for  yo  =  500s-1  is 
nearly  seven  times  that  for  y„  =  50,000  s-1,  for  dipolar  materials  the  average  strains  at  which 
the  deformation  localizes  in  the  two  cases  are  nearly  equal  to  each  other.  Also  for  nonpolar 
materials  the  shear  stress  s  becomes  essentially  uniform  throughout  the  specimen  soon  after  the 
initial  perturbation  is  introduced  and  stays  uniform  within  the  block,  such  is  not  the  case  for 
dipolar  materials.  For  dipolar  materials  the  amplitude  of  the  oscillations  in  the  shear  stress 
distribution  is  extremely  small  for  yo  =  500s-1  but  at  y0  =  50,000s-1  it  is  noticeable  and  seems 
to  grow  as  the  deformation  localizes.  Note  that  for  dipolar  materials  it  is  (j  -  lo  y)  and  not  j 
that  acts  as  a  flux  for  the  linear  momentum.  The  average  shear  stress  s  exhibits  the  same 
behavior  as  that  for  nonpolar  materials  in  the  sense  that  it  first  increases  and  subsequently 
begins  to  drop  when  thermal  softening  effects  exceed  the  combined  effects  of  strain  and 
strain-rate  hardening. 

At  y„  =  500s-1  the  initial  temperature  distribution  tends  to  become  uniform  throughout  the 
block  and  stays  virtually  uniform  until  the  deformation  begins  to  localize  near  y  =  1.0.  At  that 
instant  the  temperature  near  y  =  1.0  increases  slightly  more  than  that  at  other  places.  Of  course 
the  temperature  at  every  point  in  the  body  keeps  on  increasing  because  of  the  plastic  working. 
Why  the  shear  band  forms  near  y  =  1.0  and  not  near  y  =  0.0  is  not  clear.  At  the  higher  applied 
strain-rate  y0  =  50,000  s-1  neither  the  temperature  nor  the  plastic  strain-rate  ever  becomes 
uniform  within  the  block.  It  seems  that  the  amplitude  of  the  temperature  oscillations  keeps  on 
increasing  with  the  increase  in  the  average  strain  in  the  specimen. 

For  dipolar  materials,  in  addition  to  the  thermal  and  viscous  lengths,  there  are  three  length 
scales  [13]  characteristic  of  the  material.  Herein  as  well  as  in  two  previous  studies  [9,  11]  the 
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Fig.  4.  Evolution  of  the  temperature,  plastic  strain-rate  and  the  shear  stress  in  dipolar  materials 

(/  =  0.01)  at  y0  =  500  s’1. 

three  length  scales  have  been  set  equal  to  each  other.  The  problem  when  the  three  length 
scales  are  different  is  under  investigation  and  will  be  reported  on  in  a  future  paper.  The 
previous  work  with  yo  =  500s-'  and  with  a  single  temperature  perturbation  centered  at  y  =  0 
indicated  that  the  inclusion  of  dipolar  effects  delays  considerably  the  localization  of  the 
deformation.  The  delay  observed  here  by  comparing  results  plotted  in  Figs  3  and  5  is 
significantly  more  as  compared  to  that  found  previously  and  also  that  computed  presently  for 
y0  -  50,000s"1.  The  main  reason  for  the  difference  between  the  previously  computed  [9, 11] 
and  current  results  is  due  to  the  number  of  flaws/imperfections  considered  in  the  two  cases.  An 
explanation  for  the  difference  between  the  results  computed  at  yo  =  500  and  50.000s-'  lies  in 
the  different  values  of  the  thermal  lengths  in  the  two  cases.  At  y0  =  50Qs-1  the  temperature 
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Fig.  5.  Evolution  of  the  temperature,  plastic  strain-rate  and  the  shear  stress  in  dipolar  materials 

(/  =  (). 01)  at  y„  =  50,000s-'. 

initially  tends  to  become  uniform  because  of  heat  conduction  and  once  the  temperature 
fluctuations  die  out  the  block  deforms  homogeneously.  What  eventually  causes  the  temperature 
to  rise  at  y  =  1.0  and  not  at  y  =  0.0  is  unclear. 

We  note  that  the  CPU  time  required  to  compute  results  for  the  dipolar  materials  is  nearly 
four  times  that  needed  -for  non-polar  materials.  Thus  not  many  numerical  experiments  could  be 
conducted  for  dipolar  materials. 


CONCLUSIONS 

For  nonpolar  materials  the  initiation  of  the  shear  bands  is  significantly  delayed  at  y0  =  500  s~ 1 
when  a  periodic  temperature  perturbation  with  eleven  peaks  is  introduced  initially  as  compared 
to  the  case  when  the  initial  temperature  perturbation  has  a  single  bump  with  peak  at  the  center 
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of  the  block.  However,  for  y0  =  50,000  s'1  the  reverse  happens.  In  the  former  case  the 
deformation  localizes  at  points  where  the  initial  temperature  has  relative  minima  values,  in  the 
latter  case  the  centers  of  shear  bands  coincide  with  the  places  where  the  initial  temperature  has 
relative  maximum  values.  These  transitions  are  found  to  occur  at  y0  =  1000  s'1  and  seem  to  be 
caused  by  the  lower  value  of  the  thermal  length  at  the  higher  value  of  y0  and  the  different  time 
scale  associated  with  the  work-hardening  in  the  two  cases. 

For  dipolar  materials  the  average  strains  in  the  specimen  at  which  the  deformation  localizes 
at  y0  =  500  and  50,000  s"1  are  nearly  equal  to  each  other.  These  values  of  average  strains  are 
higher  than  the  corresponding  values  for  nonpolar  materials.  Only  one  band  forms  at 
yo  =  500s'*  but  the  number  of  shear  bands  equal  the  number  of  peaks  in  the  initial 
temperature  at  y0  =  50,000  s-1. 
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THE  INITIATION  AND  GROWTH  OF,  AND  THE 
INTERACTION  AMONG,  ADIABATIC  SHEAR  BANDS 
IN  SIMPLE  AND  DIPOLAR  MATERIALS 
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Abstract  The  problem  of  simple  shearing  of  a  block  of  simple  (nonpolar)  and  dipolar  ther- 
moviscoplastic  materials  is  studied  with  the  objectives  of  exploring  the  initiation  and  growth 
of,  and  the  interaction  among,  adiabatic  shear  bands.  A  shear  band  is  assumed  to  have  iormed 
if  the  addition  of  a  perturbation  to  the  homogeneous  fields  just  before  the  peak  stress  is  reached 
results  in  a  localization  of  the  shear  strain.  The  effect  of  adding  perturbations  of  different  sizes 
and  of  the  same  size  but  at  different  locations  in  the  slab  is  investigated.  It  is  shown  that  in 
simple  materials,  two  shear  bands  coalesce  if  the  distance  between  them  is  small  but  grow  inde¬ 
pendently,  although  at  a  slower  rate,  if  the  distance  between  them  is  large.  However,  for  dipolar 
materials,  the  two  bands  coalesce  even  when  the  distance  between  them  at  the  time  ol  their  initi¬ 
ation  is  20  times  the  material  characteristic  length. 


I.  INTRODUCTION 

Adiabatic  shear  is  the  name  given  to  a  localization  phenomenon  that  occurs  during  high- 
rate  plastic  deformation,  such  as  machining,  explosive  forming,  shock  impact  loading, 
ballistic  penetration,  fragmentation,  ore  crushing,  impact  tooling  failure,  and  metal 
shaping  and  forming  processes.  The  localization  of  shear  strain  has  been  observed 
mostly  in  steels,  but  also  in  nonferrous  metals  and  polymers.  Practical  interest  in  the 
phenomenon  derives  from  the  fact  that  progressive  shearing  on  an  intense  shear  band 
provides  an  undesirable  mode  of  material  resistance  to  imposed  deformation,  and  the 
bands  are  often  precursors  to  shear  fractures. 

Shear  band  formation  is  generally  enhanced  at  high  strain  rates  because  the  lack  of 
time  for  heat  diffusion  'P  nonuniform  straining  to  cause  nonunilorm  heating.  Non- 
uniform  temperatures  enhance  plastic  flow  in  the  hotter  regions  and  reduce  plastic  flow 
in  the  colder  regions.  Furthermore,  heat  generation  is  greatest  in  the  regions  of  high¬ 
est  strain  rate.  Thus,  the  strain  rates  in  the  hot,  high-strain-rate  regions  tend  to  become 
larger,  while  those  in  the  cold,  low-strain-rate  regions  tend  to  become  more  nonuniform 
and  may  localize  into  a  narrow  region  referred  to  as  a  shear  band.  Whether  or  not  this 
thermoplastic  instability  mechanism  leads  to  shear  bands  seems  to  depend  upon  strain 
hardening,  strain-rate  hardening,  thermal  softening,  thermal  diffusivity,  and  the  strength 
of  the  initial  inhomogeneity. 

Zener  and  Holloman  [1944]  recognized  the  destabilizing  effect  of  thermal  soften¬ 
ing  in  reducing  the  slope  of  the  stress-strain  curve  in  nearly  adiabatic  deformations.  The 
dynamic  torsional  experiments  of  Culver  (1973)  on  mild  steel,  titanium,  and  6061 -T6 
aluminum  indicate  that  the  localization  began  near  the  peak  in  the  stress-strain  curve 
for  each  material  tested.  This  observation  seems  to  be  borne  out  by  the  experimental 
work  of  Costin  et  al.  (1980).  They  found  that  in  dynamic  torsion  tests  on  short  speci¬ 
mens  of  10' S  CRS  and  1020  HRS,  the  shear  bands  appeared  in  the  1018  CRS  specimens 
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and  not  in  1020  HRS  specimens.  Extensive  experiments,  conducted  at  different  temper¬ 
atures  and  strain  rates,  revealed  that  the  shear  stress-shear  strain  curve  for  HRS  was 
increasing  monotonically,  whereas  the  curve  for  CRS  had  a  peak  in  it. 

Johnson  ei  al.  [1983]  and  Lindholm  and  Johnson  [1983]  reported  dynamic  torsion 
test  data  on  six  ductile  metals  and  six  additional  materials  of  much  less  ductility.  They 
proposed  a  constitutive  relation  that  accounts  for  strain  and  strain-rate  hardening  and 
thermal  softening.  Their  analysis  indicates  that  high  ductility  tends  to  produce  a  rela¬ 
tively  wide  shear  band  when  compared  to  other  materials.  In  their  experimental  setup, 
the  frictional  force  at  the  grips  provided  a  constraining  axial  force.  Thus,  an  axial  load 
or  stress  component  developed  with  increasing  torsional  deformation.  However,  this 
axial  load  component  was  not  measured. 

Staker  [1981]  made  use  of  an  instability  analysis  to  model  the  appearance  of  adia¬ 
batic  shear  bands  in  the  controlled  explosive  expansion  of  steel  cylinders.  He  argued 
that  because  of  the  short  times  involved  in  explosive  loading  the  deformation  could  be 
considered  adiabatic  and  there  was  no  need  to  consider  the  thermophysical  properties 
of  materials.  Rhcht  [1964]  included  heat  conduction  and  thermal  softening  in  the 
interpretation  of  shear  bands  formed  during  high-speed  machining.  In  his  investigation 
of  instability  in  the  shear  zone  ahead  of  a  cutting  tool  during  orthogonal  machining, 
he  used  a  thermal  model  incorporating  uniform,  constant-rate  heat  generation  at  a  plane 
in  an  infinite  medium.  He  showed  that  the  critical  strain  rate  for  catastrophic  shear  in 
mild  steel  is  1400  times  greater  than  that  for  titanium.  Based  on  these  data  it  can  be 
shown  that  the  difference  in  the  thermophvsical  properties  of  St  and  Ti  contributes  a 
factor  of  6,  the  ratio  of  yield  stresses  a  factor  of  4,  and  the  difference  in  mechanical 
behavior  a  factor  of  58.  This  last  factor  is  directly  proportional  to  the  material’s  capacity 
to  strain-harden  and  inversely  proportional  to  its  tendency  to  thermally  soften. 

In  a  departure  from  the  notion  of  a  criterion  based  on  a  stress  maximum,  Clifton 
[1980]  and  Bai  [1981]  examined  the  growth  of  infinitesimal  periodic  nonunilormities 
in  an  otherwise  uniform  simple  shearing  deformation.  They  included  strain  hardening, 
thermal  softening,  strain  rate  sensitivity,  and  heat  conduction.  Bai  also  included  the 
effect  of  inertia  forces.  Both  linearized  the  field  equations  about  the  unperturbed  time- 
dependent  homogeneous  deformation  state  and  found  that  the  magnitude  of  the  im¬ 
perfections  may  grow  or  decay  in  time,  depending  on  the  material  parameters,  h.  ; 
average  rate  of  strain,  and  the  fixed  spatial  wavelength  of  the  initial  imperfection.  Birns 
[1983]  used  a  dual  asymptotic  expansion  to  include  the  time  dependence  of  the  homo¬ 
geneous  solution  in  the  analysis  of  the  growth  of  an  initially  small  perturbation.  His 
work  suggests  that  initiation  of  an  unstable  shear  band,  followed  by  exponential  growth, 
occurs  after  a  critical  shear  strain  corresponding  to  the  peak  stress  in  the  homogeneous 
deformation  for  the  same  overall  strain  rate  is  reached.  On  the  other  hand.  Shawki  ei 
al.  (I983|,  by  using  both  numerical  and  perturbation  techniques,  concluded  that  expo¬ 
nential  growth  is  not  a  sufficient  condition  for  judging  whether  or  not  a  shear  band 
forms,  as  the  corresponding  homogeneous  deformation  may  also  grow  extremely  rap¬ 
idly  once  the  peak  stress  has  been  reached  and  growth  is  not  restricted  to  a  narrow  band. 

Erlich  ei  al.  [1980]  noted  that  according  to  simple  wave  theory  applied  to  one¬ 
dimensional  plastic  wave  propagation,  the  strain  level  at  which  the  shear  tangent  modu¬ 
lus  becomes  zero  propagates  at  zero  speed.  They  postulated  a  criterion  of  adiabatic  shear 
band  formation  based  on  this  "wave  trapping"  idea.  This  idea  was  applied  in  Olson 
ei  al.  [1981]  in  a  numerical  finite  element  simulation  of  plastic  shear  wave  propagation 
under  adiabatic  conditions.  The  numerical  solution  did  indeed  exhibit  a  concent  i  at  ion 
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of  shear  strain  in  a  layer  of  finite  elements  that  was  similar  in  nature  to  observed  adi¬ 
abatic  shear  bands. 

In  studying  the  growth  of  shear  bands  in  the  center  of  a  finite  slab  after  initiation 
at  a  small  imperfection,  Merzer  [1983]  concluded  that  the  final  width  of  the  band 
depends  on  the  thermal  diffusivity  and  the  overall  strain  rate.  Wu  and  Freund  [1984], 
in  studying  the  formation  of  shear  bands  at  a  moving  boundary,  concluded  that  ther¬ 
mal  diffusivity  has  little  influence  on  the  final  shape  of  the  band.  The  detailed  geome¬ 
try  and  constitutive  equations  considered  in  these  two  papers  are  quite  different,  so 
perhaps  it  is  not  surprising  to  find  apparently  contradictory  results.  In  fact,  in  both 
papers  there  are  two  natural  length  scales:  one  arising  from  the  rate  effect  in  the  con¬ 
stitutive  equation  and  one  arising  from  heat  conductivity.  In  the  latter  paper,  these  scales 
have  been  arbitrarily  set  equal  to  each  other,  whereas  in  the  former  paper  the  relative 
effect  of  heat  conductivity  has  been  examined  parametrically  for  at  least  one  specific 
type  of  constitutive  equation. 

Recently,  Wright  and  Batra  [1985a,  1985b,  and  1986]  described  the  results  of  com¬ 
putations  that  simulate  the  formation  of  a  single  shear  band  from  a  local  temperature 
inhomogeneity  in  simple  and  dipolar  materials.  A  general  theory  of  thermoviscoplastic¬ 
ity,  obtained  by  modifying  the  dipolar  theory  of  Green  ei  al.  (1968)  to  include  rate 
effects,  was  used.  Wright  and  Batra’s  calculations  for  simple  materials,  as  well  as  rhe 
experimental  observations  of  Moss  [1981],  indicated  that  peak  strain  gradients  reached 
0.2  per  (rm  or  higher.  Therefore,  they  considered  worthwhile  to  investigate  the  dipo¬ 
lar  effects.  Their  computations  show  increasing  inhomogeneity  in  the  strain-,  temper¬ 
ature-,  and  strain-rate  fields,  with  the  central  amplitudes  growing  at  an  accelerating  rate. 
The  inclusion  of  dipolar  effects  has  a  stiffening  effect  in  the  sense  that  the  rate  of 
growth  of  central  amplitudes  of  the  strain,  temperature,  and  strain-rate  fields  is  lower 
as  compared  to  that  for  simple  materials.  For  dipolar  materials  there  are  at  least  three 
length  scales:  one  is  from  viscous  stress  effects,  the  second  is  from  thermal  conductivity, 
and  the  third  is  the  material  characteristic  length.  Wright  and  Batra  [1986]  studied 
the  case  when  all  three  length  scales  are  equal  to  one  another. 

This  paper  describes  the  results  of  some  numerical  experiments  conducted  with  the 
objective  of  analyzing  the  interaction  between  two  shear  bands.  It  also  examines  the 
effect  of  the  amplitude  and  distribution  of  the  initial  inhomogeneity  on  the  initiation 
and  growth  of  a  shear  band.  It  is  shown  that  a  narrower  inhomogeneity  results  in  a 
rapid  growth  of  the  band  as  compared  to  a  wider  one  having  the  same  central  ampli¬ 
tude.  A  stronger  inhomogeneity  results  in  the  formation  of  a  shear  band  even  before 
the  peak  in  the  homogeneous  stress-strain  curve  is  reached.  The  two  shear  bands  that 
would  grow  independently  in  simple  materials  seem  to  coalesce  in  dipolar  materials  even 
when  the  material  characteristic  length  is  of  the  distance  between  them  at  the  time 
of  their  initiation. 

We  note  that  there  is  no  experimental  evidence  available  on  the  interaction  between 
two  or  more  shear  bands  in  simple  or  dipolar  materials. 


II.  FORMULATION  OK  THE  SIMPLE  SHEARING  PROBLEM 

We  study  the  simple  shearing  deformations  of  a  dipolar  viscoplastic  material  and 
assume  that  all  of  the  variables  have  been  nondimensionalized.  Thus,  the  body  occu¬ 
pies  the  infinite  slab  bounded  by  the  planes  v  =  ±1.  Referring  the  reader  to  Wright  and 
Batra  [1986]  for  details,  we  note  that  the  governing  equations  are 
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I )  =  -  (s  -  lo,v),v,  (i) 

P 

d  =  kd,yv  +  A(s~  +  a~),  (2) 

s  =  n{v,y-As),  (3) 

a  =  /ft|v,lv-  y  <rj,  (4) 


with  boundary  conditions 

t>(  ±  I,/)  =  ±  I,  <?,,.(  ±  I,/)  =  0,  o(±l,r)=0  (7) 

and  a  suitable  set  of  initial  conditions.  Equations  (1)  and  (2)  express  the  balance  of  linear 
momentum  and  internal  energy,  respectively.  Herein  v  is  the  velocity  of  a  material  par¬ 
ticle,  /  is  a  material  characteristic  length,  $  is  the  temperature  change  of  the  material 
particle  from  that  in  the  reference  configuration,  and  s  and  a  may  be  interpreted  as  the 
shear  stress  and  the  dipolar  shear  stress,  respectively.  A  superimposed  dot  indicates 
material  time  differentiation,  and  a  comma  followed  by  y  signifies  partial  differentia¬ 
tion  with  respect  toy.  The  nondimensional  variables  are  related  to  their  dimensional 
(barred)  counterparts  as  follows: 

y  =  Y/H ,  t  =  iytu  y  =  c.y.  d  -  c.yy.  y  =  y,.  +  y,,,  d  -  d,  +■ 
v  =  v  Hyn,  s  =  $/k„,  <i  —  (T//kii,  tt  =  tfpc  s„,  (8) 

7=7,  d-(lH ,  7y,  =  V  7,>.  (i,,  =  (i..H  7,,,  .l=.U„7„. 


Besides  w,  «,  and  there  are  six  other  nondimensional  parameters,  which  are  related 
to  their  dimensional  (barred)  counterparts  as  follows: 


u  =  dK„.  pc,  b  -  hy„,  k=k  pCyt,H2.  I  =  i  H. 
M  =  M'  An.  P  =  pH-yr>  a„. 


(9) 


In  cqns  (8)  and  (9),  7,,  =  r(  H,  i)/H  is  the  average  applied  strain  rate  between  the 
boundaries  Y  -  ±H,  and  *(>  is  the  yield  stress  in  a  reference  quasistatic  test. 

The  constitutive  relations  (3)-(6)  give  one  possible  model  of  viscoplastic  materials. 
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Equation  (6)  implies  that  the  plastic  pans,  .15  and  Ao/I,  of  the  strain  rate  and  the  dipolar 
strain  rate  vanish  when 


(s2  +  a2)'  2<  (l  +  (1  -ad).  (10) 

The  material  parameters  i/-,,  and  n  describe  the  strain  hardening  of  the  material,  a  the 
thermal  softening,  and  b  and  m  the  strain-rate  sensitivity  of  the  material. 

We  presume  that  the  initial  values  of  6,  s,  and  \J,  are  symmetric  and  the  initial  val¬ 
ues  of  v  and  a  are  antisymmetric  in  y  and  seek  solutions  of  eqns  (I)  through  (6)  with 
the  same  symmetry.  Thus,  the  problem  is  to  be  studied  over  the  spatial  domain  |0,1] 
and  the  boundary  conditions  become 


r(0,n 

=  0, 

e.AO.n 

=  0, 

© 

II 

o 

(ii) 

Vll.t) 

=  1, 

8. A  1,/) 

=  0, 

0(1. n  =o. 

(12) 

For  the  initial  conditions  we  take 

r(y,0)  =  \\  o(y,0)=0.  i/dy.O)  =  i£,„ 

8(y.O)  =0„  +  tf(r>.  (13) 

5(y,0)  =  5d  =  ( 1  +  ~  |  ( 1  +  b.l5„)"'(  1  -a8{\\0)). 

The  values  of  8„,  s„,  and  are  such  that  during  homogeneous  deformations  of  the 
block  the  shear  stress  5„  and  the  strain  corresponding  to  i£„  lie  on  the  shear  stress-shear 
strain  curve  for  the  material.  .1  in  eqn  (!3)<;  is  given  by  eqn  (6),  with  8  =  8tt,  s  =  s„, 
=  ij(l,  and  a  =  0.  The  function  8  describes  the  aberration  in  the. initial  temperature 
distribution  and  will  result  in  nonhomogeneous  deformations  of  the  body. 


III.  NUMERIC  AL  INTEGRATION  OE  GOVERNING  EQUATIONS 
With  the  auxiliary  variables 

u  =  v,r,  g  =  p  =  a 


we  rewrite  eqns  (l)-(4)  as 

v  =  -  (5  -  Ip). 

P 

8  =  kg,t  +  .1(5;  +  a2), 
s  =  niu  —  As). 

O  =  -  y  <7  j  . 


(14) 

(15) 

(16) 

(17) 

(18) 
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Thus,  only  first-order  spatial  derivatives  of  the  unknowns  v,  9,  s,  a ,  u,  g ,  and  p  appear 
in  the  governing  equations.  Let  H'  denote  the  space  of  functions  defined  on  (0,1)  the 
square  of  whose  first-order  derivative  is  integrable  over  [0,1].  We  approximate  the 
unknown  functions  v,  9,  s,  etc.  by  linear  combination  of  the  finite  element  basis  func¬ 
tions  {<t>i(y),  i  =  1,2,  .  .  .  ,N}  in  an  /V-dimensionai  subspace  of  H'.  For  example, 

v(y,i)  =  v,(t)<t>,(  v).  (19) 

Throughout  the  article,  a  repeated  index  implies  summation  over  the  range  of  the  index. 
Using  Galerkin’s  method  (e.g.,  see  Becker  et  at.  [1981])  we  thus  reduce  eqns  (14) 
through  (18)  to  the  following  set  of  equations: 

M,j  «,  =  ~Q„  v, ,  M„g,  =  ~  Q„  9, ,  ,V/„  p,  =  -Q„  a,.  ( 20 ) 

M„v,  =  --  A,jS,  +  -  Q0p,,  M„e,  =  -kQ„g,  +.1  ,P„,  (21) 

P  P 

s t  ftM„ u,  R,,k,  bd^b,  — p-tQ,,it,  1 , R,//..  (22) 


where 


M, 


0,0,  dv  =  M„,  Q„  =  f  0,0,  ,dy, 
Jo  Jo 


(23) 


Qh  =  Qu  -  (0,0,) id, 


(24) 


R,]k  ~ 

P;  = 


0,0,0^  dy  =  R,k,  =  Rk,n 
0,0,(s:  +  o~)dy  =  P„. 


(25) 

(26) 


We  note  that  because  of  the  nonlinear  dependence  of  P„  and  .1  upon  s.  a.  v,  and  9.  the 
coupled  set  of  ordinary  differential  eqns  (20)-(22)  is  not  that  easy  to  integrate.  The 
matrices  .V/,,.  Q,r  0m  and  P„  have  been  evaluated  by  using  the  linear  basts  func¬ 
tions.  Also,  v,(t)  denotes  the  velocity  of  node  i  at  time  r. 

We  use  the  Crank-Nicolson  method  to  integrate  eqns  (20)-(22),  with  respect  to  time 
r.  In  it,  eqns  (20)-(22),  assumed  to  hold  at  time  (r  +  Jr/2),  are  used  to  predict  the  values 
of  i\  9.  5,  a,  g,  p ,  u.  and  1 4  at  time  (r  +  Jr)  from  a  knowledge  of  their  values  at  time 
r.  This  is  accomplished  by  approximating  9,(t+Ji/2)  by  (9,(1  +  J/)  -  9,(t) ),  Jr. 
9,(1  +  Jr/2)  by  (0,(r+Jr)  +  0,(r))/ 2,  and  so  forth  and  by  first  evaluating  the  non¬ 
linear  terms  on  the  right-hand  side  of  eqns  (20)-(22)  at  time  /.  The  resulting  system 
of  linear  algebraic  equations  is  solved  for  i’,(r  +  Jr),  etc.,  the  right-hand  side  in 
eqns  (20M22)  is  now  evaluated  at  time  (r  +  Jr/2),  and  the  system  of  equations  is  solved 

again  for  v,(t  +  Jr)  etc.  This  iterative  process  is  continued  until  at  each  nodal  point, 

J v  J 9  J.v  :  Jvf  ' 

—  +  —  +  —  +  —  +  Jo  +  J g  +  Jp  +  J  it  <(  (27) 

c  9  v  i 
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where  the  subscript  /  has  been  dropped  from  v,  and  elsewhere,  A  v  denotes  the  differ¬ 
ence  between  the  newly  found  value  of  v  and  that  used  to  compute  the  right-hand  side 
in  eqns  (20)-(22),  and  t  is  a  preassigned  small  number.  The  initial  conditions  (13)  were 
used  to  find  i',(0)  and  so  forth. 


IV.  COMPUTATION  AND  DISCUSSION  OF  RESULTS 


In  order  to  compute  numerical  results,  the  following  values  of  various  nondimen- 
sional  parameters  that  correspond  to  a  typical  hard  steel  and  the  average  applied  strain 
rate  •>,>  =  500  sec-’  were  chosen. 

p  =  3.928  x  10"',  A*  =  3.978  x  10"\  a  =  0.4973,  m  =  240.3, 

n  =  0.09,  i/'n  =  0.017,  b  =  5  x  106,  m  =  0.025. 

For  homogeneous  deformations  of  the  block,  the  peak  (marked  as  point  P  in  Fig.  I) 
in  the  shear  stress-shear  strain  curve  occurs  at  a  strain  of  0.093.  The  uniform  temper¬ 
ature  6„  =  0.1003  in  the  block  when  y  =  0.0692,  corresponding  to  the  point  /  in  Fig.  1, 
was  perturbed  by  adding  a  smooth  temperature  bump 

0(  v)  =  Jtf„<  1  -  v*  -  v,fi  )pe  ,  (28) 


-95- 


R.  C.  Batra 


and  the  resulting  initial-boundary-value  problem  was  solved  by  using  the  aforemen¬ 
tioned  method.  The  domain  was  divided  into  either  20  uniform  subdomains  (usually 
called  finite  elements)  or  20  nonuniform  subdomains,  with  nodes  at  0,  0.0025,  0.01, 
0.0225,  0.04,  0.0625,  0.09,  0.1225,  0.160,  0.2025,  0.2500,  0.3025,  0.360,  0.4225,  0.490, 
0.5625,  0.640,  0.7225,  0.81,  0.9025,  and  1.0.  The  two  temperature  perturbations  for 
J0O  =  0. 1,  p  —  9,  a  =  5  and  for  J0O  =  0. 1,  p  =  2,  a  =  1  are  shown  in  Fig.  2.  For  non¬ 
polar  materials,  Fig.  3  depicts  the  growth  of  the  central  plastic  strain  rate  in  time  for 
the  two  subdivisions  of  the  domain.  Numerical  experiments  with  different  values  of  Jt 
indicated  that  Jt  =  5  x  10 gave  accurate  results.  All  of  the  results  presented  herein 
are  for  this  value  of  Jt  and  t  =  0.01.  Unless  otherwise  noted,  the  nonuniform  grid  has 
been  used.  Wright  and  Batra  [1986]  gave  a  heuristic  reasoning  to  explain  that  the  pla- 
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TIME  ELAPSED  (/isec) 


Fig.  3.  Central  plastic  strain  rate w  elapsed  time. 


teau  was  formed  because  of  the  coarse  mesh.  The  plot  in  Fig.  2  clearly  proves  that. 
Also,  for  the  wider  perturbation  the  rate  of  growth  of  the  central  plastic  strain  rate  is 
quite  a  bit  slower.  Whereas  for  the  narrower  perturbation  the  shear  band  is  formed 
when  the  average  strain  in  the  block  corresponds  to  the  point  just  past  the  peak  in  the 
shear  stress-shear  strain  curve,  the  broader  perturbation  results  in  the  formation  of  the 
band  well  past  the  peak  in  the  average  shear  stress-average  shear  strain  curve,  which 
occurs  at  i  =  47.6  jisec  after  the  temperature  field  was  perturbed.  That  a  band  is  indeed 
formed  for  the  wider  perturbation  should  be  clear  from  the  plot  of  Fig.  4.  It  is  appar¬ 
ent  that  the  rapidly  deforming  region  progressively  narrows  down  to  the  one  close  to 
the  center.  The  width  of  the  rapidly  deforming  region  is  narrower  for  the  finer  mesh. 
Whether  we  have  obtained  the  final  band  width  or  not  is  not  quite  clear  yet.  since  we 
have  not  experimented  with  other  fine  meshes.  Figure  5  depicts  how  the  central  shear 
stress  changes  in  time.  It  decreases  at  a  very  slow  rate  and  essentially  follows  the  aver¬ 
age  shear  stress-shear  strain  curve.  It  stays  nearly  uniform  throughout  the  specimen. 
Numerical  instabilities  developed  for  values  of  time  i  a  little  beyond  the  one  up  to  which 
results  arc  plotted  here.  These  can  be  due  to  the  improper  size  of  the  time  increment 
and/or  to  the  mesh  size,  among  other  factors.  Efforts  are  now  under  way  to  overcome 
these  and  to  extend  calculations  for  longer  values  of  time. 

Figures  6  and  7  depict  the  effect  of  the  amplitude  J 0„  of  the  perturbation  upon  the 
growth  of  the  central  plastic  strain  rate  for  simple  and  dipolar  materials  with  /  =  0.01 . 
Obviously,  lor  larger  amplitudes  of  the  perturbation,  the  shear  band  is  formed  well 
before  the  peak  in  the  average  shear  stress-average  shear  strain  curve  is  reached.  A  com¬ 
parison  of  the  results  presented  in  Figs.  3,  6.  and  7  clearlv  brings  out  the  stiffening 
caused  by  the  inclusion  of  the  dipolar  stresses. 
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DISTANCE  y 

Fig.  4.  Plastic  strain  rate  distribution  in  a  Unite  block  of  material  at  various  limes  as  a  shear  hand  lorm 
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In  order  to  understand  the  interaction  among  shear  bands,  we  introduced  a  pertur¬ 
bation  in  the  temperature  centered  at  different  points  along  the  thickness  of  the  slab. 
Because  of  the  symmetry  of  d  about  the  horizontal  axis,  this  amounts  to  introducing 
two  aberrations  symmetrically  placed  about  the  center  line  y  =  0.  The  hypothesis  here 
is  that  if  the  resulting  nonhomogeneous  fields,  such  as  the  plastic  strain  rate,  temper¬ 
ature,  and  the  plastic  strain,  eventually  peak  out  at  the  center  of  the  slab,  then  the  two 
bands  have  coalesced;  otherwise  they  grow  independently.  Perturbations  in  the  temper¬ 
ature  centered  at  y  =  0.025,  y  =  0.05,  and  v  =  0. 1  but  =  0. 1  were  introduced  and 
the  ensuing  initial-boundary-value  problems  were  solved.  Figures  8  and  9  show  the  dis¬ 
tribution  of  the  plastic  strain  rate  through  the  thickness  of  the  slab  for  simple  and  dipo¬ 
lar  materials,  respectively.  For  simple  materials  (Fig.  8),  the  two  bands  initially  centered 
at  y  =  ±0.025  coalesce  but  the  ones  centered  at  ±  05  and  ±.l  grow  independently  of 
each  other.  The  rate  of  growth  of  the  peak  strain  rate  (not  the  strain  rate  at  y  =  0)  is 
considerably  less  as  compared  to  that  when  either  only  one  band  appears  at  y  -  0  or 
two  bands  initially  centered  at  y  =  ±.025  merge  and  grow  as  a  single  band.  For  dipo¬ 
lar  materials,  the  bands  coalesce  in  all  three  cases.  Recall  that  the  material  character¬ 
istic  length  is  ^  of  the  distance  between  the  bands  originating  at  y=  ±0.1. 

The  distribution  of  the  dipolar  stress  across  the  slab  is  plotted  in  Fig.  10.  Because  of 
the  boundary  conditions  a(0,r)  =  a(l.t)  =0  and  the  fact  that  the  dipolar  stress  lor  the 
elastic  problem  is  proportional  to  the  curvature',  the  maximum  value  of  o  cannot  occur 
at  the  center  of  the  slab.  Note  that  a  is  very  small  as  compared  to  1 .0,  whereas  v  is  gen- 
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erally  greater  than  or  equal  to  1.0.  At  a  point  where  yp  *  0,  dp  *  0,  the  dipolar  stress 
has  negligible  contribution  to  the  effective  stress  (s:  +  ct:)'  :  that  determines  whether 
A  >  0  or  .1  =  0  at  that  point. 


V.  CONCLUSION 

It  is  shown  that  the  constitutive  model  proposed  by  Wright  and  Batra  [1986]  does 
predict  adiabatic  shear  bands  in  a  block  undergoing  simple  shearing.  The  rate  at  which 
a  band  grows  depends  upon  the  strength  of  the  inhomogeneity,  herein  modeled  as  a  per¬ 
turbation  in  the  otherwise  uniform  temperature  within  the  block.  Perturbations  of  larger 
amplitude  result  in  the  formation  of  a  shear  band  even  before  the  peak  in  the  shear 
stress-shear  strain  curve  is  reached.  Also,  a  wider  perturbation  results  in  the  shear  strain 
localization  at  a  lower  rate  as  compared  to  the  narrow  perturbation,  both  being  of  the 
same  amplitude.  The  inclusion  of  dipolar  effects  results  in  a  very  stiffening  effect  in  the 
sense  that  the  formation  of  bands  is  delayed  considerably  as  compared  to  that  in  sim¬ 
ple  materials.  In  addition,  two  bands  that  would  grow  independently  in  a  simple  mate¬ 
rial  coalesce  when  dipolar  effects  are  included.  Of  course,  the  minimum  distance 
between  two  shear  bands  to  grow  independently  of  each  other  will  possibly  depend 
upon,  among  other  factors,  the  value  of  the  material  characteristic  length  /. 
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Abstract— The  thermomechanical  problem  involving  simple  shearing  of  a  finite  slab  made  of  an 
isotropic  and  viscoplastic  material  is  studied  with  the  objective  of  finding  the  effect  of  the  strain 
hardening  parameter,  strain-rate  hardening  parameters,  thermal  softening  coefficient  and  thermal 
conductivity  on  the  initiation  and  growth  of  adiabatic  shear  bands.  The  body  is  placed  in  a  hard 
loading  device,  i  e.  the  velocity  is  prescribed  at  its  top  and  bottom  surfaces.  A  shear  band  is  presumed 
to  have  formed  if  the  addition  of  a  perturbation  in  the  temperature  at  the  onset  of  plastic  flow 
results  in  a  localization  of  the  shear  strain.  The  critical  strain  at  which  the  band  begins  to  form  is 
found  and  its  dependence  on  various  material  parameters  is  investigated. 


INTRODUCTION 

In  1944  Zener  and  Hollomon[I]  recognized  the  destabilizing  effect  of  thermal  softening  in 
reducing  the  slope  of  the  stress-strain  curve  in  nearly  adiabatic  deformations.  They  pos¬ 
tulated  that  a  negative  slope  of  the  stress-strain  curve  implies  an  intrinsic  instability  of  the 
material.  Thus  the  strain  at  which  the  shear  strain  localization  may  initiate  corresponds  to 
the  peak  in  the  stress-strain  curve.  They  observed  32  /um  wide  shear  bands  in  a  steel  plate 
punched  by  a  standard  die  and  estimated  the  maximum  strain  in  the  band  to  be  100. 
Recht[2],  assuming  that  the  instability  occurs  at  the  peak  in  the  stress-strain  curve  and  this 
curve  is  independent  of  the  strain  rate,  derived  values  of  strain  rate  necessary  to  produce 
shear  strain  localization  and  compared  these  values  for  different  materials.  Slaker[3]  used 
the  same  instability  criterion  but  included  the  dependence  of  the  flow  stress  upon  strain 
rate  also.  Assuming  parabolic  strain  and  strain-rate  hardening  laws,  he  concluded  that 
important  material  parameters  are  the  specific  heat,  slope  of  the  temperature  dependence 
of  the  flow  stress,  and  parameters  indicating  the  strain  hardening  capacity  of  the  material. 
The  thermal  conductivity,  yield  strength  and  strain-rate  sensitivity  do  not  enter  in  as  a  first- 
order  effect. 

Instead  of  presuming  that  the  material  becomes  unstable  at  a  stress  maximum. 
Clifton[4]  and  Bai[5]  studied  the  growth  of  infinitesimal  periodic  perturbations  superimposed 
on  a  body  deformed  by  a  finite  amount  in  simple  shear.  Both  investigators  included  the 
effect  of  strain  hardening,  strain-rate  sensitivity,  thermal  softening  and  heat  conduction. 
Bai  also  included  the  effect  of  inertia  forces.  Bai's  instability  criterion  is  essentially  insensitive 
to  strain-rate  hardening  parameters  and  for  a  parabolic  type  strain  hardening  material  gives 
the  same  value  of  critical  strain  as  that  derived  by  Staker[3].  Burns(6]  used  a  dual  asymptotic 
expansion  to  account  for  the  time  dependence  of  the  homogeneous  solution  in  the  analysis 
of  the  growth  of  superimposed  infinitesimal  perturbations.  He  showed  that  the  growth  rate 
of  small  perturbations  is  controlled  by  the  ratio  of  the  slope  of  the  homogeneous  stress  vs 
strain  curve  to  the  rate  of  change  of  the  plastic  flow  stress  with  respect  to  the  strain  rate. 
However,  this  growth  rate  was  not  large  enough  for  Litonski's[7]  constitutive  relation  to 
account  for  the  experimental  observations  of  Costin  et  a/.[8].  Costin  et  ai.  observed  370- 
500  jim  wide  shear  bands  do';ng  dynamic  torsion  tests  involving  strain  rates  of  500  s' '  on 
short  specimens  of  1018  co.  oiled  steel.  Similar  observations  on  twelve  ductile  materials 
have  been  reported  by  J.-'m*,  .  et  a/.[9]  and  Lindholm  and  Johnson[10], 

Merzerfllj,  by  usinc  u  material  model  due  to  Bodner  and  Partom[12],  studied  the 
problem  of  twisting  of  a  thin  tubular  specimen  having  a  notch  in  its  periphery.  He  concluded 
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that  the  thermal  conductivity  played  a  key  role  in  determining  the  width  of  the  shear  band. 
On  the  other  hand  Wu  and  Freund[13],  by  using  a  different  material  model,  studied 
wave  propagation  in  an  infinite  medium  and  concluded  that  the  thermal  conductivity  has 
essentially  no  effect  on  the  width  of  a  shear  band. 

Clifton  et  at. [14]  used  both  numerical  and  perturbation  techniques  to  study  the 
initiation  and  growth  of  shear  bands.  They  used  a  power  law  model  and  reported  that  the 
thermal  conductivity  had  virtually  no  effect  on  their  results.  They  did  not  compute  the  band 
width  explicitly.  Also  the  rate  of  growth  of  the  nonuniformity  increased  strongly  with 
decreasing  strain-rate  sensitivity. 

Recently.  Wright  and  Batra[15-I7]and  Batra[18]  described  the  results  of  computations 
that  simulate  the  formation  of  a  shear  band  from  a  local  temperature  inhomogeneity  in 
simple  and  dipolar  materials.  The  constitutive  relation  was  derived  by  modifying  the  dipolar 
theory  of  Green  et  u/.[19]  to  include  rate  effects.  For  simple  materials  this  constitutive 
relation  reduces  to  one  that  is  very  similar  to  that  proposed  by  Litonski[7]  and  Lindholm 
and  JohnsonflO].  Whereas  Litonski,  and  Lindholm  and  Johnson  suggested  constitutive 
relations  valid  for  simple  shearing  deformation  of  a  ductile  material,  that  proposed  by 
Wright  and  Batra  is  easily  amenable  to  a  general  state  of  stress.  Wnght  and  Batra[  1 5—  1 7] 
and  Batra[18]  studied  the  simple  shearing  of  a  block  made  of  an  isotropic  and  viscoplastic 
material  and  perturbed  the  homogeneous  solution  by  adding  a  temperature  bump  just  prior 
to  the  occurrence  of  the  peak  stress.  The  stress  field  was  calculated  so  that  every  material 
point  was  on  the  yield  surface  corresponding  to  its  new  temperature  when  all  other  variables 
were  held  fixed.  The  full  set  of  non-linear  coupled  governing  equations  was  solved  numeri¬ 
cally.  In  Ref.  [18]  Batra  used  a  similar  method  to  study  the  interaction  among  shear  bands 
in  simple  and  dipolar  materials. 

Herein,  for  simple  materials,  we  study  the  effect  of  various  material  parameters  on  the 
initiation  and  growth  of  adiabatic  shear  bands  by  adding  the  temperature  perturbation  in 
the  configuration  when  the  body  just  starts  deforming  plastically.  Since  this  point  is  far 
from  the  peak  in  the  stress-strain  curve,  our  results  should  reflect  the  dependence  of  the 
rate  of  growth  of  the  perturbation  upon  the  material  parameter  being  changed.  We  should 
add  that  the  complete  set  of  coupled  non-linear  equations  are  solved  numerically  by  the 
Galerkin-Crank-Nicolson  method.  The  computed  results  show  that  for  the  problems 
studied  in  which  the  non-dimensional  thermal  length  varied  from  0  to  0.063.  the  thermal 
conductivity  has  no  noticeable  effect  on  the  strain  at  which  the  shear  strain  localization 
occurs.  In  all  but  one  of  the  cases  studied,  the  stress-strain  curve  had  a  peak  in  it  and  a 
narrow  region  near  the  center  eventually  deformed  very  rapidly  with  the  rest  of  the  material 
essentially  not  deforming  at  all.  For  the  exceptional  case  noted  earlier,  the  combination  of 
the  values  of  material  parameters  was  such  that  the  stress-strain  curve  had  no  peak  in  it. 
In  this  case  no  localization  of  deformation  occurred  even  when  the  amplitude  of  the  initial 
temperature  perturbation  was  increased  to  three  times  its  value  for  other  cases.  This 
seems  to  confirm  that  the  peaking  out  of  stress  is  a  necessary  condition  for  the  shearing 
deformation  to  localize. 


FORMULATION  OF  THE  PROBLEM 

We  study  the  simple  shearing  deformations  of  a  semi-infinite,  isotropic  and  viscoplastic 
body  bounded  by  the  planes  Y  =  ±H  and  consider  deformation  fields  of  the  type 

.f  =  .V  +  u(  Y.  t).  v  =  K.  r  =  Z.  H  =  0(  Y.  D.  ( I ) 

Thus,  with  respect  to  a  rectangular  Cartesian  set  of  axes,  (  t.  r.  r)  denote  the  current 
coordinates  of  a  material  point  that  occupied  the  place  (,V.  Y.  Z)  in  a  stress-free  reference 
configuration.  The  functions  w  and  I)  give,  respectively,  the  displacement  of  the  material 
point  in  the  t-direction  and  its  temperature  change  from  that  in  the  reference  configuration. 
In  the  absence  of  body  forces  and  external  sources  of  heat,  the  balance  of  linear  momentum 
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divided  into  20  subdomains  (finite  elements)  with  nodes  at  0,  0.0025,  0.0 1,  0.0225.  0.04. 
0.0625.  0.09.  0.1225.  0.160.  0.2025.  0.2500.  0.3025.  0.3600.  0.4225.  0.4900.  0.5625.  0.6400. 
0.7225.  0.8100.  0.9025.  1.0.  The  uniform  time  increment  At  =  5  x  10  6  was  used  in  the 
Crank-Nicoison  method. 

The  following  values  of  material  parameters  that  correspond  to  a  typical  hard  steel 
were  chosen : 

p  =  7860  kg  m  \  *  =  49.216  W  m  ' 1  C - '.  p  =  80  GPa. 

k„  =  333  MPa.  a  =  0.00552' C 1 .  c.  =  473  J  kg  1  C  '.  (20) 

m  =  0.025.  n  =  0.09,  t/t„  =  0.017,  b  =  I04  s. 

For  this  choice  of  parameters,  60  =  89.6“C.  Also  we  took  y0  =  500  s  1  and  H  =  2.58  mm. 

As  pointed  out  by  Wright  and  Batra(15],  implicit  in  eqns  (16)  are  two  relative  length 
scales,  namely  a  thermal  length  (kipcryaH2)1 2  and  viscous  length  (6/7/)  (Ko/p)1  '  The  effect 
of  a  change  in  these  as  well  as  in  the  values  of  parameters  describing  thermal  softening, 
work  hardening,  and  rate  hardening  of  the  material  is  studied.  The  viscous  length  was 
varied  by  altering  the  value  of  the  material  parameter  b  while  keeping  the  values  of  all  other 
parameters  unchanged.  Figure  1  shows  the  shear  stress,  the  temperature  change  and  the 
plastic  strain  rate  at  the  center  vs  the  time  elapsed.  Since  the  average  strain  rate  in  the 
specimen  is  kept  fixed  at  500  s'  '.  the  average  strain  plotted  as  abscissa  in  Fig.  1  and 
other  figures  is  related  to  the  elapsed  time  A /  in  seconds  by  yJV(l  =  500  A/.  We  note  that  the 
shear  stress  in  the  specimen  was  initially  nonuniform.  However,  after  a  brief  interval  during 
which  the  field  variables  essentially  regain  their  balance,  the  shear  stress  becomes  uniform 
throughout  the  slab  and  stays  uniform  up  to  the  time  results  plotted  here.  It  is  obvious 
from  the  stress-strain  curves  plotted  in  Fig.  1  that  with  the  increase  in  the  value  of  b  the 
peak  stress  increases  but  this  peak  occurs  at  a  lower  value  of  average  strain.  The  amount 
by  which  the  peak  moves  to  the  left  decreases  with  every  10  fold  increase  in  the  value  of  b 
suggesting  that  eventually  an  increase  in  the  value  of  b  will  not  affect  the  strain  at  which 
the  peak  stress  occurs.  The  central  plastic  strain  rate  increases  rather  slowly  first,  but  begins 
to  increase  rapidly  as  the  peak  in  the  stress-strain  curve  is  reached  and  eventually  increases 
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Fig  I  Effect  of  the  value  of  material  parameter  b  on  the  evolution  of  the  central  stress,  central 
plastic  strain  rate  and  the  central  temperature. 
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at  an  infinite  rate  when  the  average  strain  in  the  specimen  is  well  past  that  corresponding 
to  the  peak  stress.  This  point  is  indicated  by  F  on  the  stress-strain  curves.  The  temperature 
at  the  center  first  increases  linearly  but  soon  begins  to  increase  at  a  faster  rate  and  eventually 
grows  extremely  rapidly.  Since  the  shear  stress  decreases  once  the  peak  is  passed,  the  rate 
of  increase  of  temperature  is  not  as  fast  as  the  rate  of  increase  of  the  plastic  strain  rate.  For 
the  same  value  of  average  strain,  a  higher  value  of  b  results  in  a  higher  value  of  0.  We 
should  add  that  unlike  the  shear  stress,  the  temperature,  plastic  strain,  plastic  strain  rate, 
and  particle  velocity  fields  are  not  uniform  throughout  the  thickness  of  the  slab.  Figure  2 
depicts  the  cross-plots  of  the  plastic  strain  rate,  temperature  change  and  the  particle  velocity 
for  different  values  of  b.  The  values  of  time  elapsed  or  the  average  strain  for  different  curves 
are  not  quite  the  same.  This  is  due  .to  the  fact  that  computed  results  were  printed  for 
identical  values  of  t  but  the  explosive  growth  in  central  plastic  strain  rate  occurs  at  different 
times.  Except  for  a  narrow  region  near  the  center,  the  temperature  distribution  within  the 
specimen  is  unaffected  by  a  change  in  the  value  of  material  parameter  b.  As  expected  the 
deformation  has  localized,  and  the  majority  of  the  block  away  from  the  center  moves  as  a 
rigid  body. 

Following  Wright(22]  we  define  the  width  of  a  shear  band  as  twice  the  distance  of  the 
point  from  the  center  where  the  value  of  the  plastic  strain  rate  drops  to  one-tenth  of  its 
maximum  value  at  the  center.  We  note  that  during  the  development  of  the  shear  band,  the 
plastic  strain  rate  is  maximum  at  the  center.  At  a  time  later  than  the  one  when  a  dramatic 
nse  in  the  growth  rate  of  the  central  plastic  strain  rate  occurs,  the  value  of  the  plastic  strain 
rate  at  a  node  adjoining  the  center  of  slab  becomes  greater  than  that  at  the  center.  This  is 
due  to  the  grid  being  not  as  fine  as  is  perhaps  required  for  the  proper  resolution  of  the 
deformation  field.  Notwithstanding  this  shortcoming,  one  can  still  investigate  the  effect  of 
the  viscous  length  on  the  band  width.  For  all  five  values  of  the  viscous  length,  the  band 
width  came  out  to  be  116  /un.  In  each  case,  the  central  plastic  strain  rate  had  reached  a 
value  of  80  times  the  applied  average  strain  rate  of  500  s  '.  Numerical  experiments  with  a 
mesh  that  was  finer  near  the  center  yielded  the  same  value  of  the  band  width  and  gave 
a  maximum  value  of  the  plastic  strain  rate  at  the  center.  We  are  currently  developing 
a  computer  code  that  will  refine  the  mesh  adaptively  and  give  a  much  better  resolution  of 
the  deformation  fields  near  the  center  of  the  slab. 


Fig  ’.  Fifed  of  the  value  of  material  parameter  />on  I  he  distribution  ol  par  title  velocity,  temperature 
change  and  plastic  strain  rale  within  the  specimen 
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and  interna)  energy  may  be  written  as 


pti  =  S  1 

(2) 

(3) 

Here  p  is  the  mass  density  which  stays  constant  since  the  deformations  considered  are 
isochoric  and  the  effect  of  temperature  changes  on  the  mass  density  is  being  neglected,  .v  is 
the  shearing  stress  in  the  .v-direction  on  a  plane  v  =  constant,  e  is  the  specific  internal 
energy,  q  is  the  heat  flux,  a  superimposed  dot  indicates  material  time  derivative,  and  a 
comma  followed  by  r  signifies  partial  differentiation  with  respect  to  r.  We  assume  that  the 
shear  strain  has  an  additive  decomposition 

7  =  u  ,  =  7c  +  7r  (4) 

and  that  a  loading  or  yield  function  /  exists  such  that 

fis.ii.'}  p)  =  x  (5) 

where  /  is  a  monotonicallv  decreasing  function  in  yp  and  k  is  a  measure  of  the  work 
hardening  of  the  material.  The  criterion  for  elastic  and  plastic  loading  is 

f{s.  0.  0)  $  k.  elastic  (6) 

/(.v.  0.  0)  >  k.  plastic.  (7) 

In  the  latter  case,  the  sigh  of  yp  is  the  same  as  that  of  5  and  its  absolute  value  can  be  found 
uniquely  from  eqn  (5)  because  of  the  assumed  monotonicity  of  /.  The  reader  is  referred  to 


Wright  and  Batra[l5.  17]  for  further  discussion  of  the  yield  surface. 

Here  we  make  the  following  choice  of  constitutive  functions 

pe  =  ifiy; +pc,H  (8) 

q  =  -A0.  (9) 

K  =  K"(I  +  iA„)  (l0) 

xt/i  =  vyp  (ID 

/  =  ,v/(  I  —«<))(  I  +hyp)m  (12) 


where  /i  is  a  constant  shear  modulus,  c,  is  the  specific  heat  at  constant  volume,  k  is  the 
thermal  conductivity,  is  the  plastic  strain  in  a  slow  isothermal  reference  test  for  which 
the  stress-strain  curve  (neglecting  elastic  strains)  is  given  by  eqn  ( 10).  parameter  a  describes 
the  thermal  softening  of  the  material  and  material  parameters  h  and  m  give  its  strain-rate 
sensitivity.  From  eqn  (8)  it  follows  by  using  standard  thermodynamic  arguments  [17]  that 

.r  =  /iyc.  (13) 


Therefore 


R.  C  Batra 


Before  stating  the  initial  and  boundary  conditions  we  introduce  non-dimensional 
variables  (indicated  by  superimposed  bars  in  eqn  ( 1 5))  as  follows : 


£ 

II 

u  =  HU 

s  =  K0S. 

0 

^6 

III 

*  1 

II 

i . 

/  =  —  /, 

7  o 

7  =  V, 

K  =  KnK. 

4> 

H 

i: 

ii 

a  =  ad0. 

5  =  by0. 

p  =  pH2y:oiK0, 

£ 

=  ki  (pc,  (15) 

Here  y0  is  the  average  strain  rate  imposed  in  the  problem.  Dropping  the  overbars,  the 
complete  set  of  equations  in  non-dimensional  variables  may  be  written  as 


1 

v  =  s 
P 

6  =  kO  .,.,+r/p 
.f  =  i u(t\,  -yp) 


yp  =  max 


0. 


_ i _ 

(!+*/*«)■( I -«0)/ 


(16) 


Whereas  we  have  assumed  that  all  of  the  plastic  work  is  converted  into  heat,  some  authors 
(e.g.  Sulijoadikusumo  and  Dillon(20])  assume  that  only  about  90%  of  the  plastic  work  is 
transformed  into  heat.  Farren  and  Taylor[21]  found  that  in  tensile  experiments  on  steels, 
copper  and  aluminum,  the  heat  rise  represented  86.5.  90.5-92  and  95%.  respectively,  of  the 
plastic  work.  In  eqns  (16)  i>  is  the  particle  velocity  in  the  r-direction.  The  boundary 
conditions 


t(±l.  t)  =  ±1.  0,(±  I.  r)  =  0  (17) 

ensure  that  the  overall  applied  strain  rate  is  y„  and  the  deformations  are  adiabatic.  For  the 
initial  conditions  we  take 

r(  r,  0)  =  i\  iln y.  0)  =  0.  0(  y.  0)  =  0.1(1  —  i  -')"  e  ''  (18a) 

s(  .1-,  0)  =  (I  0»(l  +b)m.  (18b) 

Thus  the  initial  perturbation  in  the  temperature  is  introduced  when  the  material  just  starts 
deforming  plastically  and  the  initial  stress  distribution  is  adjusted  so  that  all  of  the  material 
points  are  on  their  corresponding  yield  surfaces. 

We  seek  solutions  of  eqns  ( 16M 18)  such  that  v  is  antisymmetric  in  r.  and  .v  and  0  are 
symmetric  in  y.  Thus  we  can  study  the  problem  on  the  domain  [0.  I]  and  replace  boundary 
conditions  (17)  by 

r(0,  0  =  0.  0  ,(0.  r)  =  0.  r(l./)=l.  0,(1.  M  =  0.  (19) 


RESPONSE  TO  PERTURBATIONS 

The  details  of  integrating  the  governing  equations  ( 16)  under  side  conditions  ( 18)  and 
( 19)  by  the  Galerkin-Crank-Nicolson  method  are  given  in  Ref.  [18).  The  domain  [0.  I  ]  was 
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Fig.  3.  Stress-strain  curves  and  the  evolution  of  the  central  plastic  strain  rate  and  temperature  for 

different  values  of  m. 

Another  way  to  alter  viscous  effects  in  the  materials  is  to  keep  h  fixed  but  change  the 
rate  sensitivity  exponent  m.  The  effect  of  this  change  in  the  value  of  m  is  depicted  in  Figs  3 
and  4.  Whereas  an  increase  in  the  value  of  m  from  0.005  to  0.015  resulted  in  higher  values 
of  the  critical  strains  at  which  the  peak  P  in  the  stress-strain  curve  occurs,  subsequent 
increases  m  the  value  of  m  hardly  changed  the  critical  strain.  However,  the  difference 
between  the  value  of  the  average  strains  corresponding  to  point  F  when  the  explosive 
growth  in  the  central  plastic  strain  rate  occurs  and  the  point  P  increases  with  m.  Our 
numerical  results  agree  with  the  analytical  results  of  Staker[3]  and  Clifton  ci  a/.[l4]  if  we 
restrict  ourselves  to  m  $  0.015.  a  value  typical  for  many  metals.  Clifton  i'i  al.[  14]  also  noted 


Fig.  4.  Distribution  of  the  temperature,  particle  velocity  and  plastic  strain  rate  in  the  specimen  lor 

different  values  of  m. 
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Fig.  5.  Stress-strain  curves  and  the  evolution  of  the  central  plastic  strain  rate  for  different  values 

of  n. 


that  even  though  the  critical  strain  may  be  insensitive  to  the  value  of  m,  subsequent  growth 
of  the  instability  will  be  affected  by  its  value.  At  early,  times,  the  central  temperature  is 
unaffected  by  m.  Curves  representing  the  evolution  of  the  central  temperature  more  or  less 
coincide  until  the  temperature  at  the  center  begins  a  sharp  rise.  This  behavior  differs  from 
that  when  the  values  of  b  were  increased.  Figure  4  depicts  the  distribution  of  the  plastic 
strain  rate,  temperature  and  particle  velocity  through  the  thickness  of  the  slab.  Note  that 
these  curves  do  not  correspond  to  the  time  when  the  shear  band  had  fully  developed  in  each 
case.  Achieving  that  goal  would  have  required  sorting  through  a  tremendous  amount  of 
computer  output. 

Figures  5  and  6  show  results  for  different  values  of  the  strain-hardening  exponent  n. 
An  increase  in  the  value  of  n  increases  the  strain  at  which  the  peak  P  in  the  stress-strain 
curve  occurs  and  also  the  value  of  the  strain  at  which  the  explosive  growth  in  the  central 
plastic  strain  rate  occurs.  Whereas  Bai[5]  and  Staker[3]  showed  that  the  critical  strain  is 
proportional  to  n.  our  computed  values  of  the  critical  strain  divided  by  n  gave  8.67.  1 1.4, 
12.14.  and  12.67  for  n  =  0.03,  0.05,  0.07,  and  0.09.  respectively.  They  assumed  a  parabolic 
type  <7p)"  hardening  rule  and  we  have  represented  this  effect  by  (1  Since  i/>„  «  I. 

our  criterion  will  increase  the  flow  stress  more  than  the  simple  parabolic  hardening  rule  will 
for  the  same  amount  of  plastic  deformation  and  identical  values  of  n.  In  Fig.  6  are  plotted 
the  central  temperature  vs  time  or  average  strain,  and  the  distribution  of  the  particle 
velocity,  temperature  and  plastic  strain  rate  within  the  specimen.  Again  these  plots  of 
quantities  vs  the  distance  from  the  center  are  at  different  stages  of  the  shear  band  develop¬ 
ment  for  different  values  of  n.  These  depict  that  qualitatively  there  is  no  change  in  the 
way  various  field  variables  evolve  in  the  specimen  as  the  strain-hardening  exponent  is 
increased. 

In  Fig.  7  are  plotted  the  shear  stress  at  the  center,  the  central  plastic  strain  rate  and 
the  central  temperatures  vs  the  average  strain  for  three  different  values  of  the  thermal 
softening  coefficient  a.  For  a  -  0.000552.  no  peak  in  the  stress-strain  curve  occurred  for 
strains  up  to  35%,  and  the  central  plastic  strain  rate  grew  at  a  snail's  pace,  the  central 
temperature  increased  linearly  with  average  strain  and  the  particle  velocity  (cf.  Fig.  8)  had 
a  linear  variation  through  the  thickness  of  the  specimen.  The  temperature  was  slightly 
higher  at  the  center  than  it  was  at  the  edge.  To  see  if  the  increase  in  the  amplitude  of  the 
temperature  perturbation  given  by  eqn  (18a I  would  result  in  shear-strain  localization,  two 
more  cases  with  the  central  amplitudes  equal  to  0.2  and  0.3  were  tried.  In  neither  case 
did  the  deformation  localize.  When  similar  numerical  experiments  were  conducted  with 
a  -  0.00552.  the  deformation  localized  near  the  center  at  values  of  strain  well  below  the 
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Fig.  ft  Temperature  change  vs  average  strain  and  the  distribution  of  the  temperature  change,  plastic 
strain  rate  and  particle  velocity  within  the  slab  for  different  values  of  n 


ones  at  which  the  peak  stress  occurred.  This  supports  the  viewpoint  that  the  existence  of 
the  peak  in  the  stress-strain  curve  is  a  necessary  condition  for  the  deformation  to  localize. 

The  expressions  for  the  critical  strain  derived  by  Staker[3]  and  Bai[5]  imply  that  the 
critical  strain  is  inversely  proportional  to  the  thermal  softening  coefficient  a.  When  we 
halved  the  value  of  a  from  0.005552  to  0.002776.  the  peak  in  the  stress-strain  curve  occurred 
at  an  average  strain  of  0.0551  instead  of  0.15.  Whereas  we  have  solved  a  complete  set  of 
equations  both  Staker  and  Bai  approximated  the  change  in  temperature  caused  by  plastic 
working.  As  in  the  other  cases  studied,  the  dramatic  growth  in  the  central  plastic  strain  rate 
occurs  at  a  value  of  strain  well  past  the  peak  in  the  stress-strain  curve.  The  width  of  the 
shear  band  is  hardly  affected  by  the  value  of  the  thermal  softening  coefficient  a  provided 
that  the  deformation  does  localize. 
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Ftg.  7  Effect  of  the  value  of  the  thermal  softening  coefficient  a  on  the  growth  of  the  central  stress, 
plastic  strain  rale  and  temperature 
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Fig.  8.  Variation  of  the  particle  velocity,  temperature  and  plastic  strain  rate  within  the  specimen  for 

different  values  of  a. 


Finally  we  investigated  the  effect  of  the  thermal  length  by  changing  the  value  of  the 
thermal  conductivity  k.  We  note  that  for  the  values  of  p.  y„  and  H  given  in  eqns  (20),  the 
non-dimensional  thermal  length  decreases  from  0.0631  to  0  when  the  values  of  thermal 
conductivity  k  are  changed  from  49.216  to  0.  For  values  of  k  in  this  range,  the  stress-strain 
curves,  plotted  in  Fig.  9.  up  to  the  peak  stress  are  unaffected.  However,  beyond  this  peak 
the  drop  in  the  stress  is  slightly  affected  by  the  value  of  the  thermal  conductivity.  Also  the 
strain  at  which  the  explosive  growth  in  the  central  plastic  strain  rate  occurs  increases  slightly 
with  an  increase  in  the  value  of  A.  (n  Fig.  10  are  plotted  the  variations  of  the  plastic  strain 
rate,  temperature  and  particle  velocity  in  the  specimen.  Again  qualitatively  there  is  no 
difference  in  these  plots  when  the  thermal  conductivity  is  varied. 

For  all  four  values  of  the  thermal  length  used,  the  band  width  came  out  to  be  116  pm. 
In  each  case,  the  central  plastic  strain  rate  had  reached  a  value  of  80  times  the  applied  strain 
rate.  These  results  are  in  agreement  with  those  of  Wu  and  Freund[13]  but  disagree  with  the 
conclusions  drawn  by  Merzer(l  I].  Note  that  the  constitutive  relation  and  the  method  of 


Fig.  4  Effect  of  the  thermal  conductivity  on  the  growth  of  the  central  shear  stress,  temperature  and 

the  plastic  strain  rate. 
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Fig  10  Distribution  of  the  plastic  strain  rate,  temperature  change  and  the  particle  velocity  in  the 
specimen  at  four  different  values  of  thermal  conductivity. 


studying  the  shear  band  development  used  herein  are  different  from  those  employed  by 
Merzer.  and  Wu  and  Freund. 


CONCLUSIONS 

The  problem  of  shear  strain  localization  in  simple  viscoplastic  materials  has  been 
studied  by  using  a  constitutive  relation  proposed  earlier  by  Wright  and  Batra[l5].  A  shear 
band  is  presumed  to  have  formed  if  the  addition  of  a  temperature  bump  at  the  center  of  the 
specimen  in  the  configuration  in  which  it  just  deforms  plastically  results  in  the  localization  of 
the  deformation.  The  complete  set  of  coupled  non-linear  equations  is  integrated  numerically 
by  using  the  Galerkin-Crank-Nicolson  method. 

These  numerical  experiments  reinforce  the  belief  that  the  existence  of  a  peak  in  the 
stress-strain  curve  is  a  necessary  condition  for  the  shear  band  to  develop.  The  non-dimen¬ 
sional  thermal  length  with  values  between  0  and  0.063  has  little  effect  on  the  critical  strain 
at  which  the  peak  in  the  stress-strain  curve  occurs,  the  strain  at  which  the  plastic  strain 
rate  at  the  center  begins  an  extremely  rapid  rise  and  the  width  of  the  band.  The  critical 
strain  did  not  come  out  to  be  inversely  proportional  to  the  thermal  softening  parameter  as 
has  been  approximated  in  some  theoretical  studies.  Even  though  values  of  the  rate-hardening 
exponent  m  greater  than  0.015  did  not  increase  the  value  of  the  critical  strain,  the  strains 
at  which  dramatic  growth  in  the  band  development  occurs  did  increase  with  m.  Thus  from 
a  practical  viewpoint,  higher  values  of  the  strain-rate  hardening  exponent  m,  the  viscous 
length  and  the  strain  hardening  exponent  n  would  delay  the  development  of  a  shear  band. 

Finally  we  note  that  conclusions  drawn  herein  are  strictly  applicable  to  the  constitutive 
model  used. 
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